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|. INTRODUCTION 


Past selection work on quantitative characters, notably using 
Drosophila, has usually been made with equal rigour on both sexes. 
In the present account numerous lines, from the same F,, have been 
differentially selected with respect to the sexes within the same line. 
Since the effective response to selection depends upon the corres- 
pondence between genotype and phenotype, any sexual difference in 
the correlation between genotype and phenotype would result in 
selection acting with different speeds in the two sexes. 

The failure of the phenotype to mirror exactly the genotype in 
all individuals from a culture means that in many selections individuals 
are selected which are genetically lower or higher than expected. 
If a mid level is required, this inaccuracy becomes more apparent 
than if high or low selections are being considered. In the case of 
extreme selections the direction is definite and any inaccurate 
selections will not affect the next generation so markedly, for the 
top or bottom stratum of the population is already represented in 
the selected individuals. In a mid selection there may be slight 
inaccuracies in selection. For example the environment may affect 
one or other level preponderantly and thus offer for selection individuals 
which are mostly from one extreme of the population and will result 
in a gradual drift of the mean. This implies directional interaction 
of genotype and environment. In general, mid-lines may be expected 
to drift more than the extremes but to drift erratically without specific 
and continued direction. In the selections under review the best 
methods for mid level continuation are demonstrated. 

The fertility of selected lines is usually impaired by continued 
successful selection. Selection for particular genotypes will un- 
doubtedly disturb the normal good fertility to some extent by correlated 
response but as most selections also tend to inbreed the population 
these two factors must act concurrently. The difference in rigour 
among the twelve selections gives an opportunity to study the effects 
on fertility of selection, for inbreeding was constant in all twelve 
lines. 
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2. EXPERIMENTAL METHOD 


Selection was practised on the abdominal chaete of Drosophila melanogaster 
described by Mather, 1941. Two inbred lines, Oregon + (inbred for 250 generations) 
and Samarkand-+-(inbred for 150 generations) were crossed reciprocally and a 
large F, was made up from the two reciprocal F,’s. The variously selected lines 
were taken from the F,. The selections, apart from the F, when 80 flies of each 
sex were counted, were made on counts of 20 99 and 20 gg. Two cultures, one a 
reserve, were made from each line using 2 99 and 2 gg. The reserve was never 
used unless the principal culture failed, which was seldom. Variability, therefore, 
was not preserved to the same extent as in the method employed by Mather and 
Harrison (1949) and the response cannot, in consequence, be expected to be so 
great. 

The whole experiment was repeated later using the same inbred material and 
using the same technique in order to assess more accurately the effects of the different 
selective methods. So much seems to depend on the accuracy of genotype selection 
in the first few generations that for a proper assessment of the responses several 
replicates would have been preferred but, due to the prodigious labour this would 
have involved, only two were done. 

The twelve selected lines are listed below together with their symbols which 
will be used throughout this account :— 


29 selected 3d selected Symbol 
High High HH 
High Mid * HM 
High Low HL 
Mid High MH 
Mid Mid MM 
Mid Low ML 
Low High LH 
Low Mid LM 
Low Low LL 
Mid Alternating M (H and L) 

High and Low ft 
Alternating 
High and Low Mid (H and L) M 
Alternating Alternating (H and L) (H and L) 
High and Low High and Low 


* The mid selection is at the mean of the particular sex in the culture counted at that 
generation. 

t Selection alternates each generation, i.e. first generation of selection was high, second 
low, third high, etc. gg and 29 both alternated in the same direction at the same time. 


3. CULTURE VARIATION 


It has already been pointed out that inaccuracies in selection are 
likely to occur in a highly variable culture. This is due to develop- 
mental differences within a culture which affect the expression of 
chaeta number to some extent. A comparison of the variances of a 
group of F, means with the mean of the same F, variances shows, 
however, that the difference between cultures is relatively low, while 
the variation within a culture is, as expected, much higher. This 
shows that the variation due to unknown environmental causes is 
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due in the main to variation between individuals in a culture and 
not to variation between the cultures. A little variation between 
cultures does occur, as demonstrated by Mather, 1949. 

















TABLE 1 
Mean chaeta number of 11 F,’s Mean squares on 
Or. x Skd) corresponding groups 
(20 gg counted in each class) of the 20 g¢ 
I. 34°55 4°137 
2. 34°75 6°421 
3- 34°95 2°158 
4. 34°85 5°42) 
5+ 34°5 4°15 
6. 35°1 }Mean square 0-268 4°158 | Average i 
7 33°35 8-947 | “Tare 5°945 
8. 33°55 12°263 
9. 33°7 6-316 
10. 34°7 4°211 
II. 34°55 7°211 











Temperature and culture condition do not affect the expression 
of chaeta number very greatly but the unknown developmental 
differences can give considerable differences in the variances. Culture 3 
has a low variance of 2-158, but culture 8 gives an estimate of 12-263. 
The selective sieve will doubtless vary in its efficiency if the individuals 
within the cultures can vary so greatly. It is the development of the 
individuals within a culture that is, in some cases, masking genetic 
differences between individuals. Between lines, which is the com- 
parison in which we are interested, the quirks of development may 
be expected to cancel each other out and the comparison of culture 
means is therefore a reliable appraisal of the genetic conditions 
prevailing in the line. 


4. SELECTION LINES 


The complete data from the twelve selected lines together with 
the repeat experiment are presented in the accompanying graphs 
(figs. 1 to 4). These show the comparison of the duplicated selections. 


HH. The response in HH: (the first experiment) was the same until $7 (the 
seventh generation of selection) as in the line in the second experiment, hereinafter 
called HHe. But HHr rose slowly after S7 while HH2 failed to respond at all 
after $7. 

HM. These two lines responded very differently. HM fell slightly until S12 
and then stabilised itself, while HM2 rose steadily for the whole 17 generations. 
This comparison is an object lesson in the way in which two similarly selected 
lines follow different paths. This difference in response may be due either to the 
different genetic potentialities of the first parents chosen or to the difference in the 
genetic make-up of the mid males, which being most heterozygous, are perhaps 
most likely to vary in response according to the vagaries of recombination while 
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some might fail to respond owing to chance homozygosity for recombinations of 


mid effect. 


One other possibility may be mentioned, namely that chromosomal 


selection of one or other set of inbred chromosomes might have been selected before 
recombination took place in HM1 thus minimising the polygenic variability available 
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Fics. 1-4.—The response of the various selections. 


for release. 


Solid lines show the means of the sexes ; 


dotted lines the females and broken lines the males. 


The reciprocal selection, MH, also responded differently in the two 


experiments while all other lines were on the whole similar. 
HL. HL fella little until S3 and stabilised itself for the remaining 14 generations. 


HL2 on the other hand rose a little and then fell slightly. 


In these two opposed 


selections the pull between selected gg and 9? was about equal for apart from the 
first and perhaps the last few generations the lines did not respond. 

















REVERSAL OF A SEX CHARACTER BY SELECTION 157 


MH. It was mentioned in the reciprocal selection, HM, that the divergence 
of the two experiments is perhaps greatest in this type of selection. Whether this 
divergence is consequent upon the type of selection or merely fortuitous is unproved. 
MHI rose steadily until Sg and from there remained fairly stable. MH2, however, 
remained stable throughout the seventeen generations of selection. 

MM. MM1z died at S4 and had to be remade from M (H and L) 29 and 
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(H and L) M @¢ obtained from S6 cultures. These two lines were thought 
to approximate as near as possible to the MM selection. The mean of MMg2 and 
the repeat of MM: drift down steadily though slowly for the run of the experiment. 
It would seem that for a stable mid line the MM type of selection is not the most 
satisfactory. 

ML. ML responded in a low direction throughout the selection. Here the 
power of the low selection on the g¢ as compared with the mid selection on the 
292 cannot be effectively established as ML2 died prematurely at S6. However, 

L2 
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it appears from the other lines in which a mid selection was used that whenever the 
other sex was selected in any way other than in a high direction the line never 
responded in a + direction. The mid selection, therefore, would seem incapable 
of giving a high boost to the selection. 

LH. In this opposed selection there is closer agreement between the replicates 
than in the reciprocal selection, HL. There is very little departure from the point 





CHAETAE 

















35 ett entettetet 35 fectratentntntehentetenttententras 
GENERATIONS 
Fic. 3. 


of commencement and it seems that the low female selection and the high male 
selection have balanced each other. No greater release or build-up of high or low 
polygenic complexes has therefore been favoured in this set. 

LM. The LM lines both responded in a low direction and were equal in response 
to the LL selected lines. LMr1 and LMg2 both showed a curious reversal in the 
number of chaete produced by gg and 99. This phenomenon will be dealt with 
in greater detail later in the account. 

LL. Both lines responded similarly, LL2 responding a little more than LL. 
No increased response in a low direction occurred after S11 when both lines become 
stable. The low level reached by the LL lines was no more than that reached by 
many other selections under discussion. 





———— 
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M (H and L). These two lines responded in a similar fashion to the LL lines 
except that no more response occurred after $7. Here it is evident that in this 
material a low response is easier to obtain than a high. 

(Hand L) M. The first line of this selection remained at the mid level throughout 
the seventeen generations of selection with some fluctuations. (H and L) M2 
remained at a slightly lower level, except for a rise at S6 and S17. 


4s (H&L)(H&L) cll (H&L)(H&L)2 
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(H and L) (H and L). The first line remained about the mid level throughout 
while (H and L) (H and L)2 was a little higher towards the end of experiment 
following a rise at Sg. 


If we summarise the results of these various selections further 
we have those which responded and those which remained unchanged 
in both experiments, in the first experiment only and in the second 
experiment only. 
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Although it is possible to hazard;some conclusions from table 2 
it must be emphasised that they rest on two sets of data only. The 
overall responses are small and this is unfortunate for it does not 
allow an accurate assessment of line differences. If we consider the 
selections in the first experiment then we could say that the g¢ are 

















TABLE 2 
In expt. 1 and 2 In expt. 1 only In expt. 2 only 
= 
Good response ke MH | HM 
High HH | 
Slight response See dan HH 
High | (H and L) (H and L) 
Good response | LM, M (H and L) ML | 
Low LL 
Slight response 
Low MM HM, HL 
| 
Mid level . | LM, (H and L) M, | (H and L) (H and L) | HL, MH 
LH 
| 








less efficient in the production of a high chaeta number under selection 
but if we pool the data from the two experiments, as we must, then 
this sexual difference in the response to selection disappears. 

If the sexes do not differ in their ability to respond to selection 
which seems to be the case from these data, then in an opposed selection 
like LH or HL the relative ability to release plus or minus variability 
must be shown by the direction of response. Neither LH nor HL 
(in both experiments) can be said to show any striking response in 
either direction. The double alternating selection (H and L) (H and L) 
also remains near or at the mid level. These results indicate that 
in this material no more release of variability is found in one direction 
over the other and show that selection is equally effective in both 
sexes, 


5. FERTILITY OF THE LINES 


Fertility counts were made on all the lines in the first experiment. 
This was measured by taking the total number of flies emerging 
from the culture. The figures listed below are averages of flies per 
bottle of the lines over the whole seventeen generations. 


? . . - 280 ML - . . ° - 290 
HH . . + 255 (HandL)M =: ; - 210 
M (H and L) . * 230 MH : : : > 190 
LH - , ° - 9g§0 HM : , . - 180 
LM :- ° 230 HL - . ° 175 


(H and L) (Hand L) - 230 (MM not comparable due to line 
being remade so omitted) 
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The fertility of all the lines was good in both experiments and 
from the fertility figures obtained from expt. 1 it appears that the 
extremes of selection (i.e. HH and LL) are not deleteriously affected 
by selection any more than any of the other lines, in fact they head 
the list for good fertility. Any selection, therefore, is just as likely 
to upset the fertility balance as an extreme type of selection. The 
inbreeding of a selected stock may therefore have a considerable 
effect on the fertility apart from the fertility unbalance due to selection. 

A comparison of the present chaeta responses with those of Mather 
and Harrison (1949) shows the present ones to be much lower even 
though the original cross was the same. As the responses were small 
it may be that there was not such a great correlated effect on fertility. 
This difference in response between the two experiments may indicate 
a change in the genetic make-up of Samarkand and/or Oregon inbred 
since its use by Mather and Harrison some eleven years earlier. 


6. SEX DIFFERENCES IN CHAETA NUMBER 


Normally the females have more abdominal chaete on the ventral 
surfaces of 4th and 5th segments than the males. The sex difference 
in chaeta number varies according to the total number of chaeta. 
Fig. 5 illustrates this effect. The data used for fig. 5 are obtained 
from the many levels obtained in the selections of Mather and 
Harrison (1949), which originated from the same material that 
was again used in this experiment, namely Samarkand x Oregon. 
The levels of Mather and Harrison were all established by a LL or 
an HH type of selection. 

The relationship of sex difference to total chaeta number gives 
us a yard stick by which we can measure any aberrant difference. 
In fig. 6 we insert the regression line and limits obtained in fig. 5 
and superimpose the sex differences occurring in the selection lines 
under discussion. Many points lie outside the limits of fig. 5 which 
may be due to the small effects operating on the sexes consequent 
upon the various selection methods. The path of regression employed 
is too steep or too high in overall position for the present data, for 
the HH lines of the present experiment fall below lowest values in 
fig. 5 though the sex mean is well within range. However, it appears 
that the spread of the variously selected lines is much greater than 
expected from previous experience and that the LM lines are extremely 
aberrant. 

Sex difference has already been shown to vary with slight genetic 
alterations (Harrison and Mather, 1950), but never to the extent 
of that shown by the LM lines (fig. 3). The LM selections show a 
complete reversal, the females carrying fewer chaete than the males. 
This reversal, however, is due not to equal changes in both sexes 
but only to the females which carry fewer than normally. The males 
remained constant throughout both experiments. In LMz1 when the 
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females became lower than the males at $7 the emergence of the 
females became about one half that of the males. In LMga2 the 
emergence of the males and females was equal. Blackened larve 
and pupe present in LM cultures were presumably the females 
which had failed to emerge. 

If in the LMr line crossing over had occurred between the X 
and Y chromosomes (Mather, 1944; Philip, 1935) then the females 
homozygous for the recombinant chromosomes may have been lethal 
while the heterozygous recombinant females showed the reduction 


J 
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wn 


SEX DIFFERENCE 
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“a 








-5 
a 2a 28 32 6°36 40 44 48 52 56 60 64 68 
CHAETA NO. (MEAN OF SEXES) 


Fic, 5.—The relation between high and low chaeta number and the sex diflerence. (Data 
obtained from “‘ Manifold Effect and Selection,” K. Mather and B. J. Harrison, 1949.) 


in chaeta number. The Y chromosome has been shown to be slightly 
active but due to its comparative inertness the recombinant would 
be expected to be less active polygenically than the normal X. The 
low female selection in these LM lines would favour the new 
recombinant X chromosomes while the hemizygous males with a 
recombinant X were not so markedly affected, the mid selection 
on the males not being strong enough to pick out recombinant Y 
chromosomes. That the phenomenon occurred again in the LMa2 
line may be said to be possibly fortuitous for the type of selection 
favouring recombinant females existed in the LH lines but no sex 
differences of the LM type occurred in LH: but in LHa, towards 
the end of the experiment the sex reversal began to appear. In LHa, 
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however, the reversal was due to the males increasing in chaeta 
number while the females remained at the same level. It appears 
to be a possibility then that the LHe selection was favouring 
recombinant Y chromosomes while the selection had failed to pick 
out low chaeta producing recombinant X chromosomes. Had this 
occurred the sex reversal would doubtless have been considerable. 
One other explanation for the alternation of the sex differences 
remains to be mentioned. It is possible that due to the sexually 
opposed selections practised here genes have been selected which 
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Fic. 6.—The sex differences obtained from the present selection lines plotted on the path of 
regression obtained in Fig. 5. Letters and plus sign indicate results from Experiment 1 
while plain key letters indicate Experiment 2. The positions are obtained by pooling 
results of the last five generations of selection (i.e. S :3-S 17). 


manifest themselves differently in the two sexes. The fact that they 
did not show themselves in all selections favouring such a situation 
must similarly be attributed to inaccuracies of selection. The poor 
emergence of females in LM1 may be due in this case to a correlation 
of chaeta genes, sex limited in expression, and genes for poor viability. 
Where the correlation did not exist, or was not so complete, females 
emerged and exhibited the sex reversal in chaeta number. 

The sex differences vary materially between many of the lines. 
In (H and L) (H and L)1 the difference is large throughout the 
experiment, while in (H and L) (H and L)2 the difference diminishes 
after S8. HL1 also shows a large difference while the rest of the 
lines, apart from the LMs are either less constant or the difference 
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is small. It is apparent that sex differences can be altered by selection 
but the mechanism by which they are obtained can only be conjectured 
here. 

The similarity of response of means in reciprocal selections shows 
that selection is equally effective in both sexes. The changes in sex 
differences, however, show that it is possible to get sex-limited responses 
to selection. This is presumably a way sexual selection in nature 
produces changes in secondary sex characters. 


7. SUMMARY 


1. Twelve different methods of selection have been tried on the 
abdominal chaete in Drosophila melanogaster. The sexes were subjected 
to opposed, and to partially opposed, selections, e.g. 92 low and 
33 high ; 29 mid level and 3¢ high. 

2. Opposed selections on the two sexes resulted in little or no 
response to selection. The best method for a mid level continuation 
of a stock is considered. 

3. The fertilities of the lines tested did not indicate a greater 
decrease in fertility with a more rigorous method of selection owing 
perhaps to the total smallness of the responses. 

4. The difference in behaviour of the duplicated selections shows 
the caution needed in interpreting the results of an isolated selection 
when comparisons of treatments are being made. 

5. Sex differences in chaeta number vary considerably between 
lines and in some cases seem to be correlated with the type of selection. 
The relation of sex differences in chaeta number and the mean chaeta 
number of the sexes is disturbed by the different methods of selection 
when compared with an arbitrary standard. 

6. Normally, 22 have more abdominal chaete than the gg, but 
in two selections the difference has actually been reversed. Thus a 
sex reversal in a secondary sex character has been obtained by 
selection. These reversals ‘can equally well be explained by assuming 
crossing-over between the X and Y chromosomes, or selection for 
sex-limited genes. 
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|. INTRODUCTION 


WHEN an individual, homozygous for one or other of two alleles 
of a gene for which the population is heterogenic, mates with another 
bearing to it a relationship closer than the average of the population, 
the offspring will in general include a greater proportion of such 
homozygotes than would the offspring of a mating of average relation- 
ship. Heterozygotes, on the other hand, give half heterozygotes and 
half homozygotes among their offspring, no matter what their mates 
may be. Thus inbreeding leads to a fall in the proportion of hetero- 
zygous offspring of homozygotes, without any change in the proportion 
of homozygous offspring of heterozygotes. The proportion of homo- 
zygotes will therefore increase in the population and, as it increases, 
we should expect the process to continue and the chance of any 
offspring of a homozygote being itself homozygous also to increase. 
So, mutation apart, the outcome of sufficiently close inbreeding will 
in the long run be complete homozygosis. The same result will, 
of course, follow where more than two alleles are present in the 
population. 

The approach to homozygosis depends on the number of homo- 
zygotes produced by heterozygous parents exceeding any loss of 
homozygotes brought about by crossing or indeed by any other 
means. Of other means, mutation will in general be unimportant. 
It may result in the final state falling slightly short of complete homo- 
zygosis ; but the proportion of heterozygosis maintained by mutation 
pressure must be small unless the gene has an unusually high mutation 
rate. It is worth noting that any mutation will tend to cause hetero- 
zygosis in inbred populations. This will be as true where both alleles 
are equally mutable as where the mutation rate of one exceeds that 
of the other. The effect of a difference in the mutation rates will be 
seen in a drift towards the preponderance of one allele in the population, 
a preponderance depending finally on the relative mutation rates 
of the alleles and not on their initial frequencies. With equality of 
mutation rates, on the other hand, the initial gene frequencies will 
be maintained apart from random variations. The change under 
unequal mutation rates must in any case be very slow, detectable 
only over very large numbers of generations and quite negligible in 
ordinary inbreeding experiments. 
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Mutation by itself may therefore be neglected in considering the 
outcome of inbreeding, except perhaps in special cases. The effects 
of selection, on the other hand, can be far from negligible. Where 
one allele is at an advantage, but heterozygotes as such are not, 
homozygosis will ensue at a speed depending on the system of inbreed- 
ing. The homozygotes for the favoured allele will then finally oust 
those for the less favoured at a speed depending on the selection 
pressure, and usually below that of the progress of inbreeding. With 
the heterozygotes at an advantage over both homozygotes, however, 
the whole course and even the outcome of the inbreeding will be 
altered, and altered profoundly if the selective differential is high 
enough. For, if the loss of homozygotes by selection is not less than 
their accretion from heterozygous parents, the proportion of homo- 
zygotes in the population will not increase in the population. The 
progress of the inbreeding will be prevented and heterozygotes will 
be maintained in the population. Inbreeding and selection will 
come to balance and annul one another’s effect. 


2. GENERAL CONSIDERATION 


In calculating the effect of selection favouring heterozygotes on 
the constitution of inbred populations we shall confine our attention 
to single genes. The results will, of course, apply to members of sets 
of genes affecting the same character provided they do not interact 
in producing their effects. 

The problem may be divided into four stages. 


(i) The complete description of the population when a particular 
mating system is in operation. 

(ii) The investigation of the way this description changes from 
one generation to the next in the absence of selection. 

(iii) The introduction of a parameter x(0o <x <1) as the average 
survival rate of the homozygotes compared with the 
heterozygotes. Initially x is assumed to be the same for 
all homozygotes. 

(iv) The description of the equilibrium population attained for 
various values of x. 


Six mating systems will be considered, viz. selfing, sib mating, 
parent-offspring mating, half-sib mating, double first-cousin mating 
and mating at random. Random mating, though not an inbreeding 
system, is included for the purpose of comparison with the others. A 
system of mixed selfing and random mating behaves in a similar way 
under selection against the homozygotes, but, as its mathematical 
treatment is different, it will be considered in a separate paper. 

It will be convenient to deal first with some general points and 
results. 

Consider the genotypes of individuals in a population containing 
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a certain gene with any number of alleles. As the mating system 
becomes more complex, so do the groups of individuals (mating 
units) by which the population is naturally classified. Under selfing 
we are concerned with the proportions of homozygotes and hetero- 
zygotes ; under sib mating with the various mating couples ; under 
double first-cousin mating with pairs of mating couples. In a large 
population and for the first five of the above systems, the frequencies 
of these groups in a generation are linearly related to those of the 
previous generation, so that the changes in frequency from generation 
to generation may be completely specified by a matrix, called the 
generation matrix of the system (Fisher, 1949). 

The parameter x is introduced by multiplying the rows of the 
generation matrix by suitable powers of x, and dividing the columns 
by suitable compensating linear forms in x. The maximum latent 
root of the resulting matrix determines the final stable equilbrium 
population, but, as these latent roots are functions of x, different 
ones, determining differing types of equilibrium, may be the greatest 
for different values of x. 

In the absence of selection (x = 1), the final result of continuing 
with any of the five inbreeding systems is homozygosis ; the population 
splits into sub-populations of close relatives (families) each containing 
only one allele of the gene under consideration. Below some lower 
value of x the final stable state is a population containing some hetero- 
zygotes, each family containing two alleles of the gene. With more 
complex mating systems populations whose families contain three 
alleles may be stable at some still lower value of x. However, in 
any range of x in which there are n alleles in families in the stable 
state, families containing m<n alleles can exist, but only in unstable 
equilibrium. 

Stable equilibrium and unstable equilibrium have the following 
meanings. In the former, if extra alleles are introduced (through a 
mating irregularity or mutation) the family will in time split into 
sub-families containing the original number of alleles. In the latter, 
extra alleles up to n-m in number, once introduced, will tend to remain 
in the family which will in time adjust its frequencies to a new (stable 
or unstable) equilibrium. Both of these types of equilibrium are 
stable to the arbitrary introduction into a family of extra individuals 
carrying alleles already in the family—in time the same equilibrium 
will be re-established. 

Linkage has no effect on the equilibrium state of a gene, but it 
may seriously upset the rate of approach to equilibrium. Linkage 
effects are ignored in this paper. 


3. THE MATING SYSTEMS 


3.1. Selfing.—For any number of alleles a, 5, c, . . . of one gene, the 
frequencies p and g (p-++g=1) of all homozygotes and of all heterozygotes 
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become (p+4g, }¢) in one generation. The generation matrix is thus 
39, 39 g g 








Parents 
aa ab 
aa I \ 
Offspring 
ab . 4 
where aa denotes all the homozygotes aa, bb,...... and ad the hetero- 
oo ee Seca a It is assumed that the population 


is sufficiently large for random variations to be negligible. If all the 
heterozygotes have the same viability and if only x of any of the 
homozygotes survive for each one of the heterozygotes, the effects 
of reduction in viability operating between and within families with 
the same effect (e.g. competition between members of different 
progenies being as strong as between members of the same family), 
the matrix becomes 
: 
; s 


The latent roots are x and } and these determine equilibrium popula- 
tions (1, 0) and (x, 1-2x) * respectively. After many generations 
the population is that determined by the maximum latent root. 


For x $ homozygosity results. 

For x<} some heterozygosity always remains (if it was present 
in the original population. Otherwise an unstable homozygous state 
is possible). 

After n generations the frequencies (p,, ¢,) of homozygotes and 
heterozygotes are given by 


nl, 9% -n_ px 
pax" p tf) +2 I —2x 
Gn 2-"q 


If we start with a completely heterozygous population (p, g) =(0, 1), 
the frequency of heterozygotes remaining after n generations of selfing 


I—2x 
I —x —2nyntl 


=h= 





for x #4 


2 
or = —— for x = 
2-+n 


noe 





* While (p, g) is taken to be such that +g = 1, it is convenient at times to use (u, v) 
such that uw and v are only proportional to the population frequencies. This need cause 
no trouble since where necessary we can divide through by (u+v). 
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so that as n becomes very large 
ho for x4 
I —2x 


or—> —__—_ for x <4. 
I—x 


Fig. 1 depicts the approach to equilibrium starting from complete 
heterozygosity (n = 0) or from a random mating population containing 
two alleles (n = 1). Near x = } the rate of approach to equilibrium 
is very slow. 


1.0 





n=O 
9 


NO 


30 


12) | 2 3 4 5 6 7 8 9 I. 
xX 


Fic. 1.—Approach to equilibrium under selfing. 








OO MK W 


h = proportion of heterozygotes in the population. 


x = relative survival rate of the homozygotes, the survival rate of the heterozygotes 
being taken as 1. 

n = number of generations of selfing starting from an entirely heterozygous population 

(n = 0) or from a random mating population containing only two alleles 


(n = 1). 


If we put d= 100(1—x), d is the percentage disadvantage of 
the homozygotes compared with the heterozygotes. Column (i) in 
table 1 contains, for each of the five systems, the value of d above 
which heterozygosity is possible. These figures are larger than those 

M 
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given erroneously in a previous paper (Mather and Hayman, 1952) 
on the same topic. 





TABLE 1 
i ii iii 
Selfing (S) é ‘ , : 50°0 66-7 500 
Sib mating (Sb) . ‘ ; ; 23°7 44°6 34°0 
Parent-offspring mating (PO1) . ; 19°8 39°9 15'2 
moe) ls : 25°0 47°2 34°6 
Half-sib ' mating (HSb) : ‘ 18:8 as 30°9 
Double first-cousin mating (DFC) ; 14°5 ve 23°2 




















The figures in this table are percentage selective disadvantages of homozygotes as 
compared with heterozygotes. 

(i) Shows the maximum disadvantage under which the inbreeding system will 
eventually produce complete homozygosity. 

(ii) Shows the disadvantage at which the population will eventually become 50 per 
cent. homozygous (see 4.3). 

(iii) Shows the maximum disadvantage of one homozygote for which the population 
will become homogenic when the other homozygote is lethal, only two alleles 
being present initially (see 4.3). 


3.2. Stb mating.—Here the population is appropriately classified 
into mating couples. Up to four alleles may occur in such a mating 
unit of which seven distinct types are therefore possible. The genera- 
tion matrix, allowing for selection, is 


aaX aa aa xX ab ab x ab aa x bb ab X ac aa Xx be ab x cd 











aaXaa . . x 4x2 om : 4x2 

aaxab . : : . x ; We ‘ x 

abxab . "i , $ 4 I } 3 } 
aaxbb . ‘ ‘ ; 4x7 

abxac . AMON 3 $ $ 
aaxbe . : : , ; , 4x 

aes | . . | ‘cae ia a Fis 
Column divisoral I ; x+1 x+1 I x+3 I I 


To show how this is derived, take the second column for example. 
Without selection the progeny of aaxab will give rise to the three 
types of crosses aaXxaa, aaXxab, abxab in the frequencies }: 4: 4. 
If only x of the homozygotes survive, these frequencies become pro- 
portional to }x?: 4x:}. Now if x of the progeny of aa x aa survive, 
3(x-+1) of the progeny of aa x ab will survive, so that }x*, 4x, } must 
each be divided by a factor, $(x-+1), to make their total equal to 
4(x+1). When x =1 this matrix reduces to that given by Fisher 
1949). 
This matrix provides a good illustration of some general theory. 
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A mating unit containing n alleles (an n-unit) can only give rise to 
units containing m<n alleles. Hence the generation matrix contains 
sub-matrices (delimited by dotted lines) along the leading diagonal 
below which only zeros occur. The characteristic equation of the 
generation matrix is thus the product of the characteristic equations 
of these sub-matrices. Let A,; be the ith latent root of the nth sub- 
matrix in order down the diagonal (under selfing A,, = x, As; = 4), 
A,, the greatest of the A,,; and V,;, V,, the corresponding latent vectors. 
Then the elements of V,; corresponding to m-units (m>n) are 
evidently zero. 

Now two standard results, (i) that the maximum latent root of a 
matrix with non-negative elements is real and lies between the 
greatest and least of the column (or the row) totals (and is actually 
unity if each column totals unity), and (ii) that the elements of the 
corresponding latent vector are real and of the same sign, may be 
applied to any sub-matrix of the generation matrix. Hence the 
elements of V,, corresponding to n-units have the same sign, while 
if A,>A,, (all m<n) all the elements of V,, have the same sign. Any 
such A, will therefore define a real population equilibrium. If 
further A,>A,, (all m+ n) then the equilibrium is stable; but if 
A,<A,, (some m>n) the equilibrium is unstable. 

It only remains to determine the regions in the range of x in which 
each of the A,, becomes the maximum for the whole matrix. Under 
selfing, sib, parent-offspring, and (it is presumed) half-sib and double 
first-cousin mating, A, is the greatest near x = 1, A, is the greatest 
for some lower values of x, A, at some still lower values, and so on as 
far as the complexity of the mating system allows. This is not true 
in general, because for an inbreeding system weaker than first-cousin 
mating full homozygosis is never attained, 7.e. A, is never the only 
maximum latent root of the generation matrix. However, there 
may still be a sequential appearance of X,, Aj.1, Ansa - - + as the 
maximum latent root as x decreases from 1 to o. In the first three 
of the systems under consideration only families containing either 
1 or 2 alleles are possible, depending on the value of x. For half-sib 
and double first-cousin mating, we have not determined whether 
families containing more than two alleles can exist or not. 

To return to sib-mating, the characteristic equation is 


(A —x) {4.(%+1)2A3 —2(x-+1)(2x+1)A? —x?(x+1)A+23} 

{(2A —1)(2(3-+%)A+) }(4A—1) = 0 
Hence A, =x and V,=(13;0;03;0;0;0;0) corresponding to 
homozygosis. Let A,;, Ag. and A; be the roots of the cubic. Then 
for any A,, and the corresponding value of x 

2 
Var = (Er +aa-ta} 5 an(r teas aC ta){(rta)a—apr ; 
x*{(1+x)A—x}3 03; 0; o) 
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As x varies from 0 to 1, each Ay; varies monotonically over the range 
given in the accompanying table. 








x | Agi | Aaa Aes 
oO 4 o fe) 
I H1+V5) t (1-5) 











Evidently A,,; = A, 


Asi _ 3 
Aes = headed so that A,, = A 
32 2(3 +x) 31 3 
vg — } 


Both A, and A, are <$<A, for all x in o to 1 so that we are only 
interested in the relative magnitudes of A, and A,. At x = 1, Ay>Ag 
and at x = 0, A,>A,. The changeover occurs when A = x in the 
cubic, which then becomes 4x?+3x?—2x—2—=0 of which the 











1.0 
n=O 
9 
8 
7 
6 
h -5 2 
4 3 
3 4 
5 
2 
' lO 
30 
fe) °° 5E oD ae aa es -e e  e 


x 


Fic. 2.—Approach to equilibrium under sib mating. h, x and n as in fig. 1. 
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maximum root is x = 0-7633 (and therefore d = 23:7). Below this 
value of x two-allele families form the stable state ; above it homo- 
zygosis results. Fig. 2 depicts the approach to equilibrium. Near 
x = 0-76 the rate of approach is very slow ; even after 30 generations 
the population is still 12 per cent. heterozygous. 

3.3.  Parent-offspring mating.—Systems in which an individual is 
used twice, first as offspring and then as parent, will be discussed. 
This double use of an individual for two generations introduces 
certain complications : the survival rate, x, of the homozygotes com- 
pared with the heterozygotes may change over this period and the 
number of offspring may be insufficient to mate all the parents, or 
vice versa. We shall consider two cases, firstly where all parents 
are mated, and secondly where all of a given proportion of offspring 
are mated. The former needs a high birth-rate for low values of «x, 
while in the latter the proportion mated must be low if the birth rate 
is high. 

(i). All parents mated. Generally we may suppose that com- 
pared with the heterozygotes only x of the homozygotes survive 
to the first mating and » of these to the second mating. If y= 1 
the population becomes homozygous for all x. If» =x we have 
the following generation matrix where the parent is written first in 
the cross and we note that each mated couple must have at least one 
gene in common. 








| aaX aa aa x ab ab X aa ab x ab ab x ac 
aaXaa. | x ¢ x 
aaXab. ‘ : | . : x 
abxaa. . , ‘ x . x x 
abxab. ‘ 4 I . I I 
abxac. - : | P 3 2 
Cones divisors . I x+1 x+1 x+1 x+3 


The characteristic equation is 
(A—x) {A((x+1)22 —(x-+1)A—#2) }{(x-+3)A—2} = 0 


As for sib-mating, homozygosis ensues near x = 1, while for low 
values of x two-allele families form the stable state. The changeover 
occurs at the maximum root of 


x8+-2x2 —x—1 = 0 i.e. at x = 08019. 
(ii). All of a proportion of offspring mated. The generation 
matrix is 
x ° hx ° 
hx ; hx he 
$ : 4 + 
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The characteristic equation is 

(A —x) {A(4A? —2A —x) }(A—}) = 0 
and changeover from homozygous to mixed equilibrium occurs at 
x = 0°7500. 

3.4. Half-sib mating.—If a male taken at random is mated to all 
its half-sisters who are also half-sisters of one another, then the 
population may be classified into the groups of three which are the 
parents of each pair of half-sibs. The surplus males may be neglected. 
For simplicity only two alleles of the gene are considered in this 
case and even then there are 10 mating units. 

The critical value of x is 0-8122, a root of an equation of the tenth 
degree. 


1.0 


9 


R3 


R2 


2 S PO2 POl 
tl S\b 








©] | 2 3 4 5 6 BS 8 9 10 
x 
Fic. 3.—The relation of proportion of heterozygotes (h) to survival rate (x) at equilibrium 


under selfing (S), sib mating (Sb), two kinds of parent-offspring mating (PO1, PO2— 
see text), and random mating with two (R2) and three (R93) alleles. 


3.5. Double first-cousin mating.—Again considering only two alleles, 
there are now 13 mating units, being all possible types of arrangement 
of four grandparents. The critical value of x is 08552. 
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3.6. Random mating.—With / alleles, the frequency of heterozygotes 
in the population = h = oe and each allele occurs equally often 
# _ 
provided that x # 1. 
Fig. 3 shows the frequency of heterozygotes in the stable state 
for selfing (S), sib (Sb) mating, parent-offspring (PO1, PO2) mating 
and random mating with two (R2) or three (R3) alleles. 


4. TWO ALLELES WITH DIFFERENT VIABILITIES 


In the foregoing all homozygotes were assumed to have equal 
disadvantages. Restricting ourselves to two alleles (a and 5), we 
will now consider the case where the two homozygotes have unequal 
disadvantages, aa, bb, and ab having the respective viabilities x, » 
and 1 (x and y<1). It obviously now becomes necessary to treat 
the two homozygotes separately in the generation matrix. 

4.1. Selfing—The generation matrix is 





| aa bb ab 
ae - | x . he 
a d : y ty 
ab : ay . ; $ 


The latent roots and corresponding equilibrium populations are 
= % V = (1; 0; 0) 
= V = (03 1; 0) 

2 V = (a(1 —29) 5 p(t —2x) 5 2(1 —2x) (1 —29)) 


i 


When x> 4, y the population becomes homozygous aa. 
When »> 4, x the population becomes homozygous bd. 


When $>x, _y some heterozygosity may remain. In this case Hardy’s 
law is satisfied if 


(1 —2x)*(1 —ay)® = xy(1 —2y) (1 —22) 


1.e. 1 —2x—2y+2xy = 0. 

Fig. 4 shows the type of final stable population for various values 
of x and ». The areas marked A and B correspond to homozygous 
aa and bb respectively ; in C the population, though containing 
heterozygotes, has more homozygotes than a population mating at 
random without selection; in D it contains more heterozygotes ; 
on the boundary between C and D, the effects of inbreeding and 
selection exactly counterbalance. There is a smooth transition across 
all boundary lines except that between A and B. 
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4.2. Stb mating.—The generation matrix is 


aa X aa bb x bb aa xX ab bb x ab abx ab aa x bb 








aaXaa. | x , 4x2 3 4x? 
ie me ae ee er 

; | , x ; x 
bbxab. . : | : , ; y y , 
abxab. ; “a : : $ 4 I I 
aaX bb . ° ‘ | dxy 
Column divisors . | I I x+1 +1 x+y+2 I 


The characteristic equation is 


(A—x)(A—») f(A, * 9) = 0 


F(A, x,y) being a quartic in A. Putting A = «x in the quartic determines 
the boundary between A and C, and A = y the boundary between 
Band C. The Hardy law curve is given by an equation of the twelfth 
degree. 


4.3. Parent-offspring mating.—Fig. 4 gives the regions for each type 
of equilibrium both for the case when all parents are mated (POr) 
and when all of a proportion of the offspring are mated (PO2). The 
Hardy law curve has been dotted as it has been determined only at 
its endpoints and centre point. 

Column (ii) of table 1 contains, for each inbreeding system, the 
values of d (= 100(1—x)) where x = y cuts the Hardy law curve, 
i.e. at which the population is 50 per cent. heterozygous with equal 
proportions of each allele. These values may, of course, be obtained 
by solving h=} in the formule of section 3. 

Column (iii) of table 1 contains the value of d above which hetero- 
zygosity is possible when the other homozygote is lethal (y = 0). 


4.4. Half-sib and double first-cousin mating —These are too cumber- 
some for a full treatment, but the case when y=o has been 
investigated. The critical values of x are given in table 1. 


4.5. Random mating.—The frequencies of aa, 6b and ab are pro- 
portional to x(1—y)?, »(1—x)?, 2(1—x)(1—y) respectively. This 
agrees with the result given by Sturtevant and Mather (1938). 





Fic. 4.—Phase diagrams showing the types of equilibrium achieved for two alleles, the 
homozygotes of which have different viabilities, under various systems of inbreeding 
(S, Sb, PO1, POa, as in fig. 3, HSb = half-sib mating, DFC = double first-cousin 
mating). 
x = survival rate of one homozygote (aa) compared with the heterozygote (ad). 
y = survival rate of the other homozygote (bd). 

When the point (x, y) is in A the population becomes homozygous aa, in B homo- 

zygous bb. In C the population becomes less heterozygous than a population mating 
at random without selection, in D more heterozygous. 
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4.6. The following results which emerge from this section are 
additional to those mentioned in sections 2 and 3. 

(i) Strong selection against (even lethality of) one homozygote 
does not necessarily cause the corresponding allele to disappear from 
the population. Selection, not necessarily so strong against the other 
homozygote can hold the heterozygote in the population. 

(ii) On the common boundaries of the areas A, B and C of fig. 4 
two of the latent roots of the generation matrix are equal and at the 
common point of A, B and C three are equal. Near these boundaries 
the latent roots are nearly equal and two of them have a ratio equal 
to 1: 1+8 (where 8 is very small) so that the relative contributions 
of the two corresponding latent vectors after n generations will be 
in the ratio 1 : (1+8)" = 1-+n8, nearly, i.e. the latent vector corres- 
ponding to the smaller latent root disappears linearly and not 
exponentially. One example of this effect is in section 3 when, for 
x = y, it was noted that the disappearance of the heterozygote was 
very slow near the critical value of x. Another example was found 
in a more detailed investigation of mixed selfing and random mating, 
with general values of x and _», where the same phenomenon occurs 
near x = y = 1. This means that there is a critical region near the 
common boundaries of A, B and C in which, although the equilibrium 
is stable and fixed by the mating system and survival rate of the 
homozygotes, small disturbing influences may cause the population 
to drift widely from equilibrium. 


5. OTHER MODES OF SELECTION 


It has been assumed so far that there is no separation of the 
progeny of each mating, or at least, that this separation, if it exists, 
does not cause any difference between the relative survival rates of 
individuals in the same or in different families. This will be true 
of many plants in the wild or in populations of cultivated crops, but 
it may be only partly true of animals, depending on the stages in the 
animal’s life cycle at which selection acts, or of experimental families 
of both plants and animals. We shall consider briefly two other 
schemes, one in which the same number of offspring is selected at 
random from each mating of a generation to provide matings for the 
next generation, and another in which it is assumed that all homo- 
zygotes having some heterozygous sibs show a survival rate differing 
from that of homozygotes all of whose sibs are homozygotes. Only 
the first three mating systems will be discussed in this section. 

5.1. Constant number of progeny—The generation matrices of 
sections 3 and 4 will now contain unity instead of x and _» in their top 
left-hand corners. This is the maximum latent root of the matrix 
except at x = 0. Hence the population eventually becomes com- 
pletely homozygous for all values of x. 

5.2. Differential survival rates. Let z be the survival rate (compared 
with the heterozygotes) of homozygotes with only homozygous sibs, 
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and let x be the survival rate of the other homozygotes. The generation 
matrices are those of section 3 with the x in the top left-hand corner 
replaced by z. Then A, = z, and A, is, as before, the maximum latent 
root of the second sub-matrix down the leading diagonal. A graph 
of this maximum latent root against x shows, for any given x, the 
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Fic. 5.—Phase diagrams showing the types of equilibrium when survival rates of homo- 
zygotes differ according to their origin. 
Z = survival rate of homozygotes with only homozygous sibs. 
x = survival rate of the other homozygotes. 
When the point (x, z) is in Ho the population becomes homozygous, in M some 
heterozygotes remain. Along the broken line z = x. 


values which z must exceed for the population to become completely 
homozygous. Below this value some heterozygosity will always 
remain. 

Fic. 5 shows phase diagrams for x and z for the three mating systems. 
The diagrams for half-sib mating and double first-cousin mating are 
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similar to those for Sb and PO. In the areas marked Ho the values 
of x and z are such that the population will eventually become homo- 
zygous ; in M it will remain mixed homozygous and heterozygous. 
As we expect, for low values of x, a value of z considerably greater 
than x is necessary for complete homozygosity, while for values of x 
near 1, z needs to be less than x to produce a mixed population. In 
fact, if z is near 1 no amount of selection against the other homozygotes 
can stop the population from becoming homozygous sooner or later. 
The broken line has only been included to divide the diagrams into 
the two parts in which either x>z or x<z. The value of x at which 
this line meets the boundary between Ho and M is that obtained in 
section 3 and tabulated (as d = 100(1 —x)) in column (i) of table 1. 


6. PRACTICAL IMPLICATIONS 


Perhaps the most striking outcome of these calculations is the 
size of the disadvantage which can suffice to prevent the attainment 
of full homozygosity under any system of mating other than selfing. 
With sib mating, the closest system of inbreeding possible in many 
species, a disadvantage of only 24 per cent. in the homozygotes will 
serve to maintain heterozygotes in the population, and even at 
disadvantages lower than this the progress towards homozygosis will 
be very slow. 

Natural selection in favour of individuals heterozygous for structural 
differences of the chromosomes has already been shown to have a 
profound effect on the outcome of inbreeding in Campanula persicifolia 
(Darlington and La Cour, 1950) and Drosophila subobscura (Philip 
et al., 1944). Evidently differences associated with these structural 
alterations are of the order of magnitude in question. Few com- 
parative figures are available for uncomplicated genic differences, 
but disadvantages large enough to be effective seem to be not unlikely 
in most if not all species of naturally outbreeding plants and animals. 
The relative weakness and incapacity of inbred lines is proverbial in 
such species, and indeed it is often considered desirable to select 
for vigour during inbreeding programmes with domestic animals 
and cultivated plants. In so far as it is successful, such selection 
would, it would seem, be likely to achieve its object largely by 
delaying or even in some circumstances preventing progress of the 
inbreeding. 

The vitiating effects on the progress of inbreeding of the dis- 
advantage, which homozygotes may suffer, may of course be reduced 
by suitable management in experiment. Any step which reduces 
the intensity of competition between individuals even within families 
should assist towards this end. In plants, the avoidance of crowding 
in seed pans, in seedling boxes and in the field, will help to reduce 
the intensity of such competition. The individuals to be taken for 
breeding can also be chosen without reference to their vigour or 
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fertility so that their representatives in the next generation will be 
independent of their phenotypes. Even so, however, there must always 
remain the possibility of competition between pollen grains and 
tubes in the process of fertilisation, and there is the further possibility 
of competition between developing embryos, especially where many 
seeds are normally developed in the capsule. It is thus likely that 
although the effects of disadvantage on the attainment of homozygosity 
may be reduced, they will not be eliminated entirely. 

In the same way intra-uterine competition and post-natal com- 
petition during suckling will always be potential sources of disturbance 
in mammals, as will post-hatching competition in birds or larval 
competition in insects, except where each individual is raised virtually 
separately—a course which must in general be impracticable. It is 
not therefore surprising that Dziigiines (1950) finds evidence of slow 
progress of inbreeding in poultry. Progress towards homozygosis 
may be aided even more by ensuring that the offspring of homozygous. 
parents are not at a disadvantage with respect to those of heterozygotes, 
that is by ensuring that, so far as possible, all pairs of partners contribute 
equally to the next generation. This precaution may, however, only 
be partially successful since even in the absence of competition from 
heterozygotes the homozygotes may die out differentially by reason 
of reduced vigour and disease resistance. 

In so far as the fitness of homozygotes depends on polygenic 
systems the present findings must be applied with caution, since it 
is characteristic of such systems that new and superior balances can 
be achieved given sufficient time. This is, however, more of a long 
term than a short term possibility. The linked polygenic combinations 
will be expected to behave more or less as units, effective factors as 
they have been called (Mather, 1949), for short periods, their properties 
of change and readjustment becoming important only as the generations 
accumulate. The effects of disadvantage of homozygosis for polygenic 
combinations may be expected to be felt with full force in the average 
inbreeding experiment, even though this disadvantage can be mitigated 
in the long run by readjustments such as must have occurred during 
the rise of inbreeding species from outbreeding ancestors in nature 
(Darlington and Mather, 1949). 

In designing inbreeding programmes, especially in outbreeding 
species, the consequences of disadvantage of homozygotes must be 
faced. Two implications of such disadvantage are clear. First, the 
inbreeding should be 23 close as possible: there is nothing to be 
gained from slow inbreeding. If the lines produced by close mating 
are too poor to survive or to be of value, anything better which results 
from a looser mating system will be better only because it has stabilised 
at a higher level of heterozygosity. Such a line may be of value 
where an equilibrium at which lower homozygosity is combined 
with a reasonable fitness is preferable to high homozygosity accom- 
panied by low fitness, as may well be the case in crop and stock 
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improvement work. But let such a situation be recognised for what 
it is, and any delusion dispelled that the same level of homozygosity 
has been produced by the slow inbreeding, albeit at the sacrifice of 
time, as would have been achieved by a closer mating system. 

The second point concerns the interruption of a system of inbreed- 
ing. In the absence of selective disadvantage, the interruption of a 
system by the substitution of a looser mating for a generation or so 
will in general merely slow down the progress towards homozygosity 
(Fisher, 1949). When the homozygotes are at a disadvantage, 
however, an actual loss of homozygosity might occur. Such will be 
the case when the proportion of homozygotes, though lower than 
that at which equilibrium would be reached with the closer system of 
inbreeding, is higher than that characteristic of the looser mating 
type. Since these various equilibria are all stable, the proportion of 
homozygotes will fall back towards the equilibrium value characteristic 
of the looser system should it have become too high by reason of the 
closer systems applied earlier. The same will be true when the 
proportion has risen too high as a result of sampling variation, except, 
of course, where the drift has led to complete homozygosis in all the 
breeding lines. In such a case only mutation could lead to a reapproach 
to the normal equilibrium. 

Drift may have another effect. Where the two alternative homo- 
zygotes are equally, or nearly equally, fit, they will, apart from drift, 
tend to appear in the population with frequencies relative to one 
another which reflect the initial gene frequencies of the outbred 
ancestors rather than the relative fitness. The relative frequencies 
of the homozygotes would in such a case be liable to marked change 
by drift in small populations. Where, however, the two homozygotes 
differ sharply in fitness, their relative frequencies will be determined 
chiefly by this difference and selection will correct any changes which 
may arise from drift. 

One last point remains to be made. Agreement of the proportions 
of the two homozygotes and the heterozygote in a population with 
the Hardy law is apt to be taken as indicating that mating is at random 
within that population. Certainly, such an agreement will follow 
from random mating, other things being equal. We can now see, 
however, that the Hardy proportions may appear even with close 
inbreeding, provided the homozygotes are at the necessary dis- 
advantage, and this disadvantage, though high in a system of complete 
selfing, is not impossibly great with looser inbreeding (table 1 (ii)). 
Indeed disadvantages of a magnitude which experiment would 
suggest as not impossible in a natural population, may combine with 
a moderately inbreeding system to give proportions of homozygotes 
lower than expected from Hardy’s series. Balanced polymorphism 
may therefore accompany inbreeding. The extreme case of this is, 
of course, found in Oenothera where regular selfing is accompanied 
by permanent hybridity because homozygotes are either lethal or 
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are prevented from occurring by special kinds of gametic lethality 
and competition. Given a knowledge of the selective properties of 
the various genotypes, conclusions may be drawn about the system 
of mating prevailing in a population from the frequencies of the 
genotypes. But in the absence of such information, neither random 
mating nor any other system can be inferred with full confidence. 


7. SUMMARY 


The effect of selection against homozygotes on the progress of 
inbreeding in large populations is considered. At disadvantages 
above a critical value which varies with the system of inbreeding 
(see table 1), complete homozygosis is never obtained, the population 
continuing to carry heterozygotes at equilibrium. 

Diagrams are given to show the progress of inbreeding and the 
equilibrium reached for the case where any number of alleles are 
present, but all homozygotes have the same relative disadvantage, 
and for the case of two alleles, the homozygotes for which may have 
different disadvantages. 

Consideration is extended to situations in which the homozygous 
offspring of homozygous parents are at a lesser disadvantage than 
homozygotes at least one of whose parents was heterozygous. Where 
all parents contribute equally to the next generation, the final state 
of the population under the inbreeding systems discussed will be one 
of homozygosity. 

The practical implications of these findings are discussed and 
it is emphasised that with homozygotes at a disadvantage experi- 
mental inbreeding need not lead to homozygosis ; that inbreeding 
should be as close as possible if homozygosis is being sought ; that 
interruption of the inbreeding may cause loss of homozygosis ; and 
that agreement with the Hardy law need not imply random mating. 
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1. INTRODUCTION 


THis investigation arose as part of a_ general investigation 
(Hayman and Mather, 1953) of behaviour under inbreeding when 
homozygotes are at a selective disadvantage compared with hetero- 
zygotes. The inbreeding systems there discussed were handled by 
matrix methods not applicable to mixed selfing and random mating, 
for which an elementary algebraic treatment is necessary. This 
mating system occurs in plants which have an imperfect selfing 
mechanism combined with either widespread pollen dispersal or 
seed mixing by an external agency (as in cultivated crops). 

One gene with only two alleles, a and 4, is considered, so that a 
population containing this gene may be classified into the three 
genotypes aa, bb, and ab with frequencies p, g and r. Any such 
population may be represented by a point in a trilinear diagram, 
based on an equilateral reference triangle, and such that the lengths 
of the perpendiculars from the “‘ population point ” to the sides of the 
triangle are proportional to p, g and r. We shall investigate how a 
large population changes from generation to generation under certain 
conditions, and represent these changes by curves in a trilinear 
diagram. Further, we shall obtain phase diagrams, similar to fig. 4 
of the paper quoted, showing broadly the equilibrium populations 
for various magnitudes of the selective disadvantages of the homo- 
zygotes aa and bb. The populations are assumed to be large enough 
for random fluctuations in gene frequencies to be unimportant. 


2. EQUILIBRIUM 


Let p, gq and r (p+g-+r = 1) be the frequencies of aa, bb and ab 
respectively and let x, y and 1 (x and _y <1) be their respective viabilities. 
Let s and ¢ (s+¢ = 1) be the frequencies of selfing and random mating. 

Starting with a population (f», go, 7) and no selection, the 
frequencies in the next generation under selfing are 


{Pot+4ro, Jot+4o 37o} 


while under random mating they are 


{(Pot+310)?s (Go+310)*, 2(Po+3%o) (Go+3"o)} 
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Under the mixed system and in the presence of selection they become 
(P1» Ya» 1) where 
Pi xx{s(Potdro) +t(potdro)*} = x{Pot+d (sro —lo)} 
11 X{S(Jotdro) +H(Got+470)*} = Go td (so leo) } » (1) 
1X $579 +2t( Potro) (Yo+so) = To—4(57o —leo) 


and e = 4pq-—r 


s = 1.00 s=.75 
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Fic. 1.—Phase diagram showing types of equilibrium with various proportions of selfing 
(s) and random mating (¢ = 1—s). 
x = survival rate of one homozygote (aa) relative to that of the heterozygote (ad). 
y = relative survival rate of the other homozygote (bd). 

When the point (x, y) is in A, the population becomes homozygous aa, in B, homo- 
zygous bb. In C, the population becomes less heterozygous than a population mating 
at random without selection, in D, more heterozygous. 

The points a, 6, c, d and e in the diagram for s = 0-50 are those referred to in 4.2 
and fig. 4. 
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0 Fo "o 
x{1-+ a (sr —te)} = y{1+ re (sr —te)} = 1 — = (sr—te) . (2) 


These equations determine p, g and r in terms of x, y and s. 

Fig. 1 shows the type of equilibrium population for various values 
of x, y and s. The areas marked A and B correspond to homozygons 
aa and 66 respectively ; in C the population, though containing 
heterozygotes, has more homozygotes than a population mating 
purely at random without selection ; in D it contains more hetero- 
zygotes ; on the boundary between C and D, the effects of inbreeding 
and selection exactly counterbalance each other. (The small letters 
in the diagram for s = } will be referred to in 4.2.) 

The inner boundaries of A and B are obtained by making r tend 
to zero in the equilibrium equations (2). The boundary between 
A and C is 

2sxy —2x*-+2x —sx —sy = 0. 
Interchanging x and _» gives the boundary between B and C. 

Hardy’s law is the statement that e = 0. Eliminating p, q and r 

from this and equations (2) produces 


(2-+5s)xy —2x —29+2—s = 0. 
This is the boundary between C and D. All three of these curves 
are hyperbolas. 
A point to note is that strong selection against (even lethality of) 
one homozygote does not necessarily cause the corresponding allele 


to disappcar from the population. Selection against the other homo- 
zygote can hold the heterozygote in the population. 


3. EXACT SOLUTIONS 


3.1. When x = _y = 1, equations (2) show that the equilibrium 
state satisfies sr—te = 0. This is a parabola in the (jf, q, r) trilinear 
diagram passing through the lower vertices of the reference triangle. 
Fig. 2 contains equilibrium parabolas for s = 0, 4, ? and 1. 

3.2. When x = y <1, then from (2) p = q, and therefore 


e = (2p—r)(2p+r) = 2p—r. 
Eliminating p, q and ¢ from (2) in this case leaves 
xr(1+s—sr) = (1 —r)(1—s-+sr), 


giving the frequency, 7, of the heterozygote in terms of the survival 
rate, x, of the homozygotes. Graphs of r against x for s = 0, }, 4, ? 
and 1 are shown in fig. 3. 


In the trilinear diagram the equilibrium point is 
(3(1 —r), 3(1 —r), r) 


on the perpendicular bisector of the base of the reference triangle. 
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3.3. General values of x and_y. Eliminating sr—de from equations 
(2) leaves 


aq(r+2p)=pirt+2q)- - ~« «+ (3) 
where ay(1—x) = x(1—y). 


Equations (3) are, for various values of a, a set of hyperbolas through 
the vertices of the reference triangle. 











Fic. 2.—Trilinear diagram showing approach to equilibrium in the absence of selection. 
Distances from the sides labelled p, g and r represent the proportions of aa, bb and 
ab respectively in the population. s is the proportion of selfing. 

The points of the parabolas represent possible equilibrium populations for the 
corresponding values of s. The approach to equilibrium is shown for three cases by 
arrowed lines, the dots indicating the position of the population in successive generations. 


Now there will be an equilibrium parabola through each 
equilibrium point. If it has equation 
(1—f)r—Be=o . ; ‘ « “ia} 


it must be satisfied by the solution of (2). The value of f for this 
to be so can be shown to be the positive root of 


B?(2xy —x —y) +B{(4 —2t)xy —(3 —t) (x+y) +2}+#(2 —x—y) = 0 
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Once a and f have been evaluated, we can solve (3) and (4) for 
the equilibrium point. The solutions are 
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with d-! = B(1+)(1+a)?. The appropriate solution is that corres- 
ponding to the area (A, B, C or D) in which (x, y) lies in fig. 1. When 
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Fic. 3.—Equilibrium with equal viabilities of the homozygotes. 

r = proportion of heterozygotes in the population. 

x = survival rate of both homozygotes compared with the heterozygote. 
s = proportion of selfing. 


s is given and (x, y) lies in C or D the equilibrium point lies above the 
parabola sr—te = 0, since the population must eventually become 
more heterozygous than it would in the absence of selection against 
the homozygotes. 
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Evidently when x and » are less than unity, equilibrium is 
independent of the initial state (unless one allele is initially absent) 
and is determined by the two parameters a and B. As these are 
functions of the three population parameters x, y and s, it is clear 
that different sets of values of x, y and s can determine the same pair 
of values of a and f, and hence the same equilibrium. This means, 
for example, that an increase in the proportion of selfing can be 
exactly counterbalanced by a decrease in the viability of the homo- 
zygotes. The diagrams of fig. 1 are sections for given s of a 3-dimen- 
sional graph of x, y and s, and in this graph each line of intersection 
of a pair of surfaces, given by equations (3) and (4) determines one 
equilibrium point. 


4. THE APPROACH TO EQUILIBRIUM 


4.1. When x = y = 1, the statements (1) lead to the well-known 
results for the case of no selection. The proportionalities now become 
equations and further sr, —te, = 45(sr9 —teg) 


“+ Pn = Pot im, > pot im = u?+suv/(2—S) 
In = Tothm, > Iotdm = v?+suv/(2—S) 


Tr, =To—-M, > %>—m == Quv(2—2s)/(2—s) 
where mM, = (ST9—teo) (1 —($5)”) /(2 —5) 
and m = (sf9—teg) /(2—S) 


m determines the limiting frequencies at equilibrium. The second 
form for the limiting frequencies indicates departures from a random 
mating population since u = py+4ry and v = gy+419, which are the 
gene frequencies. 

In 3.1 it was shown that the equilibrium point lies on the parabola 
sr—te =o in the trilinear diagram. Now p,+47r, = fpotdry and 
In+4%n = To-+4%, Which also follow directly from the constancy of 
gene frequencies (assuming no mutation). Hence the point (f,, ¢n) Tn) 
lies on the perpendicular line from (fo, qo, 7) to the base of the 
reference triangle. Evidently as n increases the population point 
proceeds along this line until it reaches the equilibrium parabola, 
the distance from equilibrium being reduced to 4s of its previous 
value every generation. Fig. 2 illustrates three cases when s = 0, 
4 and 1, the first attaining equilibrium in one generation. The dots 
indicate successive population points. 

4.2. For general values of x and y, successive applications of the 
proportionalities (1) enable the population point to be plotted from 
generation to generation, while the limit point may be determined 
independently by 3.3. Fig. 4 illustrates the approaches to equilibrium 
starting with two different populations F = (0-6, 0-3, o-1) or 
G = (0-1, 0°3, 0-6) under an equal mixture of selfing and random 
mating (s = 4). Commencing at F, the first few and some subsequent 
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generations have been indicated by numbered dots. Five sets of 
selection values (x, y) giving five equilibrium points have been chosen 
as in the table. 























Equilibrium points 
a b c d e 

ee —| | 
Selection values x . ; 0°259 0°750 0500 | 0600 0°950 
Sf - : 0°259 07500 | 0°750 | 0:g00 0-900 
Equilibrium frequencies p . O°125 0°634 | 0°057 0-000 0°724 
@. 0°125 0°057 0°634 | 1-000 0°072 
Fr 0750 | 0309 0-309 | 0-000 0°204 





These have also been labelled in fig. 1 for s = }. 
Some points may be noted. Most of the initial disequilibrium 








r 


Fic. 4.—Trilinear diagram of the approach to equilibrium when the proportion of selfing, 
$=}. 

F and G are two possible starting points while a, 5, c, d, and e are five equilibrium 
points corresponding to the values of (x, y) shown in fig. 1, s = 4. The solid lines 
show the paths of the population from F or G to a, b,c, d or e. The dots on the lines 
from F show the position of the population after the number of generations indicated. 

The broken line represents possible equilibria in the absence of selection. 
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between heterozygotes and homozygotes is rapidly removed: the 
disequilibrium between the two homozygotes disappears very slowly. 

With strong selection against the homozygotes the equilibrium 
is attained quickly and fairly directly. With weak selection against 
the homozygotes, the equilibrium parabola through the equilibrium 
point is first approached quickly and nearly vertically, thus correcting 
any heterozygote-homozygote disequilibrium. Then the population 
point moves slowly round close to this parabola to the final equilibrium 
point, thus correcting any disequilibrium between the two homozygotes. 

We recall that when there is no selection against the homozygotes 
the equilibrium depends on the initial gene frequencies. Hence, if 
fluctuations in these frequencies occur (by chance or by mutation) 
the equilibrium point shifts around the equilibrium parabola (deter- 
mined by s) since there is no controlling force to hold it at any 
particular point on the parabola. 

When there is only weak selection against the homozygotes a 
similar situation prevails. In case (e) and starting from (0-6, 0-3, 0-1) 
the population point is still distant 0-125 units from the equilibrium 
point after 30 generations and the next generation only reduces this 
gap by 2-1 per cent. This slow approach to equilibrium shows that 
the force working towards equilibrium is weak. This means that, 
though in general the equilibrium is stable and is fixed by the mating 
system and the survival rate of the homozygotes, there is a region 
near x = y = 1 in fig. 1 in which small disturbing influences may 
cause the population point to drift widely away from equilibrium 
along the equilibrium parabola. 


5. SUMMARY 


The equilibrium is investigated of a large population breeding 
under a system of mixed selfing and random mating with selection 
against the homozygotes. Diagrams are produced giving the types 
of equilibrium and illustrating the approach to equilibrium in 
particular cases. 
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1. INTRODUCTION 


NATURAL populations of most species of Drosophila are variable 
with respect to the gene arrangement in the chromosomes of different 
individuals. In some species, a considerable proportion of flies found 
in nature are inversion heterozygotes. The degree of the variability 
of the chromosome structure varies however in geographic races of 
the same species, as well as in different species. To account for the 
different variability of different races and species, da Cunha, Burla 
and Dobzhanksy (1950), and Dobzhansky, Burla and da Cunha 
(1950) have advanced the following working hypothesis. Since the 
chromosomal polymorphism in Drosophila populations is adaptive 
and balanced, the amount of this polymorphism present in a race 
or a species is likely to be a function of the number and variety of 
the ecological niches which the population has mastered. Other 
things being equal, species which exploit many different kinds of 
food, or which occur in many different biotopes, are likely to exhibit 
a greater degree of adaptive polymorphism than do species which are 
specialised to exploit only a limited variety of adaptive opportunities. 

The data reported in the present paper have been collected in 
order to submit the above working hypothesis to further tests. We 
have chosen several species of Drosophila which are favourable for 
studies on the chromosomes in the larve salivary glands. The amounts 
of the chromosomal variability were then compared either in pairs 
of closely related species, or in different populations of the same 
species. The data so obtained can be evaluated against the background 
of the ecological characteristics of the respective races and species 
described by Dobzhansky and Pavan (1950), Pavan (unpublished), 
and Cordeiro (unpublished). 


2. METHODS 


Single females collected in nature were allowed to produce off- 
spring in individual cultures with banana-agar nutrient medium. 
The salivary glands of mature larve produced in the cultures were 
stained with acetic orcein.. The chromosomes of a single larva in the 
offspring of each wild progenitor were then examined. 
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3. DROSOPHILA POLYMORPHA AND DROSOPHILA CARDINOIDES 


D. polymorpha Dobzhansky and Pavan and D. cardinoides Dobzhansky 
and Pavan are closely related species belonging to the cardini species 
group of the Subgenus Drosophila. They are sympatric in a large 
area, which extends, as far as known, from northern Argentina to 
Central America. The metaphase chromosome groups of both species 
are alike; they contain two pairs of metacentric autosomes (the 
second and the third), a pair of acrocentric X-chromosomes, and a 
pair of dot-like autosomes. The inversions found in these chromosomes 
are described by da Cunha and Salzano (in press). 


TABLE 1 


Frequencies (in per cent.) of various inversions found in heterozygous condition, and the mean 
numbers of heterozygous inversions per individual in D. polymorpha 



































Chromosome Inversions A A Pirassununga —o 
a —— -_ | — $$$ ee | — 
II { none 71°6 988 98°7 
( Aa 28°4 12 13 
none 92°5 100 100 
III Aa 75 
Bb 7°5 
none 43°3 63°6 82°5 
x Aa 43°3 30°4 17°5 
| Bb 56°7 36°4 17°5 
Ce 43°3 she 17°5 
{ im x 0°46+0°08 O-or o-or 
Mean per individual 
lo 1°62+0'18 r10+0-16 0°54+0'14 
: 3 67 83 93 
Flies tested 
67 88 63 











The chromosomal polymorphism in D. polymorpha is greater than 


that in D. cardinoides. 


In the former species, six different inversions 








have been found in 461 wild individuals collected in three localities, 
while in D. cardinoides three inversions were found in 141 flies, also 
from three localities. More important are the data on the frequencies 
of heterozygosis for inversions in the populations of the two species, 
summarised in tables 1 and 2. 

Sufficient samples of D. polymorpha were available from three 
localities in the State of Sao Paulo, Brazil, namely Vila Atlantica, 
Pirassununga, and Mogi das Cruzes. The mean numbers of hetero- 
zygous inversions, per female, were 1-82, 1-10 and 0-54 in the three 
localities. ‘The mean numbers in the males are lower, since no 
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inversion heterozygosis can be detected in the X-chromosome in the 


males. 


Samples of D. cardinoides were studied from Vila Atlantica (State 
of Sao Paulo), from Cantareira (same State, about 50 km. from 
Mogi das Cruzes), and from Catuni (State of Bahia). 
samples the inversions proved to be rare. 
zygous inversions per individual were 0-06, o and 0-05 in the three 
localities respectively. 

D. polymorpha is a species which is relatively more common, at 


least in southern Brazil, than D. cardinoides. 


In all these 


The numbers of hetero- 


At Pirassununga, the 


frequencies of D. polymorpha vary from 0-2 to 22 per cent. of the total 


catch of all species of Drosophila, depending upon season. 


D. 


cardinoides is much less frequent, varying from 0-03 to 0-5 per cent. 


TABLE 2 


Frequencies (in per cent.) of various inversions found in heterozygous condition, and the mean 


numbers of heterozygous inversions per individual in D. cardinoides 








| Cantareira | 


Vila 





























Chromosome Inversions Braganga Atlantica Catuni 
| | none 83 100 100 94 
II | <u ae 17 fs 
| | Bo a 6 
} 
f none 100 100 93 100 
an \ Aa “i 7 res 
Mean per individual 0-16 7) 0:06 -+0:04 0°05 +003 
Flies studied 6 50 46 55 








of the total catch. At Mogi das Cruzes, D. polymorpha varies from 
0-9 to 22 per cent., and D. cardinoides from 0-01 to 3°5 per cent. of all 
At Vila Atlantica, from o-2 to 3 per cent. of the flies 
are D. polymorpha and 0-1 to 4°3 per cent. are D. cardinoides. It should 
be noted that the high proportions of D. cardinoides have been observed 
close to human habitations, and the proportions in undisturbed 
localities are always below 1 per cent. 
at Catuni, D. cardinoides has been collected only in orchards, near 


Drosophila. 


human habitations. 


collected under such conditions. 


Similarly, at Cantareira and 


D. polymorpha is a species which we have not 


4, DROSOPHILA GUARAMUNU AND DROSOPHILA GRISEOLINEATA 


D. guaramunu Dobzhansky and Pavan and D. griseolineata Duda 
are closely related to each other, and both of them belong to the 
guarani species group of the subgenus Drosophila (King, 1947). Their 
known distribution regions extend from northern Argentine to central 
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Brazil. In both species, the metaphase chromosome plaics contain 
five pairs of rod-shaped and one pair of dot-like chromosomes. A 
total of 16 kinds of inversions have been found in the offspring of 
896 wild individuals of D. guaramunu from seven different localities. 
A description of these inversions can be found in the paper of Brncic 
(in press). 


encountered in the offspring of 612 wild flies from six localities. 


inversions are to be described by Salzano (unpublished). 


TABLE 3 


In D. griseolineata, 5 kinds of inversions have been 


These 


Frequencies (in per cent.) of various inversions in heterozygous condition, and the 


mean numbers of heterozygous inversion per individual of D. guaramunu 























Chromo-| Inver- | Ponta . | Mogi das Piras- 
scene ee Chapada | Grossa Itapoan Feliz Braganga Tees sununga 
I { None 60:8 ras 56°4 60-0 48°5 45'0 50°8 47°5 

Aa 39°2 43°6 40°0 515 55°0 49°2 52°5 
( None 43°4 41°8 50°0 515 45°0 49°2 47°9 
III (4 Aaor Bb 506 | 471 50°0 41°0 50°0 43°4 47°9 
| Aa Bb bse | 10°9 os 75 5'0 74 41 
None 32°6 | 74°5 30°0 30°3 30°0 49°2 40r! 
Aa 19°5 3°6 10°0 212 25'0 25°4 20°6 
Bb I5r | 10°9 ana 18+1 25°0 19°7 15°2 
Ce ake - oat ce hin o'2 
Dd mE Pe ses 4°5 50 08 2°4 
Ee wer 606 | 14'S 30°0 36-2 30°0 08 13°9 
Iv |) & se | we Hr cay ee O'4 
Gg 4°3 3°6 ae: I° | 10:0 6°5 3°7 
Hh 14'I 14°5 30°0 31°8 =| 30:0 0:8 12°6 
li 5°4 os se 20°0 1°97 56 
i 5°4 20°0 1°7 56 
Kk 5°4 | 200 1*7 5°6 
Ll 5'4 20°0 1°97 56 
Mm ase | 1*7 “ 
ab. S| es eee el eee el ES - | — = 
Mean per 213013 | 1°65+0°14 | 1°60+0°38 | 2-18+.0'17 | 3:20-+-0°46 | 1°66+-0:09 00+0:06 
individual 
—_—— ——| ——$—$—$— | —__ ———_—_- - _ —| — 
Flies tested 92 | 55 10 66 20 192 | 461 





























The frequencies of the various inversions in the population samples 














studied are summarised in tables 3 and 4. It can be seen at a glance 
that the chromosomal variability is much greater in D. guaramunu 
than in D. griseolineata. In the former species, the localities in the 
State of Rio Grande do Sul (Chapada, Ponta Grossa, Itapoan and 
Feliz) have between 1-6 and 2-2 inversions per individual. In the 
State of Sao Paulo (Braganga, Mogi das Cruzes, Pirassununga) the 
inversion frequencies are about the same as in Rio Grande do Sul 
(1: 7 to 3:2, the latter figure being subject to a large error). Since 
no inversions have been detected in the X-chromosome of these 











species, the data for females and males can be combined. 
griseolineata, the State of Rio Grande do Sul (Ponta Grossa, Itapoan, 
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In D. 


Emboaba) shows perhaps fewer inversions per individual than does 
the State of Sao Paulo (Vila Atlantica, Cantareira, and Mogi das 
Cruzes). But even the highest recorded inversion frequency per 
individual of D. griseolineata (0-47, cf. table 4) is lower than the lowest 


frequency found in D. guaramunu (1-6, cf. table 3). 


According to Dobzhansky and Pavan (1950), D. guaramunu is 
more common and widespread in Brazil than is D. griseolineata. The 
former species proved to be one of the commonest species of Drosophila 
in 5 out of the 35 localities studied by these authors. The latter species 


TABLE 4 


Frequencies (in per cent.) of various inversions found in heterozygous condition, and the mean 


numbers of heterozygous inversion per individual of D. griseolineata 


























| | 
Chromo- Inver- | Ponta | Vila : Mogi das | 
some sions ar Grossa Itapoan | Emboaba Atlantica Cantareira | Cruzes 
none 96°5 97'I 986 80°8 go°6 | 829 
II * | - ag se I'l was 2°6 
3°5 29 1°4 14:0 | 8-0 | g'2 
acy va 4°1 | 14 | 53 
| 
— 94°7 1000 100°0 86-4 | 81:0 =6|)~— 868 
II! 53 a5 ies 13°6 19'0 140 
53 13°2 190 14°0 
ae | eee = es Sr 
Mean per individual 0'14+0:06 | 0-02+0°06 | o-or-+0-04 | 0-47+0°04 | 0°47+0°09 | 0-46+0-09 
se oat Saas TR a Metaiahe 
Flies tested 57 69 7 265 | 74 76 
| | | 














has been found in 6 out of the 35 localities, but in all of them it was 
one of the rarer forms. D. guaramunu is, indeed, one of the dominant 
species both in the savanna of central Brazil and in the subtropical 
forest of the State of Rio Grande do Sul. At Pirassununga, which 
represents ecologically a southern outlier of the savanna of central 
Brazil, the lowest frequency of D. guaramunu in 18 different collections 
made during four years of observation was about 2 per cent., and in 
9 of these collections it was above 18 per cent., of the total catch of 
Drosophila flies. At Mogi das Cruzes, its frequency varied from 
0 to 23 per cent. of the total, and at Vila Atlantica from o to 6 per 
cent. of the total. In Rio Grande do Sul, D. guaramunu is encountered 
with frequencies from 0-5 to 5 per cent. on the coast, and from 11 to 
51 per cent. in the interior of the State. 

D. griseolineata varies greatly in frequency, but it is on the whole 
much less common than D. guaramunu. The recorded frequencies of 
the former species vary from o to go per cent. at Vila Atlantica and 
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from o to 16 per cent. at Mogi das Cruzes. In Rio Grande do Sul, 
this species is relatively more common on the coast (2 to 30 per cent.) 
than in the interior, where it is quite rare. Furthermore, the frequencies 
of this species are subject to violent fluctuations at different seasons. 


5. DROSOPHILA NEBULOSA 


Drosophila nebulosa Sturtevant is a member of the willistoni species 
group of the subgenus Sophophora. Its distribution region extends 
from Texas in the North (Patterson and Wagner, 1943) to at least 
Buenos Aires in the South (Pavan, unpublished). No Drosophila 
collecting has ever been made south of Buenos Aires. In Brazil, 
D. nebulosa is one of the commonest species, particularly in drier 
savanna and desert regions, less common in the rain forests. Further- 
more, the frequencies of D. nebulosa relative to those of other species 
of Drosophila tend to be higher during the warmer and drier than 
during the wetter and cooler seasons (Pavan, unpublished). 

The metaphase plates of D. nebulosa show two pairs of metacentric 
and one pair of acrocentric chromosomes. Pavan (1946) examined 
the salivary gland chromosomes of a series of strains, and described 
9 inversions, all but one of which were found in a single chromosome, 
the third, which is the acrocentric autosome. In our material, we 
have observed 7 out of the 8 inversions described by Pavan in the 
third chromosome, and found 3 new inversions, namely fi, Jj, and 
Kk, also in the third chromosome. The single inversion found by 
Pavan in the X-chromosome was not found in our material. 

The frequencies of the inversion heterozygosis in the populations 
examined by ourselves are given in table 5. Only the gene arrange- 
ments HH, Hh and hh were identified both when homozygous and 
when heterozygous. The inversions Aa, Bb, Cc, Gg and Hh occur 
frequently in all the regions studied ; the inversions Ee, Ff, li, Zj 
and Kk are less common or restricted to only small fractions of the 
species distribution area. 

It may be noted that a group of populations comprising Mucajai, 
Monjolinho, Pirassununga, Iguassi and Tucuman show high frequencies 
of the gene arrangement f, while Rio Negro, Imperatriz, Salitre 
and Catun{ are characterised by high frequencies of the alternative 
gene arrangement, H. This distribution of the frequencies of h and 
H shows no clear correlations with the ecological characteristics of 
the localities in which the respective samples were collected. Thus, 
Mucajai is more similar to Monjolinho, some 2500 km. away, than 
it is to Rio Negro, less than 500 km. away. Mucajai, Iguassi and 
Rio Negro are covered with rain forest, while the other localities 
have savanna or desert vegetation. 

The mean frequencies of the numbers of heterozygous inversions 
per individual show, on the contrary, a fairly clear correlation with 
the ecological characteristics of the regions in which the population 
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samples were taken. The highest frequencies are observed at Mucajai, 
Rio Negro, Monjolinho, and Imperatriz. The first and the second 
of these localities are covered with exuberant rain forests, but are 
close to margins of savannas. On Rio Negro in particular, D. nebulosa 
proved to be commoner in drier localities than in wetter ones (up 
to 50 per cent. of the total catch of Drosophila flies on the beach of 
Tatuaquara, near Manaus). Monjolinho and Imperatriz lie in 
ecologically highly complex regions, on or near the ecotones between 
the savanna and the rain forest, or gallery forest, vegetations. By 
contrast, Salitre and Catuni, which show low frequencies of inversions, 
are in a desert (caatinga) region, which is dry and nearly lifeless 
during more than six months per year. Although D. nebulosa is the 
commonest species of Drosophila in this desert, it certainly can occupy 
there only a very limited number of ecological niches. It should be 
noted that D. willistont has also shown the lowest numbers of hetero- 
zygous inversions per individual in this region (da Cunha, Burla and 
Dobzhansky, 1950). Finally, Tucuman lies at the foot of the Andes, 
presumably near the margin of the distribution area of the species. 
The population sample was collected there in an isolated forested 
region having a limited variety of plant species. 


6. DROSOPHILA BANDEIRANTORUM AND DROSOPHILA PALLIDIPENNIS 


Drosophila bandeirantorum Dobzhansky and Pavan and D. pallidipennis 
Duda are species with no clear affinities among the known species of 
Drosophila. The former has a rather restricted known geographic 
distribution area in southern Brazil. The latter has been recorded 
from the State of Sao Paulo, in southern Brazil, as well as from Mexico, 


but has so far not been encountered in the Amazon Basin. Both . 


species are fairly rare, and appear to be specialised to occupy a 
restricted variety of ecological niches (Dobzhansky and Pavan, 1950). 
Their salivary gland chromosomes are among the best found in any 
species of Drosophila, and this served as a stimulus to study them 


for inversions. 
In 441 examined individuals of D. bandeirantorum, two kinds of 


inversions were found, one in each, the second and the third chromo- 


somes (autosomes). In D. pallidipennis, Patterson and Dobzhansky 
(1945), Freire-Maia and Engel (1949), and ourselves found two 
inversions, both of them in the same autosome (overlapping inversions). 
The mean frequency of heterozygous inversions per individual of 
D. bandeirantorum proved to be 0:49+0-02. The material examined 
came from 4 populations in the State of Sao Paulo. In the 70 
individuals of D. pallidipennis, from Pirassununga, in the State of 
Sao Paulo, examined by Freire-Maia and Engel (1949) and ourselves, 
the mean number of heterozygous inversions per individual proved 
to be even smaller, namely 0-:07-+-0-04. 
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7. DISCUSSION 

The data presented above are in a general agreement with the 
hypothesis that the amount of genetic variability present in a popula- 
tion is positively correlated with the number of ecological niches 
which this population exploits. Thus, the common and widespread 
species, D. polymorpha, has more heterozygous inversions per individual 
than its closest but less common relative, D. cardinoides. Similarly, 
the common D. guaramunu has more inversions than its less frequent 
relative, D. griseolineata. In D. nebulosa the populations which inhabit 
the ecologically limited and highly exacting desert environments 
show fewer inversions than do the inhabitants of richer and more 
diversified rain forest and savanna environments (cf. da Cunha, 
Burla and Dobzhansky, 1950). 

Several reservations must, however, be made. In the first place, 
only closely related species, or races within a species, can profitably 
be compared. In different groups of organisms the genetic variability 
takes different forms, and chromosomal inversions are characteristic 
of the adaptive polymorphism in some sub-divisions of the animal 
and plant kingdoms but not in others. Even within the genus 
Drosophila, some species groups seem to be richer in inversions than 
others (thus, the subgenus Sophophora seems to have more inversions 
than the subgenus Drosophila). Some species have inversions in all 
chromosomes (e.g. D. willistont and relatives, D. robusta, D. virilis), 
while other species have them concentrated in a single or in few 
chromosomes (e.g. D. pseudoobscura, D. persimilis, D. nebulosa). 

Furthermore, given the present imperfect understanding of the 
ecology of Drosophila, our ability to evaluate the relative numbers 
of the adaptive niches controlled by a population is limited. By and 
large, an abundant and widespread species is likely to be ecologically 
more versatile than a rare one restricted to only few habitats. It 
stands, however, to reason that an ecologically specialised species 
may occur more commonly than an ecologically versatile one, provided 
that the special environment of the former is available in many places 
within a certain geographic region. Situations apparently contra- 
dictory to the hypothesis are, then, not unexpected, but, provided 
that the hypothesis is correct, they should be less frequent than 
situations conforming to the hypothesis. Such an apparent contra- 
diction may be seen in table 1. D. polymorpha has significantly more 
inversions per individual at Vila Atlantica than at Pirassununga or 
at Mogi das Cruzes. And yet this species reaches at times much 
higher frequencies at the latter than at the former localities. It may 
be significant in this connection, that the population size undergoes 
violent fluctuations at Pirassununga and at Mogi das Cruzes, reaching 
up to 22 per cent. of the total Drosophila catch during the summer, 
but becoming very rare during winter. At Vila Atlantica the 
population of D. polymorpha does not undergo such sharp expansions 
and contractions at different seasons. 

Oo 
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8. SUMMARY 


Natural populations of several species of Drosophila were examined 
for the incidence of inversion heterozygotes among the individuals 
composing them. Inversion heterozygotes proved to be relatively 
more frequent in Drosophila polymorpha than in its close relative, 
Drosophila cardinoides, and more frequent in Drosophila guaramunu than 
in the neighbouring species, Drosophila griseolineata. Drosophila poly- 
morpha and Drosophila guaramunu are more common and widespread 
in southern Brazil than are Drosophila cardinoides and Drosophila 
griseolineata. In Drosophila nebulosa, populations which inhabit eco- 
logically rich and diversified environments of savannas and rain 
forests have more inversions than do geographically or ecologically 
marginal populations. These facts stand in accord with the hypothesis 
that the amount of adaptive polymorphism present in a population is 
positively correlated with the variety of environments which this 
population controls. 
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1, INTRODUCTION 


THE use of monosomics, 6x —1 plants, for genetical studies and for the 
transference of desirable characteristics from related genera into 
the hexaploid wheats (Triticum vulgare) has become an important 
part of the cereal breeding programme. Some particular cytological 
studies of the monosomes (the deficient chromosomes) have been 
made: Sears (1939, 1944, 1952, 1952a), Unrau (1950), Unrau et al. 
(1950) and Morrison and Unrau (1952). As suggested in an earlier 
report there are several major problems requiring clarification. 

The B-series of speltoid wheats as classified by Nilsson-Ehle in 
1921 has been associated with the loss of one chromosome (Winge, 
1924). A review of work on this strain of wheat has been given by 
Huskins and Sander (1949). In the Sear’s monosomic series the 
speltoid chromosome is number nine (mono-IX). Much of the 
behaviour of the other monosomes has been inferred from the reported 
behaviour of this monosome. The project reported in this paper was 
undertaken with the view of comparing the meiotic behaviour of all 
the monosomes. 

One phase that is not yet clear is the behaviour of the chromosomes 
in the pollen grains (P.G’s.). From time to time inferences have 
been made that P.G’s. with varying numbers of chromosomes are 
produced or that telocentrics or isochromosomes may be transmitted 
through the male gamete, but no serious attempt has been made to 
study this phase of development. By the use of a new technique it 
was possible to make an extensive record of the behaviour of the 
chromosomes in the P.G’s. 

The morphology of the various monosomes was also compared. 
This was done to ascertain if any homologies were present and also 
as an aid to identification of telocentric or fragmented chromosomes 
in the P.G’s. or in the somatic complement of root-tips. 


2. METHODS AND MATERIAL 


All investigations of P.M.C’s. were made from acetocarmine squashes after 
whole spikes were fixed in Carnoy. Root-tips were pretreated with monobromo- 
naphthalene, fixed in 2 BD, and then prepared as Feulgen squashes. Pollen grain 
preparations were made using the technique outlined previously (Morrison, 1953). 
For studies of megasporogenesis, ovules that had been fixed in Carnoy were embedded 
in paraffin and stained with crystal violet. Sections were cut at 12 microns. 
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The 21 monosomic lines of wheat are those developed and distributed by Sears, 
i.¢. monosomics of T. vulgare var. Chinese Spring. The chromosomes are numbered 
from I to XXI and the deficient monosomes abbreviated to mono-I, mono-II, 
etc. Mono-XIV was not available for all studies. 


3. CHROMOSOME MORPHOLOGY IN T. VULGARE 


With the exception of the satellite chromosomes I and X it was 
not possible to identify each of the 21 different chromosomes by 
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Fic. 1.—The morphology of the chromosomes of Chinese Spring; XIV is omitted. 
Drawings were made from univalent laggards at TII. Magnification x 1300. 


























morphology studies of the root-tip metaphases of the monosomics. 
A comparison of the chromosomes was made however from the lagging 
univalents at telophase II (TII). The sizes and shapes (fig. 1) agree 
in general with the haploid set observed in the P.G’s. and with the 
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root-tip complement. It has long been known that IX has unequal 
arms. According to Sears and Camara (1952) VII is also unequal. 
Mono-I is sometimes visible at MI as a satellited univalent. 

The methods used in the past for studying the morphology of the 
chromosomes have been laborious and not always fruitful. Differences 
in sizes are not outstanding and divergence in arm length occurs to 
the same degree in too many chromosomes. As well, all reports 
are conditioned by cytological procedures. In the earlier works 
(Kagawa, 1929; Bhatia, 1938; and Levitsky e¢ al., 1939) two or 
three constrictions per chromosome were noted. Three pairs of 
satellite chromosomes have been recognised by Bhatia (1938), Pathak 
(1939) and Camara (1944) and two pairs by Levitsky e¢ al. (1939) 
and Tjio and Levan (1950). 

At diplotene-diakinesis I have observed from one to three bivalents 
attached to the nucleolus. This indicates that some organising 
regions are weaker than others and only operate occasionally. In 
like manner only rarely are six nucleoli seen in a resting nucleus. 
The constriction in chromosome X is never visible at MI and is not 
constant at AII. This variability in appearance and function of the 
nucleolar arm may thus be an expression of the evolutionary change 
in this polyploid (cf. Darlington, 1937, p. 228 and Woods, 1937). 
If Larson (1952) is correct and I and X are both from one genome 
then the satellite chromosomes from the other genomes must be 
considerably changed in appearance. 

From the compensating effects of tetrasomics, Sears (19525) has 
suggested several homceologous series. The limitation of attempting 
a comparison by shape of chromosomes is immediately apparent in 
the series I, XIV, XVII. Homeeologies are uncertain because of 
evolutionary changes. 


4. MEIOSIS IN POLLEN MOTHER CELLS 


Metaphase pairing —At MI the 41 chromosomes usually form 20 
bivalents and one univalent (table 1). Contrary to Sanchez-Monge 
and Mac Key (1948) no higher associations were found. There is 
no evidence for any serious degree of variation among the monosomes. 

The asynaptic condition observed in four cells of mono-IX is the 
result of either an environmental stimulus or a mutation in the sporo- 
genous tissue. It is not comparable to the desynapsis in nulli-III 
plants or to the asynapsis reported by Huskins and Hearne (1933) 
or Li et al. (1945) where all the cells within an anther showed a varying 
amount of failure in pairing. 

Univalents at AI and TI.—At AI, the univalents may lag or they 
may not; they may divide or misdivide, and at TI they may be 
included or excluded. Frequently after division both chromatids 
are swept to one pole. Division of the monosome occurs in only 
slightly more than half of the cells. Mi§sdivision of the centromeres 
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resulting in telocentrics or isochromosomes takes place in 12 per cent. 


of the cells (table 2). 


With the exception of mono-XII which is 


slightly low, there are no great differences among the monosomes. 























TABLE 1 
Chromosome pairing at metaphase I in eight monosomic lines 
MI pairing 
Monosomic No. of cells 
examined 
201 31 | 1g! 3} 18! 51 
| 
I . 194 4 wee 198 
II 72 5 aba 77 
IV 385 6 ae 391 
VIII 88 5 3 96 
im 161 4 di 165 
XV 129 2 I3I 
XVI 210 4 I 215 
XIX 168 I 169 
Total 1407 3r 4 1442 
Per cent. 97°6 | 21 03 














* In another sample from this plant in which most of the cells were at interphase, 


there were three delayed cells at MI where pairing had completely failed. 
cell pairing had been reduced to two bivalerrts and 37 univalents. 


TABLE 2 
Frequencies of misdivision and formation of micronuclei at 


interphase in eleven monosomic lines 


In one other 



































Misdivision of univalents 
= - wid Micronuclei 
| in dyads 
TI Tit 
Monosomic Ran =a . is : eee | 
| | 
No. of cells | | Per cent. No. of cells oe No. of cells P yo 
examined | misdivision examined misdivision examined micronuclei 

I 64 9 4 60 33°3 357 389 

III cel 52 27°0 286 4o°6 

V 87 149 23 17°4 327 43° 
VIII 152 13°8 46 13°0 325 31°7 

IX 13°2 73 219 293 341 

pid 120 10°'0 28 25°0 534 33'°9 

XI 255 13°3 135 mr 640 43°9 

XII 103 58 vee tes 475 42°3 
XIII 112 116 106 14°2 438 49°1 

Ye . eee | eee eee eee 225 53°7 
XI 53 | 17°0 70 II'4 ape 

Total 1014 | 12‘ | 593 17°77 3900 | 4r°0 | 





The frequencies of misdivision previously reported for mono-IX 
are 1*7 per cent. (Sanchez-Monge and Mac Key, 1948) and 39-7 per 
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cent. (Sears, 1952). It is evident that different methods of scoring 
have been used. These workers as well reported that the univalent 
divided in nearly all the cells (96 and 97 per cent.) which is not in 
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Fic. 2.—The chief variations in behaviour of the univalent chromosomes at meiosis and 
the resulting chromosome numbers in the P.G’s. Each phase is described in the text. 
Micronucleus (mn) formation in the tetrads has been described previously (Morrison 
and Unrau, 1952). 
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agreement with my results. In any event misdivision frequencies 
cannot be estimated too accurately because of questionable phases 
and observations at critical times of division. 
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Interphase.—Univalent laggards or fragments from misdivision not 
included in the TI nuclei form micronuclei (table 2). One or two 
and less frequently three or four micronuclei per dyad were observed. 
There is no correlation between the frequencies of misdivision and 
the proportion of dyads with micronuclei among the monosomes. 

There is, however, a positive correlation between the number 
of dyads with micronuclei and the number of tetrads with micronuclei 
(of. Morrison and Unrau, 1952). This confirms my original assumption 
that the inclusion or exclusion of chromosomes from the telophase 
nuclei is dependent upon the specific behaviour of the univalent. 

Second division—Chromosomes at AII were counted in some cells 
(table 3). Actual counts were frequently limited to one sister half 


TABLE 3 


Distribution of the chromosomes at anaphase II 


























Anaphase II nuclei with various | 
chromosome numbers —- 
Monosomic er 
groups 
19 | 20 20+f | 20+2ff 21 22 counted 
Mono-III . ee 49 I I M ass 58 
All others. 2 49 3 ais 29 I 84 
Total . 2 98 4 I 36 I 142 
Percent. . | I°4 | 69:0 28 07 25°4 0°7 | 
| | 




















because of non-synchronisation in both cells. In some cells the 
chromatids or fragments from misdivision of univalents were in- 
corporated into the polar groups. More frequently the lagging 
univalents were lying between the two poles in which case they 
were not counted. It is impossible to tell if they will be included 
or excluded from the microspore nucleus. Although AII is not a 
satisfactory stage at which to make a great many counts or a com- 
parison between monosomes it does given an accurate picture of the 
distribution of the chromosomes into microspore nuclei. 

Univalent chromatids that arise through equational division at 
AI cannot divide again in the second division. Like the products 
of the delayed first division they may be included, excluded or they 
may misdivide. The increased misdivision at TII (table 2) will 
result in a general increase in telocentrics. Isochromosomes formed 
at first division may be separated into telocentrics. However, reunion 
of broken ends in the resting microspores may restore both iso- 
chromosomes and monosomes if both telocentrics are included within 
one nucleus. ; 

At MII some of the micronuclei formed at the first division become 
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active. They behave as lagging chromatids. Other micronuclei 
which are isolated near the periphery of the cell are delayed in 
development. They have the appearance of prophase chromosomes 
when the others are undergoing division. They are not included in 
the telophase nuclei but again form micronuclei. Most of the micro- 
nuclei have disappeared from the P.G’s. by the time the cytoplasm 
becomes vacuolated. Presumably the micronuclei are absorbed by 
the cytoplasm. 

Whole chromosomes or fragments that are delayed in division 
are sometimes isolated from the tetrads as microcytes (cf. Frankel, 
1949). These microcytes develop into small pollen grains. Although 
only one-eighth as large, they are similar in appearance to the normal 
P.G’s. but have no divisions. 


5. POLLEN GRAINS 
Metaphase.—Although varying chromosome numbers in the P.G’s. 
have been postulated to account for progeny results this is the first 
report of actual observations (table 4). The number in any P.G. 


TABLE 4 
Distribution of the chromosome numbers at 1st division in the 
P.G’s. for twenty monosomic lines 















































P.G’s. with various chromosome numbers 
| | No.of | 
Monosomic 19 | 20 2: P.G’s. 
Beet & a |__| | examined 
18 — sams 22 | 23 
| | | | 
Jor| 24] 0 of | aplse| of |xs| | 
mn TRE) BR erie ORE VS, FS Sa 
ie : sabeae me gy v Peyprye eae cae 25 
im. ; Pts | casi Wexad gr | 1 Bip icee* GEE Pcait 50 
III » | one | 40 | ral eee be Loe acaba 51 
Iv* »| 3 | I 92 By ous | 29 | a, 128 
Vv. tis a 8 gy Bree ee a 2 Slbtedel 28 
VI fa Nica 16 I L #3). “eh Peet 30 
VII Be id ae | ag} ult | mail 29 
VIII * eS 102 | 8/ 1 Poel Sab Rak oe 
IX * Bis I 78 Re hI 5! ie ere ee 14! 
x* ae 39 | 4| 3 35 1] 83 
XI ea | as acre I a - |e 41 
XII ois 34 iad 5 ; | I 40 
aC m Pisael iene Bces 35 I I 34 Oe a 72 
ave. « Uscas sb oe oaall -oAee |p NOD lever tees 16 Oy nereal 75 
bo!) ae A al We | ae rs ene aes 3 Rie | 25 
XVII*. AR REY Ge 72 [al ey re | go | P | : 107 
XVIII . a) he Gree 18 2 10 | ia oe 3I 
> a yh Rech UB haa 52 6 15 1 | 76 
XX 4 o | cee | coe | one OF TIT OH uc 15 | 93 
> 4 ~ Phace PORE vas 24 I I | | | 28 
| 
| 
Total . 1 4 | 8 | 2 8x12 | 65 | rz | t | 393 | 4 | 6 | 3 1309 
Per cent. of Total | 0-3 | 0-6 | o-r | 62:0 50 08|o-1r| 30:0 | 0-3} 0-5 | o-2 
| | | | } | 








Mean number of chromosomes per P.G. = 20°3. 
* Combined results for two years. 
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depends upon the inclusion or exclusion of lagging univalents in that 
particular microspore nucleus. 

The 1951 results showed significant differences among the mono- 
somes for the formation of P.G’s. with 20 or 21 chromosomes. From 
more extensive studies in 1952 it is evident that results for the two 
years vary considerably (table 5) and that no specific ratio can be 
established from the results of the small samples. 


TABLE 5 


Distribution of the chromosome numbers at 1st division in the P.G’s. 
Sor two years, for mono-IX and mono-X 





P.G’s. with various chromosome numbers 


Mono. | ‘sie a: Mee No. of | No. of | Mean 

year 19 e P.G’s. | chrs.* | per P.G. 

So al ane eau I... | 23 
of|1flof\1f\efiof\:f| | 
| 7 fon |_ = “| 4 a 








es eee Pee ae ae wee kag Pee ee 41 842 20°5 
IX-1952 | 1 | 1 ee oe. eS Ry ae See eee 100 2026 20°3 
| 
X-1951 | ... | ae | oe POT 2 1 | 7 | 4 | I 33 672 20°4 
| | 
Matone| ous bose |csent Oe 1 50 1028 20°6 














* Omitting fragments. 


Some of the apparent heterogeneity among the different mono- 
somes may be due to differences in the stage of anther development. 
If the P.G’s. with 20 chromosomes are sufficiently unbalanced they 
will be delayed in coming to mitosis. Then if the anthers are examined 
late in the cycle a bias in the favour of the 20-chromosome P.G’s. 
will be introduced. No estimate of this bias could be made by 
counting the number of uni-nucleate and bi-nucleate P.G’s. present. 

Disregarding fragments and using the total for all monosomes 
(table 6) the ratio of the various numbers agrees with results for AIT. 
As expected, the frequency of 18, 19, 22 and 23 is low, but these 
P.G’s. are produced in sufficient numbers to warrant attention as 
possible gametic sources of variation. 

Because of the loss of chromosomes as micronuclei, more P.G’s. 
with 20 than with 21 chromosomes are expected. Even in the tetrads 
with no micronuclei two of the microspores will be deficient. As 
50 per cent. of the tetrads have one or two micronuclei (Morrison 
and Unrau, 1952) the ratio of P.G’s. with 21 to 20 chromosomes 
should be approximately 25 : 75. My results are not greatly divergent 
from this ratio. 

Deficient chromosomes were tabulated as fragments because it 
was not always possible to discriminate between telocentrics and 
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chromosomes with ruptures near the centromere. Of all fragments, 
60 per cent. were certain telocentrics, another 15 per cent. were 
inferred and 25 per cent. were doubtful. There is no correlation 
between misdivision frequencies and number of fragments observed 
for the various monosomes. 


TABLE 6 


Distribution of the chromosome numbers and a comparison of frequency of 
fragments at 1st division in the P.G’s. 
























































P.G’s. with various chromosome numbers | 
No. of 
g | No. of Gf 
18 19 20 21 22 23 P.C’s. 

; — fectine 
o 4 8 812 393 6 3 1226 o.4 
I wat 2 65 4 a aes 71 i 
2 “ts II sia de2 Ss II ae 
3 I eee eee eee I 3 

Total . , 4 10 889 397 6 3 1309 96 | 
No. of cells with f'| ... 2 77 4 re a 83 | 
ff per cent. ‘ Ae 20°0 8-7 10 ses mae 6-3 | 

TABLE 7 


Anaphase and telophase of 1st division in P.G’s. in 
six monosomic lines 


























Anaphase or telophase in the P.G’s. 
Monosomic | 
Normal with : Tota 
no laggards With laggards examined 
me , ° 68 | 2 70 
VII : ; 7 I 8 
VIII ‘ ; 65 I 66 
XV ‘ ‘ 18 fC) 18 
XX . é 34 | 2 36 
XXI ; : 65 ) 65 
Total . : 257 | 6 263 
| | 
97°7 percent. | 2-3 per cent. 








No actual isochromosomes were recorded although two probables 
were seen. It would be difficult to distinguish them from chromo- 
somes with median centromeres except those duplicated for the 
trabant arm of I or X. In three P.G’s. chromosomes deficient for 
only a portion of an arm were seen. 
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An examination of the second mitotic division showed that P.G’s. 
with 20 and 21 chromosomes were undergoing division. Results 
were not tabulated because the numbers of observations were limited. 
The smallness of the chromosomes also prevented an accurate analysis. 





Fic. 3.—Four pollen grain mitoses, odd chromosomes stippled. The P.G. wall is not 
drawn to scale. Magnification x 1300. 


Anaphase and Telophase.—Studies were made of anaphase and 
telophase in six monosomic lines (table 7). No attempt has been 
made to compare the various monosomes statistically because of the 
small numbers for some monosomes. The frequency of abnormal 
behaviour calculated from the total is probably near correct for all 
monosomes. 

The six laggards can be classified into two types :— 

(1) Cells with extra nuclear material out of phase. Micronuclei 
carried over into the P.G’s. become active but are delayed 














WHEAT MONOSOMICS 213 


in development (fig. 3p). The action of the micronucleus 
is comparable to the action it frequently undergoes at 
interphase. The two laggards in mono-XX were of this 
kind. 

(2) Cells with centric fragments or chromosomes that have 
sub-effective centromeres (some may be acentric fragments) 
which were included in the microspore nucleus at TII. 
The other four were of this kind (fig. 4). 


Darlington (1940) proposed non-disjunction of telocentrics at 
P.G. mitosis to account for the production of isochromosomes in 
Fritillaria. Sanchez-Monge and Mac Key (1948) have proposed a 
similar mechanism in wheat. No support for this theory was found. 
Nor is there any support for the second postulate of Rhoades (1940) 
whereby isochromosomes produced at meiosis become micronuclei 
which are carried over to the P.G. and then become active. As 
Sears (1952a) has suggested, isochromosomes and telocentrics arise 
directly through misdivision at meiosis. 


6. MEIOSIS IN EMBRYO SAC MOTHER CELLS 


Observations of meiosis were made from sectioned ovules of 
mono-VII, XVI, XX and XXI. The behaviour of the univalent 
(fig. 4D) is similar to its behaviour in the P.M.C’s., 2.e. it lags, divides 
and forms micronuclei, etc. The plane of division in one dyad is at 
right angles to the other however, and this results in a T-shaped tetrad 
of megaspores (fig. 45). 

The wheat gametophyte is monosporic and it is always the inner- 
most celi (cell number 4) that develops into the embryo sac. Micro- 
nuclei are found in all four megaspores. Out of a total of 48 tetrads 
79°2 per cent. had one or more micronuclei. In mono-XX out of 
27 tetrads 81-5 per cent. had micronuclei. 

If no micronuclei are present we cannot be certain whether the 
terminal megaspore or the first two have 21 chromosomes. Again, 
if there is one micronucleus in the third megaspore it does not indicate 
the deficiency or completeness of cell number four. Probably the 
best check on the formation of deficient or normal gametes can be 
made from crossing and selfing experiments. Sears (1944) has shown 
that in general 75 per cent. of the female gametes are deficient. 


7. PROGENY 


A record was kept of the chromosome number of all plants grown 
to maturity (table 8). Results for individual monosomics are not of 
sufficient magnitude to warrant a comparison. Of the total, the 
ratio of normal to monosomic plants is slightly higher than expected 
but the number of nullisomics compares favourably with previous 
reports (Sears, 1944 ; and Smith et al., 1949). 


——<— See 
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As the amount of zygotic elimination and embryo abortion is 
practically nil we must assume excessive certation, partial pollen 
sterility at later stages or slower pollen tube growth of defective P.G’s., 
to account for the ratio. With a ratio in the female of monosomic 
to normal as 75 : 25, then P.G’s. with 20 chromosomes can be only 
one-thirtieth as efficient as the P.G’s. with 21 chromosomes or there 
would be many more nullisomics produced. 


TABLE 8 


Progeny grown from monosomic seed 














Type of plant and chromosome number | 

Blasi Monosomic | Nullisomic Other | Total 
eee, orn puiniiicinial didetaaatiais | Aberrants * | 
| - | 

| No. of plants . 68 101 2 6 | 177 
wt | 
Per cent. ; 984 571 rr 34 | 











* These plants were as follows :— 


i, 1 with 40 plus an isochromosome. 
ii. 2 with 41 plus an isochromosome. 
iii, 2 with 41 plus a telocentric. 
iv. 1 with 40 plus a telocentric plus a chromosome with an interchange that arose 
through fusion of telocentrics, a process to be described later. 


With 6-3 per cent. of all P.G’s. having fragments it is not surprising 
that a proportion of all progeny from selfed monosomics have either 
telocentrics or isochromosomes. Deficiencies can be transmitted 
through the egg as well. Of chief importance to the plant breeder 
is the phenotypic expression of small or large deficiencies or duplications 
which can arise through the processes of misdivision of the univalent 
chromosome at meiosis. 


8. SUMMARY 


1. The 21 expected monosomics in 7. vulgare can be identified 
by crossing. Only two or three of them can be distinguished by 
chromosome structure. 

2. In all monosomics 20 bivalents are formed in 97 per cent. of 
the P.M.C’s. Chromosome mosaics, both hypoploid and hyperploid 
cells, were found in monosomic and disomic lines (see Appendix). 

3. The univalents either divide or misdivide and either lag or do 
not lag at either first or second telophase. 

4. Isochromosomes in wheat arise directly through misdivision 
at meiosis and not indirectly through non-disjunction in the P.G’s. 

5. P.G’s. with 18-23 chromosomes are formed and undergo mitosis. 
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Six per cent. of the P.G’s. contain fragments, most of which are 
telocentrics. 

6. Over twice as many P.G’s. with 20 as with 21 chromosomes 
are produced but the former are only 1/30th as effective as the latter 
in fertilisation. 


APPENDIX 


Chromosome mosaics at metaphase I.—At MI 32 cells with a reduced or aneuploid 
number and two cells with a hyperploid number were observed (table 9). Following 
the terminology of Sachs (1952) who found aneuploid cells in Triticum hybrids, 
these abnormalities are classed as chromosome mosaics. 

Frequencies are not intended to show differences among the monosomes because 
in other slides from the same spikes no mosaics were found. Normally there are 
over 500 P.M.C’s. in a squash of three anthers but because of non-synchronisation 
only about 300 at full MI are expected. In those slides where the number of cells 
at MI is small the predominating stage was telophase or interphase. Presumably 
the unbalanced cells are less efficient and thus are delayed. 




















TABLE 9 
Distribution of chromosome mosaics in pollen mother cells 
Chromosome mosaics 
Total no. of Sa Geen ——— al 
Plant type egal g a No. of Chromosome we 
ganthers | aberrant pr thon the Pairing 
“ells aberrant arrangement | 
a cells | 
Disomics 
A 2 : ‘ 36 I 29 84 131 
B ~ : — 43 2 29 g4 13 
MonosomIics 
I e Ps or 198 I 77 tT 374 3} | 
i ae P wt 83 | 6 30 10! rot | 
VE. : , 222 3 2 12! gt 
| {1 23 5" 131 
Vig. se] 96 | li 1a he 
XV . : ‘ 200 2 19 64 71 
| uu 1 
XVII : . 73 14 32 {ie be | 
XVIII ‘ : 1021 * | 2 32 122 gt | 
Monosomic | 
XIV+ISO—XIV 93 I 73 + 357 3 | 
* Count from 3 florets. t+ Hyperploid cells. 


In mono-XVII evidently the abnormality must have occurred at least four 
cell generations prior to meiosis. Premeiotic mitoses were examined in two plants 
but out of 275 cells at anaphase or metaphase no spindle irregularities were observed. 
There were three doubtful cells, two with a bridge and one with a micronucleus. 
With abnormal divisions occurring so early in the anther development it is difficult 
to get cytological proof. It is probable, however, that the aneuploid cells are caused 
by multipolar spindles as suggested by Sachs (1952). Hyperploid cells can be 
explained by fusion of two spindles. 
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Because the chromosome mosaics appeared in disomic plants it may be an 
indication of the instability of Chinese Spring. It may also suggest that the chromo- 
some mosaics may be the result of either environmental stress or gene mutation 
and that close scrutiny would reveal their presence in other non-hybrid material. 
With the rare occurrence of chromosome mosaics and with the competition they 
would suffer it is unlikely that many aneuploid cells would form viable P.G’s. 
However, it is a possibility. They thus may prove a source of variation for plant 
breeders. In any event they offer an explanation for some of the abnormal chromo- 
some numbers reported from time to time in cereals. 
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Fic. 4.—Photomicrographs of meiotic stages and mitosis in the P.G’s. D and E at x 500, 
all others at x 1000. 


A. 
B. 
C. 


D. 
E. 


MI in P.M.C., 20% 11. 
MI in aneuploid P.M.C. from mono-VIII, 5" 13}. 


AIl, one side only, from mono-III, 20 chromosomes towards each pole with 
the misdividing univalent between them. 


MI in M.M.C., 207 13, 


Late TII in M.M.C. No laggards and therefore no micronuclei formed in the 
megaspores. 


P.G. mitosis, 21 chromosomes. 


P.G. mitosis, 20 plus a telocentric chromosome, in the group at two o’clock. 


. P.G. mitosis, 20 chromosomes. 


P.G. mitosis, anaphase from mono-VII, sub-efficient chromosome at six o’clock 
left. 20 chromosomes at each pole. 
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1. INTRODUCTION 


GERM CELLS with unreduced nuclei arise in plants and animals from 
irregularities in meiosis, owing to genotypic causes or to structural 
hybridity. Many types of upset in the meiotic cycle have been 
reported and we are now coming to know the kinds of errors by which 
these variations arise. 

Genotypically controlled irregularities are often associated with 
reduced precocity of the early meiotic stages. A newtypeofirregularity, 
intermediate between normal meiotic and completely asynaptic 
behaviour, has now been found in Chrysanthemum atratum. 

C. atratum exists in two forms, diploid 2n = 18 and autotetraploid 
2n = 36. All nine chromosomes of the basic set have nearly median 
centromeres. Meiosis in the diploid form is normal. The meiotic 
abnormalities were found in four tetraploid plants, probably of one 
clone, obtained from the Royal Horticultural Society, Wisley. The 
similarity in behaviour of these plants suggests that they have arisen 
from a clone.* 


2. MEIOSIS 


Normal cells ——Pairing and chiasma formation in normal cells is 
like that in other autotetraploid species of Chrysanthemum. Associations 
of two, three and four chromosomes are observed with varying numbers 
of univalents ; bivalents each had 1-3 chiasmata (fig. 1). Both first 
and second meiotic divisions follow the usual course, the univalents 
dividing at AI. There is the expected variation in chromosome 
number of the resultant pollen grains. 

Precentric cells—In 63 per cent. of pollen mother cells there were 
regularly 18 bivalents at MI. These were abnormal in that they 
showed a precocious splitting of their centromeres as well as a reduced 
chiasma frequency. In fact each had a single terminal chiasma. 
Cells behaving in this unusual way have been termed “ precentric 
cells’ and the chromosome associations they contain “ precentric 
bivalents.” Seven per cent. of cells were intermediate and had both 

* A modified Feulgen technique was used for anther squashes with acetic alcohol 
(1 : 3) as fixative. A hydrolysis time of 4 mins. was sufficient as opposed to the 15 mins. 
for somatic root-tip squashes (Dowrick, 1952). Anthers were also fixed in a solution of 
absolute alcohol and saturated ferric acetate in glacial acetic acid (3:1) and stained in 


acetocarmine. Extra iron was unnecessary. Both methods gave satisfactory results. 
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normal and precentric chromosome associations in the same cell 


(table, fig. 3). 


Fic. 1.—(a) First metaphase in a normal pollen mother cell of C. atratum. (6) First 
anaphase in a normal cell. x 1200. 


Abnormalities are confined to pollen mother cells: root-tip 
mitoses are normal. The abnormal divisions were present in similar 
proportions not only in all four plants, but also in all five anthers of 

TABLE 
Frequency of precentric pollen mother cells at MI 




















| Rif d Normal Mixed | Precentric 
ee : | es Sees D cncasiiegi 
Plant . oI 26 3 62 
Plant 2 . 2 79 23 | 7 | 49 
Plant 3 . ; . 52 18 4 30 
Plant 4 . : ms) 107 3I | 10 | 66 
aa ‘ | 329 98 | 24 207 
a hetactint | = ; 
Percentage. ‘ nae 30 7 | 63 
the individual florets (table). There is a random distribution 


of precentric cells. They are scattered throughout the anthers, both 
in groups and isolated. 
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The division of the centromeres does not occur until the bivalents 
are already on the metaphase plate. The arrangement of the bivalents 
is then normal. Their appearance, with divided centromeres and 


completely terminalised chiasmata, is quite characteristic (fig. 2 and 
Plate). 


ny hs Yt ne 


Fic. 2.—(u) and (4) first metaphase in precentric cells. (c) First anaphase with precociously 
split chromosomes resulting from such metaphases. X 1200. 


The half-bivalents separate from each other normally, each 
passing to its respective pole without active repulsion between the 
divided centromeres of each chromosome. However, having already 
divided precociously, the centromeres of the chromatids tend, during 
anaphase movement, to fall apart with the result that at first telophase 
each nucleus contains 36 chromosomes. The first meiotic division 
of the abnormal cells does not, therefore, reduce the chromosome 
number. Nor, since crossing over is near the ends, can it greatly 
change the genetic constitution of the cell. 

In normal cells at the conclusion of first division the chromosomes 
are, with the exception of the centromere, structurally double. But, 
as they have single centromeres, they are mechanically single. 
Consequently, at MII, when the centromeres divide and the two 
halves of each chromosome separate, the number of centromeres in 
each daughter cell is half that in the mother cell. 

In precentric pollen mother cells, however, where the chromosomes 
have divided at MI and the half chromosomes have separated at AI, 
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the interphase chromosome threads are single throughout. They 
are similar to daughter univalents and, like them, do not divide 
again during a second meiotic division (Darlington, 1939). Where 
this applies to all the chromosomes it appears that the second division 
of the nucleus itself is precluded. This observation lends support to 
the view that the occurrence of a second meiotic division is determined, 
or precipitated, by the undivided nature of the centromere at the 
end of first division. Otherwise it might be expected that a second 
division would be found in the present case. In fact meiosis ends 
with the formation of dyads, the cells of which retain the somatic 
chromosome number. Both diploid and tetraploid pollen grains are 
therefore found in the same anther. 


h- I? 
WahC 


Fic, 3.—First metaphase in cell with both normal and precentric chromosome associations. 
x 1600. 


In the 7 per cent. of cells with mixed chromosome associations 
their proportions vary (fig. 3). I have been unable to trace their 
later history but the interphase nuclei arising from them will contain 
both single and double chromosome threads and, therefore, no doubt, 
be capable of undergoing a second division. The daughter chromo- 
somes from the precentric bivalents will, presumably, be scattered 
through the four resultant nuclei. 

In 3 of the 207 precentric cells (table), there were bivalents 
with chiasmata at both ends. Their behaviour is similar to those 
with a single chiasma. 

In normal cells, where pairing and coiling are complete, the 
chiasmata are situated near the chromosome ends. Pairing must, 
therefore, start at the ends if, as has been suggested previously by 
Darlington (1935), the degree and distribution of pachytene pairing 
and coiling determine the number and position of chiasmata. It 
must then be assumed that, in the precentric cells of C. atratum, the 
paired chromosomes usually pair and coil at only one end and 
consequently only a single terminal chiasma is formed. 
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3. THE ORIGIN OF THE ABNORMALITY 


These peculiarities of meiosis arise, as we should expect, before 
first metaphase. The precentric cells come into prophase later, for 
while normal cells are at diplotene the abnormal cells are still in a 
pre-pachytene condition. Later they catch up and both types of 
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Fic. 4.—The two types of meiosis in pollen mother cells of C. atratum. The reduced pairing 
and precociously split centromeres in precentric cells results in the formation of dyads. 
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cell enter first metaphase together (fig. 4). Pachytene and later 
prophase stages must, therefore, progress more rapidly than in normal 
cells. Chromosome pairing and relational coiling in the precentric 
cells is, it is assumed, lesscomplete than in the normal cells. This 
will result in a reduced chiasma frequency and, if coiling starts at 
the chromosome ends, more complete terminalisation. It will not, 
however, explain the precocious division of the centromeres. 





















































wm 30F RM.CELLS 
Z he ae | 
© 20} 
° RGRAINS 
$ (OF 4 
2 40} b 
| L. < 
‘ 12 s § 
ire 
a RM. NUCLEI ° 
P| ° 
S2 2 20} 
z 
w 
Oy 
fe) 
z 
6 10 12 14 16 18 20 +2 
RADIUS MICRONS 


Fic. 5.—The relative sizes of normal and precentric cells and nuclei. 


Correlated with these timing abnormalities are significant size 
differences between the different cell types. At early prophase, 
pollen mother cells which will show an irregular meiosis are more 
than twice as large as normal cells while their nuclei are approximately 
four times as large (fig. 5). This size advantage of the precentric 
cells is retained throughout the first division. Hence mature pollen 
grains arising from dyads have a volume, not twice, but more than 
three times as great as normal grains. 


4. DISCUSSION 


The abnormalities in C. atratum are, presumably, brought on by 
the genetic control of physiological differences in the cells at the 
onset of meiosis. 

C. atratum and the parthenogenetic earthworms described by 
Muldal (1952), form a special group, in regard to abnormal meiosis, 
where there is precocious splitting of the centromeres and the 
production of unreduced gametes. This precentric behaviour is 
intermediate between normal meiosis and “ asynapsis”’ where 
chiasmata are absent and only univalents are found. In all other 
respects the abnormality in C. atratum is analagous to that found in 
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Antennaria (Stebbins, 1932), Hieracium (Gustafsson, 1942), Calamagrostis 
(Nygren, 1946) and Scilla (Rees, 1952). 


The similarities are :— 


(i) Only some germ cells, not all, in each anther are affected. 
The proportion of abnormal cells varies from as little as 
I per cent. in some samples of Scilla to 63 per cent. in 
C. atratum and gg per cent. in Datura. 

(ii) The distribution of abnormal cells is always at random. 

(iii) Prophase begins later in abnormal cells, metaphase some- 
times, but not always, catching up. 

(iv) There is increase in the size of the retarded cells in relation 
to the normal. 


Male sterile Lathyrus odoratus (Upcott, 1937) is, like C. atratum, 
also intermediate between asynapsis and normal meiosis. It is also 
similar in that prophase in abnormal cells begins later. It differs 
in being a total defect of all cells and in not being precentric. The 
delay in early meiosis brings about, in this case, a greater terminal 
localisation of the chiasmata at first metaphase. Carried to its 
extreme this increased localisation would result, as in Scilla, in 
asynapsis. 

The behaviour of E.M.C’s. in these various plants, may be similar 
to, or different from, the pollen mother cells. The cells arising from 
the abnormal divisions may either give rise to fertile gametes or, in 
other cases, degenerate at varying stages before reaching maturity. 
In Antennaria, Hieracium and Calamagrostis, meiosis may follow yet a 
third course : it may be replaced by mitosis in the abnormal pollen 
mother cells and embryo sacs. 

Many other examples of asynapsis are known at meiosis in plants 
where the properties that are discussed here are not described, as 
in Datura (Satina and Blakeslee, 1934) and Pyrus (Thomas, 1942). 


5. SUMMARY 


1. An autotetraploid clone of C. atratum has unreduced pollen 
as a result of abnormal meiosis in 68 per cent. of pollen mother cells. 

2. The abnormal cells lag at prophase, have a prolonged growth 
phase, and are less precocious in chromosome development than 
normal cells. They are consequently larger in both cell and nuclear 
size throughout all stages of division. 

3. In these cells no quadrivalents are formed. The bivalents 
show a reduced chiasma frequency and precocious splitting of the 
centromeres at MI. The chromatids fall apart during anaphase. 

4. The first division gives rise to cells with 36 chromosomes. 
A second division is omitted. Thus, it appears, it is the capacity 
of the centromeres for division which determines or precipitates the 
occurrence of the second meiotic division. 
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Plate 
. Normal first metaphase with four quadrivalents and ten bivalents. x 2600. 


. First metaphase in precentric cell with eighteen bivalents and precociously divided 
centromeres. X 2600. 


. Two normal first anaphases with eighteen daughter bivalents in each. X 1000. 


. First anaphase in a precentric cell ; 36 chromosomes passing to each pole. X 1000. 
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1, INTRODUCTION 


THE study of mammalian chromosomes, despite their importance, 
has been greatly neglected. This is in part due to the idea that they 
are unsuitable for a proper cytological study (¢g. White, 1949). 
The best mammalian chromosomes are to be found in the marsupials, 
many of which have a small number of large chromosomes. Most 
other mammals on the other hand, have a large number of small 
chromosomes so that an accurate study is certainly a more difficult 
task than in many other organisms. But the extreme difficulties 
which have often been attributed to the study of mammalian chromo- 
somes are also in part due to a lack of technique. This has led to 
confusion on even such matters as the chromosome number of a 
particular species (Sachs, 1952b) and it now seems desirable to 
re-investigate much of the earlier work using more suitable 
techniques. 

One of the many interesting points in mammals is that the animals 
have sex chromosomes which may be morphologically distinguished 
from the autosomes. These sex chromosomes are of importance both 
from a genetical and from a cytological point of view. Their genetical 
interest is connected with the existence of complete and partial sex 
linkage and their cytological interest is connected with their special 
behaviour during mitosis and meiosis. Some species of mammals 
show little or no size differences between the sex chromosomes and 
the autosomes. In these cases, therefore, identification of the sex 
chromosomes has not been possible at all the stages of division. But 
in Microtus agrestis, the short-tailed bank vole, the sex chromosomes 
are enormous and there is no difficulty in distinguishing them at 
mitosis and meiosis. An investigation of these giant sex chromosomes 
has been one of the purposes of the present study. 

The existence of giant sex chromosomes in Swiss animals of M. 
agrestis had previously been reported by Matthey (1950) whose study 
was almost entirely based upon sectioned material. It has, however, 
now become possible, by the use of suitable squash techniques, to 
obtain much additional information not only on the sex chromosomes 
but also on general problems of chromosome behaviour and the 
sequence of meiosis. 
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2. MATERIAL AND TECHNIQUES 


The animals used in the present study were derived from two sources. Some 
were caught on the grounds of Bayfordbury, and others were obtained from Dr D. 
Chitty of Oxford. Animals from the two locations gave the same results. 

The main technique was fixation in 1 : 3 acetic-alcohol followed by Feulgen 
staining and squashing (Sachs, 19526) and this gave uniformly good results. In 
addition, fixation in aceto-carmine or acetic acid followed by the normal method 
of Feulgen staining, and fixation in 1 : 3 acetic-alcohol with a small amount of 
ferric chloride followed by aceto-carmine staining, also proved useful. As is shown 
by their own photographs, the staining procedures for mammalian testes recom- 
mended by Slizynski (1949) and Makino and Nishimura (1952) show no observable 
advantage over the normal method of making Feulgen squashes. This also applies 
to the section cutting method which is used by Matthey (1950). 


3. HETEROCHROMATIN AND THE SEX CHROMOSOMES 


The English animals of M. agrestis from Hertfordshire and Oxford- 
shire, like the Swiss animals studied by Matthey (1950) have a 
chromosome number of 2n = 50. This number and the giant sex 
chromosomes thus appear to be characteristic for the species. The 
giant sex chromosomes can be easily identified in spermatogonial 
mitosis. The size of each chromosome, however, varies according 
to stage and fixation, so that it is useless to make absolute chromosome 
measurements. At mitotic metaphase, however, after acetic-alcohol 
fixation and Feulgen staining, the autosomes are about 1-2" long 
whereas the X chromosome is about 9’ and the Y chromosome 
about 6” long (plate I, fig. 5). 

Spermatogonial cells of the testes show that the sex chromosomes 
contain heterochromatin and this can be seen in some resting nuclei 
as two distinct chromocenters. But other resting nuclei, in different 
tissues, show no chromocenters. The detection of heterochromatin 
in the sex chromosomes by their appearance as chromocenters is 
thus not possible in all types of cells. 

This presence or absence of chromocenters in the so-called resting 
nucleus is probably due to the different metabolic activities of the 
cells concerned. That the cytological appearance of heterochromatin 
can be modified by different metabolic activities of the cell was first 
found in studies on the effects of low temperature (Darlington and 
La Cour, 1938, 1940, 1941 ; La Cour 1951, e¢ al.). This has shown 
that heterochromatin can, in certain species, show different differential 
behaviour under different environmental conditions. Heterochromatin 
which was intensely stained in resting nuclei could, after the effects 
of low temperature, be less intensely stained than the euchromatin 
at metaphase. Thus, as Darlington (1947) has pointed out, this 
behaviour of heterochromatin can be explained by the presence of 
genes which show allocycly under certain conditions. It can now 
be seen from the giant sex chromosomes that under some conditions 
the heterochromatin completely loses its differential stainability even 
in the resting nucleus. 
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The allocyclic behaviour of heterochromatin, including the loss 
of differential stainability, can also be seen in the behaviour of the 
X chromosomes in some grasshoppers and crickets (White, 1940). 
In these insects the X chromosome shows what White calls “‘ negative 
heteropycnosis ” during the early spermatogonial divisions ; later it 
loses its differential stainability ; and it becomes “ positively hetero- 
pycnotic ” during the prophase of meiosis. This behaviour of the 
X chromosome can also be related to the different appearance of 
heterochromatin under the differing conditions of the cells. Thus, 
in some cells, the heterochromatic genes react in the same way as 
the euchromatic genes and cannot be cytologically detected. In 
other cells, with the appropriate metabolic conditions, they can be 
detected by their differential metabolic cycle and a collection of such 
neighbouring genes can then be seen as distinct blocks of hetero- 
chromatin. So it is that, although some resting nuclei have no 
chromocenters, the appropriate conditions show that the X and the 
Y in M. agrestis contain heterochromatin. 


4. THE SEX CHROMOSOMES AT MITOSIS 


In spermatogonial cells the sex chromosomes at mitosis are, in 
relation to the autosomes, precocious during prophase and they lag 
during anaphase. But at mitotic metaphase both the sex chromosomes 
and the autosomes are arranged upon the spindle. During anaphase 
the X and the Y start to separate from the centromeres, but they are 
still held together when the autosomes have completely separated 
to the two poles (plate I, figs. 4 and 6). This slower separation of 
both X and Y may be in part due to the large size of these chromo- 
somes moving in a restricted space. But it is more probably also due 
to the heterochromatin holding together the two chromatids. 

In spermatogonial cells the spiral structure of the sex chromosomes 
is especially clear during the early prophase in the pre-meiotic mitoses. 
What appears to be an occasional case of somatic pairing, between 
a small segment of the X and a small segment of the Y, has also been 
seen at this stage (plate I, fig. 2). 


5. TERMINAL CENTROMERES 


The mitotic cells in the ovary show that the female also have 
two enormous sex chromosomes, and that they both have submedian 
centromeres. The giant chromosomes with submedian centromeres 
are, therefore, the X chromosomes. The position of the centromere 
in the X and Y of the male has been most clearly seen at spermato- 
gonial mitosis (plate I, figs. 3 and 5). The anaphase figures in these 
cells confirms that the X has a submedian centromere whereas the 
centromere of the Y appears to be terminal. Examination at mitotic 
metaphase again confirms the submedian centromere of the X but 
it also shows that the Y has, in addition to the long arm, a small 
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structure at the centromere end of the chromosome. This does not, 
however, necessarily imply (¢f. White, 1951) that telocentrics in 
reality have a small second arm. 

The existence of stable telocentrics, i.e. chromosomes with stable 
terminal centromeres, has often been doubted (Matthey, 19514; 
White, 1951). This is partly due to the fact that the newly-formed 
telocentrics studies in plants such as Zea, Tulipa and Fritillaria, were 
found to be unstable. But this unstable condition does not apply to 
all telocentrics. Thus Darlington and La Cour (1950) have found 
that telocentrics in plants of Campanula may be unstable or stable, 
and that stable telocentrics can maintain themselves in nature. A 
further point which is often overlooked by those who doubt the 
existence of stable telocentrics, is that the centromere is not just a 
gap or a single point in the chromosome. The centromere is a 
multiple structure which can be broken at different places (Darlington, 
1939). The functional telocentric will thus still possess at least some 
and probably all of this centromere structure so that the centromere 
end of the telocentric is not just an abrupt end to the chromosome arm. 
The existence of this centromere structure in telocentrics has in fact 
been seen in plants of Phleum (Ellerstr6m and Tjio, 1950). The Y 
in M. agrestis may also be a telocentric showing this centromere 
structure at metaphase. It can thus be concluded that stable telo- 
centrics can exist, and their existence must be reckoned as a factor 
in chromosome evolution. 


6. A SEX VESICLE AT MEIOSIS 


An examination of testes containing cells at meiosis shows that 
the autosomes in many of the cells are at pachytene. A most striking 
first impression of these cells is that the nucleus contains, in addition 
to the pachytene autosomes, a large oblong body which is strongly 
Feulgen positive. Closer examination shows that this large Feulgen 
positive body consists of two Feulgen positive chromosomes, the X 
and the Y, surrounded by a somewhat lighter stained vesicle (fig. 10, 
and plate II, fig. 4). Thus when the autosomes are at pachytene, 
the X and the Y are included in a special vesicle which can be 
designated as a sex vesicle. The existence of a sex vesicle surrounding 
the giant sex chromosomes in another species of Microtus has been 
reported by Muldal (1949). 

The development of this sex vesicle in M. agrestis has already 
begun when the autosomes are at leptotene. At this stage the sex 
chromosomes are still separate, the X and the Y each forming what 
appears to be its own vesicle. The X and the Y then move together 
until they are included in the one large vesicle at pachytene (plate II, 
figs. 1, 2 and 4). The vesicle surrounding the sex chromosomes still 
persists at diplotene (fig. 1A and B) and diakinesis, but it has in most 
cases disappeared when the autosomes reach 1st metaphase. The 
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formation of this sex vesicle is characteristic of meiosis and it has 
never been seen at mitosis. 


The sex chromosomes are thus inside their vesicle during those 
stages of meiosis where chiasmata are normally formed. This raises 
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Aand B. The X and Y inside the sex vesicle x 2000 at diplotene. 

C. The X and Y inside the sex vesicle at pachytene. 

Dand E. The X and Y connected at 1st metaphase. 

F, The X and Y not connected at early 1st anaphase. 

G, Hand J. The X and Y connected and lagging at late 1st anaphase. 


N.B.—The arrangement of the sex chromosomes inside the sex vesicle is comparable to 
their arrangement when the sex vesicle has disappeared. 
In D, E and F autosomes and spindle are drawn diagrammatically. 


the question whether chiasma formation between the X and the Y 
in M. agrestis can and actually does take place. Matthey (1950), 
on the basis of his observations at 1st metaphase, at first concluded 
that chiasmata had been formed between the X and Y in M. agrestis. 
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Later, however (1951a), he concluded that this was doubtful. Both 
these opposing conclusions were based on the same set of observations. 
This confusion is probably caused by Matthey’s own observations 
on what he supposed to be the meiotic development of the sex chromo- 
somes before rst metaphase. These had led him to conclude that the 
X and the Y flow together into a large homogeneous mass of material. 
He considers this meiotic development of the sex chromosomes as 
analogous to the fusion of two drops of mercury. The union of two 
masses of material to give one large homogeneous mass at pachytene 
thus seemed a contradiction to any normal pachytene pairing and 
chiasma formation between the sex chromosomes. 

By the use of suitable techniques I find, however, that the early 
meiotic development of the sex chromosomes is not at all analogous 
to the fusion of two drops of mercury to form a homogeneous mass 
at pachytene. It can be seen in suitable preparations that the oblong 
Feulgen-positive body at pachytene is not a homogeneous structure, 
that it consists of two sex chromosomes inside a sex vesicle, and that 
there is no inherent contradiction between the early meiotic develop- 
ment of the sex chromosomes and the configurations seen at Ist 
metaphase (fig. 1). 


7. THE FORMATION OF CHIASMATA 


A study of pachytene shows that the X and Y are arranged in 
characteristic ways inside the sex vesicle. Owing to the faint differentia- 
tion between the materials of the sex chromosomes and of the vesicle, 
it is difficult to be sure of the exact mode of pairing. It appears, 
however, that X and Y are connected inside the vesicle end to end or 
by a small segment. ‘This segment could be about the same size as 
the small somatically paired segment between the X and Y which is 
occasionally seen at the prophase of mitosis. 

The X chromosome lies curved inside the vesicle at pachytene 
and diplotene and it maintains this curved shape at Ist metaphase 
and even 1st anaphase. The Y which usually lies in a straighter 
position inside the vesicle also maintains its straighter position at 
ist metaphase and anaphase. The sex vesicle often lies on one side 
of the nucleus and the sex chromosomes at Ist metaphase are often 
on the edge of the plate. 

The pachytene and 1st metaphase configurations can both be 
explained by the existence of a small pairing segment between the 
X and the Y. The configurations at 1st metaphase and 1st anaphase 
show that this pairing segment would be at the end furthest from the 
centromere in both chromosomes, and (as is to be expected) only 
reductional separation of the sex chromosomes occurs at Ist anaphase. 

The X and Y chromosomes are not always connected at Ist meta- 
phase and in these cells the sex chromosomes have presumably been 
unpaired. 
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Darlington (1934) was the first to suggest a cytological explanation 
which could account for complete and incomplete sex linkage in 
mammals. This showed that both types of sex linkage could be 
explained by the existence of a pairing segment and a differential 
(non-pairing) segment in the X and the Y. Genes on the differential 
segment would be completely sex-linked. Genes on the pairing 
segment would be partially sex-linked by the formation of chiasmata 
within the pairing segments. In the absence of any knowledge of 
the genetics of M. agrestis, nothing can be said about the presence or 
absence of incomplete sex-linkage in this animal. But from the 
cytological point of view the pachytene, 1st metaphase and 1st anaphase 
configurations can all be explained if the sex chromosomes have a 
small pairing and large differential segments. Chiasmata can then 
be formed within the pairing segment of the X and the Y. Remnants 
of the sex vesicle may, however, still be present after diakinesis. ‘This 
could conceal the chiasmata at Ist metaphase and cause the sex 
chromosomes to lag at 1st anaphase. 

The presence of a sex vesicle during the prophase of meiosis which 
is so prominently seen in M. agrestis, is not a unique phenomenon 
restricted to this species. I have seen similarities to this type of meiosis 
in a number of other mammals. As genes for sex determination are 
presumably located on the sex chromosomes, the formation of a sex 
vesicle during the prophase of meiosis could have distinct evolutionary 
advantages. The sex vesicle prevents any possibility of crossing-over 
between the sex chromosomes and the autosomes which might other- 
wise arise owing to errors of pairing at pachytene. The sex vesicle 
could thus ensure that these sex-determining genes are always main- 
tained on the sex chromosomes. This mechanism could then prevent 
the occurrence of gene combinations which might upset the normal 
balance for sex determination. 


8. ORIGIN OF THE GIANT CHROMOSOMES 


There is very little cytological data on the chromosomes of other 
species in the genus Microtus. Out of more than 200 different forms, 
only about half a dozen have been cytologically examined. Even 
among these there is some disagreement as to the chromosome number 
and the presence of giant sex chromosomes. It is, however, clear 
that species with higher chromosome numbers than M. agrestis have 
been found in Microtus and in the allied genus Clethrionomys. ‘Thus 
M. nivalis (Matthey, 19515) and C. glareolus (Sachs, 1952b) both have 
2n = 56 and no giant sex chromosomes. 

M. agrestis with 2n = 50 and giant sex chromosomes could have 
originated from these or similar species with higher chromosome 
numbers. Thus in the plant Phyllospadix japonica, Harada (1944) 
reports that the female plants have 20 small chromosomes and the 
male plants have 15 small and one enormous chromosome. This 
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large chromosome contains heterochromatin, is designated as the 
Y chromosome, and appears to have originated from the combination 
of 5 short chromosomes. 


9, LAMPBRUSH CHROMOSOMES AND THE 
SEQUENCE OF MEIOSIS 


Meiosis in the male of M. agrestis includes a modified form of 
what appears to be diplotene. At this stage the autosomes have a 
*‘ woolly ” outline with small protrusions connected to the chromo- 
somes (plate II, fig. 3). This type of diplotene, described as diffuse 
diplotene, has been observed in a number of plants and animals 
(Darlington, 1937). Its most extreme development has been found 
in animal oocytes, where the chromosomes are referred to as lamp- 
brush chromosomes. There is a distinct similarity between these 
and what can be called the lampbrush chromosomes in M. agrestis. 
In M. agrestis the lampbrush structures begin to develop at pachytene 
and show their maximum development at what is apparently diplotene. 
In the animal oocytes there is also a maximum development of the 
lampbrush at what is apparently diplotene (Dodson, 1948). 

It can be shown that lampbrush development is associated with 
chemical changes in the chromosomes. For example, in the oocytes 
development of the lampbrush chromosomes is associated with a 
decrease in Feulgen staining. ‘This decrease is associated with an 
increase in length of the chromosomes and in fully grown oocytes 
the chromosomes, except in a few small regions, are completely 
Feulgen negative (Serra, 1947). The change from a Feulgen-positive 
to a Feulgen-negative chromosome associated with this extreme 
lampbrush development may be due (a) to the increased growth of 
the chromosomes with the same amount of DNA (desoxyribose nucleic 
acid) decreasing the intensity of Feulgen staining and/or (4) to the 
transfer of DNA from the chromosomes. 

Ris (1951) maintains the former explanation and denies the 
existence of DNA transfer from the chromosomes. There are, however, 
several lines of evidence which show that nucleic acid transfer does 
take place. Thus, the lampbrush development in M. agrestis is not 
associated with an increased size of the chromosomes. In this species 
the chromosomes at diplotene are if anything shorter than those at 
preceding stages of meiosis. ‘The same applies to the lampbrush 
stages in grasshoppers (Srivastava, 1951). 

Evidence for the existence of nucleic acid transfer can also be 
obtained from studies on mitosis. Howard and Pelc (1951) have 
concluded that DNA is synthesised during the resting stage in cells 
which are going to divide. Distinct Feulgen-positive chromosomes 
are however only found at prophase, so that it is reasonable to assume 
that the DNA synthesised during the resting stage is transferred to 
the chromosomes at prophase. In the same way DNA can be trans- 
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ferred from the chromosomes at telophase when the chromosomes 
lose their distinct Feulgen positive structure. The transfer of DNA 
to the chromosomes can also be inferred from results on the effects 
of low temperature (Darlington, 1947). 

Studies on RNA (ribose nucleic acid) have likewise shown the 
transfer of RNA to and from the chromosomes. Thus Jacobson and 
Webb (1952), who assume that nucleic acid in the cell always occurs 
in conjunction with a protein, conclude that ribonucleoprotein is 
transferred to the chromosomes during prophase and transferred 
from the chromosomes during anaphase. This latter process can be 
inferred by the presence of the transferred ribonucleoprotein between 
the two groups of anaphase chromosomes. These different lines of 
approach all show that nucleic acid transfer to and from the chromo- 
somes can and does take place. 

The synthetic processes of mitotic cells, such as protein synthesis, 
take place during the resting stage and it can be seen that they are 
preceded or associated with the transfer of nucleic acid from the 
chromosomes. Those stages of meiosis which show evidence of nucleic 
acid transfer from the chromosomes, as in the case in the formation 
of lampbrush chromosomes, are thus in this respect analogous to the 
synthetic part of the resting stage in mitosis. This analogy between 
lampbrush and resting stages can further be seen in the behaviour of 
heterochromatin. It has been shown that the sex chromosomes in 
M. agrestis contain heterochromatin which can be seen as chromo- 
centers in some resting nuclei. These parts of the chromosomes 
therefore can maintain their DNA content even in these resting stages. 
In the same way the sex vesicle at diplotene, which contains the 
sex chromosomes, does not show the lampbrush formation which is 
shown by the autosomes: the heterochromatin during this stage of 
meiosis maintains its DNA content. The same differential behaviour 
of heterochromatin has thus been found both in resting stages of 
mitosis and in lampbrush stages of meiosis, and this also emphasises 
the similarity of these two stages. Lampbrush formation during 
certain stages of meiosis may thus be considered as a sign of synthetic 
processes inside the cell. 

The existence of synthetic processes after chromosome pairing 
enables us to extend our knowledge of the sequence of meiosis. It 
has been shown in a number of organisms that the beginning of 
meiosis is connected with special circumstances. These include the 
supply of materials from tissues surrounding the germ cells, and the 
existence of special types of mitosis before meiosis. The supply of 
materials from tissues surrounding the germ cells (Darlington and 
La Cour, 1946; Painter, 1943; Montalenti et al., 1950) includes 
the possibility of nucleic acid secretion associated with the beginning 
of meiosis. The direct evidence for special types of mitosis before 
meiosis comes from M. agrestis and some insects (White, 1940) which 
show that meiosis is preceded by rapid mitotic cycles. In these 
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insects it has been seen that there is almost no resting stage between 
the pre-meiotic mitoses. There is thus not a complete DNA transfer 
from the chromosomes during telophase so that the autosomes have 
an irregular outline at the beginning of the next prophase. This 
irregular outline is also found in the autosomes of M. agrestis, but it is 
not shown by the sex chromosomes which can maintain their DNA 
content even during resting stages. The special character of pre- 
meiotic mitoses can also be inferred from the existence of mitotic 
abnormalities which only occur before meiosis (Sachs, 1952a). 

The sequence of events leading up to meiosis can explain the 
presence of lampbrush chromosomes. It can reasonably be assumed 
that the rapidity of mitosis before meiosis excludes synthetic processes 
which would also be sufficient for chromosome reproduction during 
meiosis. The supply of substances from outside the germ cells may 
allow these cells to start on meiosis and proceed up to pachytene. 
Judging from the large proportion of cells at pachytene in the spermato- 
cytes of M. agrestis, the cells remain at this stage for a long time. The 
continuation of meiosis is however dependent upon the synthesis of 
new chromosome material: this would take place after chromosome 
pairing and be made apparent by the lampbrush chromosomes. 

The existence of this synthetic stage after chromosome pairing at 
meiosis allows us to understand one of the fundamental differences 
between mitosis and meiosis. This difference lies in the time of 
chromosome duplication (Darlington, 1937). Chromosome duplica- 
tion at mitosis takes place at the resting stage and chromosome 
duplication at meiosis takes place after chromosome pairing at 
pachytene. It can now be seen that the resting stage of mitosis and 
these lampbrush stages after pairing at meiosis can both be considered 
as stages for the synthesis of chromosome material. The differences 
in the times of chromosome reproduction at mitosis and meiosis 
may thus be due to the differences in the time of synthesis of chromo- 
some material. 


10. SUMMARY 


1. The rodent Microtus agrestis, with 2n = 50, has enormous sex 
chromosomes which can be easily identified at mitosis and meiosis. 
The males have an X Y sex mechanism. 

2. Both X and Y contain (or entirely consist of) heterochromatin 
which appear as chromocenters in some cells. This variation in the 
appearance of heterochromatin is probably due to different metabolic 
activities in the different cells. 

3. At spermatogonial mitosis the sex chromosomes are precocious 
during prophase and lag during anaphase. ‘There is occasional 
somatic pairing at mitotic prophase between small almost terminal 
segments of X and Y. 

4. The X and Y are separate at the beginning of meiosis and are 
included in a special sex vesicle by pachytene. The arrangements 
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Plate I 


Microphotographs of mitosis in the testes of Microtus agrestis. Feulgen xX 3000. 


Fic. 1.—Early Prophase. Early development of the X and Y. 

Fic. 2.—Later prophase. Somatic pairing between a small segment of the X and Y. 
Fics. 3 and 5.—Metaphase. 

Fic. 4.—Early anaphase. X and Y beginning to separate. 

Fic. 6.—Late anaphase. Delayed separation of the X and Y. 
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Plate II 
Microphotographs of meiosis in the testes of Microtus agrestis. Feulgen Xx 2000. 


1.—Stages in the development of the sex vesicle. The X and Y still separate (left 
hand cell) ; lying next to one another (centre) ; nearly included in the one sex vesicle 
(right hand cell). 


2.—Sex vesicle at pachytene. Unsquashed nucleus. 

3.—Sex vesicle at diplotene. Lampbrush autosomes. 

4.—Sex vesicle at pachytene. Squashed nucleus. 

5.—Early anaphase I. Curved arrangement of the X and straighter position of the Y. 


6.—Late anaphase I. Lagging of X and Y. 
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of the sex chromosomes inside the sex vesicle are comparable to their 
arrangement at 1st metaphase when the vesicle has usually disappeared. 

5. The configurations of the sex chromosomes from pachytene 
to 1st anaphase can be explained by the existence of a small pairing 
and large differential segments. The pairing segments in both X and 
Y would then be at the ends furthest from the centromere. 

6. M. agrestis and its giant sex chromosomes could have originated 
from related species with higher chromosome numbers and no giant 
sex chromosomes. 

7. Meiosis in the male M. agrestis includes stages with lampbrush 
chromosomes showing similarities to resting stages between mitoses. 
The development of lampbrush chromosomes at meiosis is interpreted 
as a sign of synthesis of chromosome material. Differences in the 
time of chromosome duplication at mitosis and meiosis can thus be 
explained by differences in the time of synthesis of chromosome 
material. 
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1. INTRODUCTION 


Maize varieties in the field vary in their self-pollination. This 
indicates that control of the inbreeding-outbreeding mechanism is 
not so simple as is usually supposed. Usually a population’s mean 
outcrossing or selfing has been recorded, with neglect of plant to 
plant variation. Behaviour of individual plants within a population 
has not been considered. The variation has been reduced to a mean 
which conceals or rejects some of the most useful information. 

Maize is a monoecious (with male and female flowers separate on 
the same plant) outcrossing species with the possibility of inbreeding. 
The table gives the average inbreeding for different populations under 
different environments. These values, however, depend on similarities 
between the tested and tester variety for characters, such as relative 
flowering times, affecting dissemination and reception of pollen. 

Jones and Brooks (1950) have listed nine factors affecting natural 
outcrossing in cross-pollinated grasses. They have omitted protandry 
and protogyny, and also plant to plant variation in flowering time, 
presumably assuming these to be of little interest. Yet these properties 
must be of cardinal importance for the breeding system of maize. 
Even extensive contamination studies have failed to determine the 
amount of inbreeding and outbreeding between and within plants 
of one variety. 

Recently I have made a multiple regression analysis on silking 
time (protrusion of styles), tasselling time (opening of anthers with 
consequent shedding of pollen) and dichogamy, or the difference of 
the two in time, in order to determine their influences on out-crossing 
of sweet corn inter-planted with a flint variety. All effected con- 
tamination (Haskell, 1953). However, besides knowing behaviour 
of varieties as a whole, knowledge of individual behaviour within 
each is also valuable, as dichogamy relations of the breeding system 
are important in the evolution of species (cf. Darwin, 1888). A 
graphic method, suggested by Dr Darlington, is given here of rapidly 
evaluating individual behaviour within a population. 

Dichogamy consists of protandry and protogyny, each of which may 
be partial or absolute: absolute protandry is when stamens are 
ripe and all pollen has been shed before the stigmas are receptive, 
so that no self-pollination occurs within the flower ; partial protandry 
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is when much of the pollen is shed prior to the styles ripening, but 
there is some overlapping so that some of the earliest styles may 
receive pollen shed by late anthers. Protogyny is equivalent to 
negative protandry : with absolute protogyny the styles and stigmas 
are receptive and are already fertilised, or inactivated, before pollen 
starts to shed ; or some of the styles are still receptive when the pollen 
is being liberated (partial protogyny). 

There is a further complication in maize that the ears, really 
female inflorescences, do not open at the same time: the main ear 
usually opens ahead of the secondaries. In the flint types of both 
maize and sweet corn, with which we are here dealing, one ear alone 
predominates. In dent types there are several ears opening con- 
secutively, so that the amount of inbreeding or outbreeding is further 
complicated, being dependant on the position of the female inflorescence 
on the plant. 


2. EXPERIMENTAL DATA 


Professor K. Mather’s data were used. ‘These recorded first 
pollen shedding, and date silks first appeared on main ears of plants 
in two maize varieties, a flint and sweet corn, sown 24th May 1945. 
Fifty plants of each were raised at random in two adjacent plots. 
Data from both were pooled. There were altogether 94 flints and 
87 sweets, as a few plants died. 


TABLE 
Self-fertilisation in maize populations 





| 
| 
| 


| 
Comment Selfing Author 
|—_—____—— seh a ene ee eee mere) ae Peer 
| Dent corn in Nebraska. Both | 0-70 per cent. Kiesselbach, 1922 
| varieties with the same flowering | 
periods | 
; Dent corn in Maryland. Average f f 3°70 Be | Kemp and Rothgeb, 1943 
of five years | \ 5:00 re | Hayes, 1918 
| | 
| Maize crops generally throughout |{ 5:00 ‘i | Wallace and Bressman, 1949 
5 \ 5°13 rp | Waller, 1917 








Figures represent means of isolated plants surrounded by plants of another variety. 
Outcrossing was detected with marker genes, e.g. Y-y. 


3. THE TWO-DIMENSIONAL GRAPH 


Fig. 1 gives a conventional plant breeding method of estimating 
earliness and uniformity. The percentage of flints and sweets com- 
mencing to flower during July are given, with silking and tasselling 
plotted separately. The graph indicates broadly that flints silked 
considerably earlier than they tasselled ; the crop burst into pollen 
shedding. A few later flowering stragglers are depicted along the 
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bottom right of the graph. The sweet variety gradually came into 
flower, and was somewhat later than the flint. Several late silkers 
are also shown. 

This method has several disadvantages. Sunday measurements 
are absent, so introducing irregularities. More important, silking 
and tasselling relations of individuals within a variety are not shown. 
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DAYS AFTER FIRST JULY 


Fic. 1.—Conventional two-dimensional method of plotting earliness and uniformity. 
Open circles refer to sweet corn plants ; black circles refer to flinty maize. 


4. THE THREE-DIMENSIONAL GRAPH 


The same data are plotted in fig. 2. An arithmetical scale replaces 
actual flowering date in July. Silking date is plotted along the 
Y-axis, and tasselling date along the X-axis. Individual behaviour 
can now be charted ; as several plants were alike, circles of relative 
sizes are used to indicate how many have the same flowering relations. 
Percentages need not be calculated. Full and open circles respectively 
indicate flint and sweet corns. Circles on the dotted line represent 
plants silking and tasselling on the same day. Circles below the 
line are protogynous plants, while those above are protandrous ; 
distance from the line is proportional to the amount of protandry or 
protogyny. 

The breeding behaviour in both populations is now readily 
evaluated. The earliest flints are highly protogynous and represent 
a large proportion of the population. Slightly later ones are more 
homogamous. Latest ones tend to be more and more protandrous ; 
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however, some of the last flints to flower are also highly protogynous. 
The earliest sweet corn plants show little dichogamy but within the 
population there is a definite trend: later plants have a greater 
tendency to be protandrous. Only a few plants in the middle flowering 
period of the populations are protogynous. The best fitting straight 
line for the sweet corn data was calculated. This was Y = 1:16X —2:71, 
and is plotted in fig. 2. 
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Fic. 2.—Three-dimensional method of plotting earliness and uniformity. Open circles 
refer to sweet corn plants ; black circles refer to flinty maize. Size of circles is pro- 
portional to number of plants. Y = 1:16X—2°71 is the best fitting straight line for 
sweet corn values. 


5. FLOWERING TIME AND THE BREEDING SYSTEM 


Fig. 1 describes population behaviour without considering the 
individuals within it. It assumes that these are generally the same, 
except for their time of flowering. The breeding system of two 
populations in the same environment is also described in fig. 2, using a 


three-dimensional graph. This provides information not previously 
apparent. 
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The amount of selfing in a plant depends on availability of its 
own pollen. Hence earliest sweet corn plants have less chance of 
selfing, as they have available a pollen mixture before and after 
tasselling commences. Outbreeding is also more likely to be higher 
when plants are protogynous: this is mostly among earliest plants 
in both populations. Seeds from the earlier portion of the population, 
therefore, possibly include more heterozygotes than those from later 
plants. More vigorous offspring would thus come from earlier plants, 
because vigour in maize partly depends on heterozygosity. The 
stragglers are not only later, but pollen choice is limited. Being more 
protandrous, their own pollen is shedding when their silks emerge. 
Latest plants therefore have greater chances of selfing, especially 
as silks do not emerge simultaneously but often spread over several 
days. More of their seeds would then be inbred. Now inbred maize 
is generally later than the original open pollinated parent, hence the 
next generation from these seeds would be still later. Again they are 
more likely to be protandrous and to inbreed both within and between 
themselves. Thus, genetically, a vicious circle is set up, so that besides 
being temporally isolated from the bulk of most of the crop’s pollen, 
homozygosity is even further forced on them. 

Duration of pollen shedding, which has not yet been recorded in 
England, influences the inbreeding-outbreeding mechanism, so that 
absolute protandry does not actually occur in maize. Even though 
absolute protandry discourages self-pollination, it is still not fool-proof 
(Mather, 1940). The floral mechanism of maize has, however, arisen 
by artificial selection and is a particular condition : here, protandry 
has not been so important as sex separation within the plant. Duration 
of pollen shedding in the U.S.A. varies from two to fourteen days 
according to variety and most silks under favourable conditions are 
pollinated on day of emergence with ears completely pollinated 
within three days (Kiesselbach, 1949). Hence all degrees of inbreeding 
and outbreeding are possible. 

The genetical aspects of dichogamy in maize were earlier 
investigated by Kempton (1924). He found that although maize is 
protandrous with occasional exceptions, when a protogynous Spanish 
popcorn was crossed with a protandrous maize, the F1 were pro- 
tandrous. The F2 segregated to give a greater mean protandry and 
a few protogynous plants. Kempton believed protogyny was a 
reversion to an ancestral condition, being normal for Tripsacum and 
Euchlaena ; it was probably not due to photoperiodism, but might 
have been a manifestation of male sterility which may postpone 
anther opening. No male sterility occurred in the flint variety used 
in the present experiment, so that male sterility is not a controlling 
factor here. Other cereals are also known in which the dichogamy 
relations vary according to the genotype: in rye, Vavilov (1951) 
found both cross- and self-pollinated forms, with it normally being 
regarded as outcrossing in Asia. 
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6. METHODS OF REPRESENTATION 


Some consideration may now be given to other methods used by 
maize breeders for illustrating the flowering of maize. Shaw and 
Thom (1951) have plotted silking time-rate (days from 50 per cent. 
silking), positive or negative, against percentage of plants silked. 
In this way they compared a variety’s rate of silking over several 
years. This method could be extended for studying environmental 
and other effects on flowering. Compared with fig. 2, it involves 
additional calculation of differences from the mean, and percentages. 
Yet it quickly demonstrated that for hot, dry years in the mid-West 
tassels rapidly appeared and had earlier pollen shedding, while there 
was a slower silking rate. Tassels appeared more slowly yet with 
more rapid silking in cooler, moister years. Thus, in the present 
data, climate might have influenced the earliest flint plants, which 
suddenly burst into pollen shedding. This is unlikely, however, since 
pollen shedding and silking for varieties sown at monthly intervals 
in England show high correlation regardless of climatic differences 
(Haskell, 19492). ‘That climate does influence the inbreeding-out- 
breeding system of some species is seen in wheat which is 99-9 per 
cent. self-pollinated in England, but shows increased outcrossing in 
drier climates (Vavilov, 1951). 

When studying the effects of sowing date on the rate of sweet corn 
tasselling (1949) I suggested a graph in which probit transformation 
of tasselling time is plotted against the logarithm of the flowering 
date. This helps to straighten the sigmoid-like curves so that uni- 
formity and rate of tasselling can more readily be assessed. Cornfield 
(1949) objected to this on the grounds that probits could be used only 
for discrete values. The effect of the probit transformation, however, 
is to straighten sigmoid curves when ordinates are measured on a 
linear scale of probits instead of percentages (cf. Finney, 1947). 
Hence their use is legitimate here. 

Another graphical way of displaying flowering time data is that 
of Singleton (1948) who graphed the percentage of silks which had 
emerged and of plants shedding pollen against date for four sowings 
of many sweet corn inbreds in Connecticut. Scaling was such that 
100 per cent. silking equalled 50 per cent. shedding pollen ; hence 
critical sowing date for producing hybrid seed is rapidly determined. 
This allows a rapid visual test of whether two inbreds will flower 
together, important when planning cross-pollination for hybrid seed 
production. Singleton’s graphs show that pollen shedding normally 
starts before silks appear. The 100 per cent. tasselling is always 
reached ahead of 100 per cent. silking, and pollen shedding always 
continues for several days after silking is completed. The present 
graph should similarly be useful for comparing the breeding behaviour 
not only of ordinary maize varieties, inbreds and single-cross and 
double-cross populations, but of other plant species. 
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7. SUMMARY 


1. Individual silking date, time of first pollen shedding and 
number of plants of maize together plotted in a three dimensional 
graph give more information on the breeding system than the con- 
ventional plotting of percentage flowering against time. 

2. The breeding behaviour of individuals within a population 
and between two inter-breeding populations is also readily assessed. 
Protandry in maize is not so important for out-breeding as being 
monoecious. 

3. The relation of dichogamy with flowering time is also seen. 
Earliest plants in both populations were more protogynous and 
hence more outbred than later plants, which were more protandrous. 
These shed their pollen while silks are emerging and therefore have 
greater chance of inbreeding. ‘The genetical consequences of this 
breeding system is evaluated in relation to heterosis. 

4. The three-dimensional graph is compared with other graphic 
methods of representing flowering in maize. It is useful for comparing 
the uniformity and breeding systems of ordinary maize, inbreds and 
single-cross and double-cross populations. 
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In these notes I indicate methods applicable to the interpretation 
of multiple linkage tests involving three or more linked loci. The 
methods advocated are based on two extensions of Kosambi’s simple 
and extremely useful theory of linkage relations, the purpose of the 
extensions being to widen the range of usefulness without losing the 
powerful simplicity of that theory. 

In the first place Kosambi’s theory (1) may be supplemented by 
the introduction of a quantity K whose value varies with position 
on the chromosome arm. I have shown elsewhere that we may 
expect high values of K to be realised for sets of loci in the neighbour- 
hood of the centromere, and low values to occur near the termini 
of the chromosome, K taking the value unity at some point related 
to the middle of the arm. In a survey (2) of data on Chromosomes II 
and III of Drosophila melanogaster (not yet fully published), I found 
excellent agreement with these predictions. Thus it seems likely 
that in other organisms the value of K as derived from experiments 
on multiple crossing-over can be of some use in suggesting the location 
on the chromosome arm of the loci concerned. In Part I, I deal 
with the statistical problems involved in the estimation of K and of 
the trend of K in the cases of three and four linked loci, and discuss 
some actual data obtained for the Agouti-pallid region of Chromosome 
V in the mouse. 

Part II is of more abstract interest and extends Kosambi’s theory 
(without introducing K) to multiple recombination with as many 
as six linked loci. Here the mathematical method used is not identical 
with that of Kosambi but proceeds from the general theory of linkage 
relations developed by Sir Ronald Fisher and myself (2-7). 

Part III introduces K into Kosambi’s theory and gives an 
elementary way of generalising Kosambi’s formule to the special 
but useful cases of several fairly closely linked loci such as are treated 
in I. 4. 


|. ESTIMATION OF THE KOSAMBI COEFFICIENT 


1. INTRODUCTION 


The Kosambi coefficient, denoted by K, is defined (5) in terms 
of the recombination frequencies of 3 linked loci, A, B and D by means 
of the equations 





oo Jit J2—Jite ee 
41) 2) 1+2 Ai V2I1+2 
Here the suffixes 1 and 2 refer to the segments AB and BD, while the 
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suffix 1-+-2 refers to the combined segment AD. Thus any three loci 
are characterised by a value of K. It will be seen that K is equal to 
C/2y,,2, where C is the genetical coincidence for the two segments, 
defined by Stevens (8) as 19/712. Kosambi (1) obtained the addition 
formula 


ite 
I+4)iJ2 


by an assumption equivalent to setting C equal to 27,,,; ie. to 
taking the value of K as unity. Regions of the chromosome in which 
K is not significantly different from unity have been called Kosambi 
regions (5), and in such regions interference is said to be at the 
Kosambi level (5). I have given (2, 5) theoretical reasons for supposing 
that in regions near the centromere K has a high value greater than 
unity, and diminishes steadily with increasing distance from the 
centromere, becoming relatively small near the terminus of an arm. 
Consequently there is a Kosambi region, where K = 1, situated 
medianly on the arm. 

These expectations have been well fulfilled in the case of Chromo- 
somes II and III in Drosophila melanogaster (2). Therefore, in the 
absence of evidence to the contrary, it seems reasonable in other 
organisms to use the observed K values as providing some indication 
of the distance of the loci concerned from the centromere, and of the 
orientation of segments relative to the centromere. This note is 
concerned with estimation of K values when data on 3 or 4 loci are 
available. 


Jing = 


2. THE USE OF THREE POINT DATA 


Backcross experiments provide the best determinations of the 
recombination fractions, and therefore of C and K. In such an 
experiment involving two segments AB and BD, suppose that out of a 
total of n organisms bred the observed class frequencies are 


) I 2 12 
2 a a, 442 
Here a, show recombination in neither 1 nor 2, a, show recombina- 


tion in 1 but not in 2, a, show recombination in 2 but not in 1, while 
4,2 show simultaneous recombination in 1 and 2. Thus 


N, = 4,+4 4, show recombination in 1, 

Ny = A+ 4, show recombination in 2, 

1142 = 4,+a, show recombination in the combined segment 
1+2 (AD), 

Nig = 42 show simultaneous recombination in 1 and 2. 


The estimators of ),, 957442) Jig are 


Jy = y/n, Jo = Ng/N, Py49 = My49/N, Fro = Myq/N. 


The sampling variance of any such function is V(») =_y(1—y)/n. 
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The estimator of C is G = NN1>/nNz, and Stevens (Joc. cit.) has shown 
that its sampling variance is 
V(G) = C(1 —Gy, —Cyy— Cry 7. + 209192) RII 


This may be estimated by inserting the estimates ,,...., 6, and a 
crude test of the significance of the difference of C from unity (i.e. 
a test for the presence of genetical interference) may be based on the 
estimated standard error or on the variance. £.g. with the data 


Ce) I 2 12 n 
1116 182 375 24 1697 


we have the percentage values 


$1 = 12°139, Sa = 23°512, CG = 49°551 
Estimated v(6) = 84-2176 (per cent.)? = 0-008,421,76 


A 


Estimated s5(C) = 91770 per cent. 


Thus = 49°55+9:18 per cent. which is significantly less than 
unity and indicative of genetical interference. 

The maximum likelihood estimator of K is simply 

K = C/2f1,_ = n'ny,/2nyngny +9. 
This statistic is biassed, the bias being given approximately by 
p= 8(_8 —3 shen) 
Mig «= n® 

The estimated sampling variance is given by 

V(K) = kt + 9: pM Bd os — MatieS} 

Mig «= My+g = MYMg My MgO n® 

or may be got by inserting the estimates in 
K? I 


: > os pe Ses -4-4K(, 179) 


Dia Sasa a. 
In the example we have 
K = 75°48+15:26 per cent., V(R) = 232°8934 (per cent.)? . 
The bias is not seriously large, since we have 
B = 0:004,041 = 0-40 per cent. 
To test departure of K from unity we calculate 
Xa) = (K—1)?/V(K) = 45317, 
where V(K) is now calculated from equation (1) with K put equal 


to unity. 
K therefore deviates just significantly from unity. 
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In cases such as the present example when K is not very significantly 
different from unity a more refined procedure ought next to be used. 
We assume ),, 2 as free parameters and take 9,9, )42 aS given as 
functions of »,, _¥, by Kosambi’s formule. A maximum likelihood 
scoring technique gives amended estimates which are then used to 
calculate expected class frequencies for a x*,,) fit to the observations. 
If the fit is good, i.e. if x? is not significantly large then we may retain 
the hypothesis K = 1. If the fit is bad we may then estimate K as 
above, together with its estimated sampling variance and standard 
error. 


3. THE USE OF FOUR POINT DATA 


When data on four loci A, B, D, E are available, the appropriate 
treatment is determined by the closeness of the observed linkage 
between them. When the segments are sufficiently long for doubles 
to be at least moderately frequent, the best use of the data requires 
that K be calculated separately for the two triads ABD, BDE, the 
data being reduced to three point data in respect of each triad. If 
the values K,, K, are significantly different then the direction of the 
trend may be taken as presumptive evidence of the orientation of 
the segment AE relative to the centromere O. E£.g. if K,>K, the 
preferred order is OAE. 

To test the significance of the difference (K,—K,) we may 
compare it with its sampling standard deviation using either a y-, 
x°- or t-test. The variance of (K,—K,) is 


V(K,) —2C(K,, K,)+V(K,), 


where C(K,, K,) is the sampling covariance of K, and Kg, given 
approximately by 





ee | "12 No N13 a, a as 
Rik, (2 + 2 4 Bey 
NN = MgNg MyMg— MgMy 49g = MgNgig = My +9Mo+3 


Mies Ming agp dT *) 
NyoNog NyNog = Mglyg =e n 


This expression has been derived by means of a formula cognate to 
that given by Sir Ronald Fisher for the approximate sampling variance 
of a statistic (9). Ifm,, my,..... , m, (totalling n) are the frequencies 
observed in r classes subject to a multinomial distribution, and if 
T,, T, are homogeneous functions of degree 0 or I in m,, mMg,....- , 
m, and n, then we have as the approximate sampling covariance of 
T, and T, the expression 


C(T,, T;) = 2m oT, eT, 9 fk éT, 
tal @m, 6m, én éan- 











ANALYSIS OF MULTIPLE LINKAGE DATA 251 


When T, = T,, this becomes Fisher’s formula 
aT\2 aT \2 
V(T) = Em(Z7)" n=) 
om on 
Each formula holds exactly when T, and T, are linear functions of 
an rae » m,/n. 


4. THE USE OF FOUR POINT DATA FOR SHORT SEGMENTS 
When A and E are close together K, and K, will differ 


significantly. Double recombination classes will occur with small 
frequencies and triples will be absent. For if p, g, r are the three 
recombination fractions the expected frequencies of doubles and 
triples will be given to good approximation by the expressions given 
below (part III), namely :— 


Jia = 2Kpq(p+9), ox = 2Kgr(g+r), dis = 2K pr(p+29+1), 
Dios = 4K%pqr{(p+9)(9+1) —39°}- 
If the observations are 
Oo I 2 3 23 13 12 123 Total 
a 4 j&@G@ 43 43 413 @12 0 n 


then it may be shown by inspection of the maximum likelihood 
equations tor the four unknowns , qg, r, K subject to joint estimation, 
that the estimators approximate closely to the “ intuitive ” ones 


b = (4,+42+443) /n, ete. 
K = (dg3+4,3+4 2) /D = (observed doubles) /(expected doubles). 


Here D is the expected total frequency of doubles calculated as 


D = n( Ji.) 23 +)13) 
= §p+9+r)? —-3(P+9+7°) 
estimated by inserting f, g, 7. Without serious error we may for this 
purpose take for f, g, f# the quantities a,/n, a,/n, a;/n. Thus we get 
the approximate estimator 
K = 3n?(d9g-+413-+442) /2(51° —S3) 
where S, = €;" +65" +4’, (fF = 1, 2, 9, «+ +) 


For the variance of K, we get by the approximate method of the 
last section 


vik) = Re{ Le | 
Ao3+443+44.9 (533 —53) n 


The first term is usually dominant so that to a rough approximation 


V(R) a K2/(observed doubles). 
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Example-—In a four point backcross experiment recently com- 
pleted, of which a full account is not yet published, the four loci 
Agouti, undulated, wellhaarig, and pallid in the house mouse were 
employed. The following class frequencies were obtained. 





























I 2 3 23 13 12 123 n 
From female heterozygotes . | 117 | 185 | 104 I te) te) ° 2503 
From male heterozygotes .| 115 | III 56 I to) I Oo | 2499 

Total . . - | 232 | 296 | 160 2 to) I te) 5002 
( 




















On the combined data for males and females we have 
p = 4°66 per cent., ¢ = 5:98 percent., # = 3:24 per cent. 
Observed doubles = 3. S, = 688, 5,2 = 47°3344 x 104 
S3 = 42°5175 X 108, 533 = 325°6607 x 108, 
Ss = 304°9229 X 107, $5 = 1°5415 X rol 


Thus K = 39°76+23°1105 per cent., V(R) = 0°055,263, and K is 
significantly less than unity, its value being suggestive of a region 
towards the chromosome terminus or at least between mid-arm and 
the terminus. The segment Agouti-pallid is thus between 30-100 cM 
from the centromere according as to what is assumed for the map 
length of the arm. 

The treatment is however rather inaccurate on account of the 
marked sex-difference of two of the recombination fractions, and is 
best carried out separately for the two sexes. We obtain 








K per cent. Information 
Females . ; . 16-24+16°41 37°15 
Males . . ‘ 96-63 +68-29 2°30 
Combined . ‘ 20'92+15°09 39°45 

















Thus in the females K is significantly low, but in the males does not 
differ significantly from o or unity. On the basis of the male and 
female data considered jointly we might legitimately assume the 
segment to be moderately distant from the centromere, since the 
difference between male and female K values is not itself significant. 


5. SUMMARY 


Methods of estimation are given for K and its sampling variance 
in the case of triads of loci, and for the covariance of K values in 
overlapping triads. It is shown how to estimate K in the case of 
three short segments with special reference to the Agouti-pallid region 
of Chromosome V in the house mouse. It is suggested that Agouti 
is moderately far from the centromere. 
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Il, THE FREQUENCIES OF SIMULTANEOUS RECOMBINATION OF 
SEVERAL LINKED LOCI 


1. INTRODUCTION 
Kosambi (1) proposed the formula 
y = #tanh 2x 


as an approximate representation of the relation between the genetical 
map-length x of the segment joining two loci and » the amount of 
genetical recombination between those loci. For 3 loci defining 
adjacent segments x, and x, this formula implies that 


JitJe 
I1+4)i1J2 


where y,, is the frequency of recombination in both simultaneously, 
and j,,, is the frequency of recombination over the combined 
segment x,-+x,. I have indicated (5, 7) the conditions under which 
these approximations are valid. If, in particular, the interference is 
represented by the function 4uexp(—2u) employed in previous 
work (4-7), then, for a segment near the (metrical) mid-point of a 
long chromosome arm, we obtain the relation 


J12 = WiD2i+e = Wi)e2 


y = 4(1—e-* cos 2x). 


For segments less than 35 cM in map-length this function differs 
from 4 tanh 2x to an extent which in practice is quite inappreciable, 
the values of » given by the two formule being discrepant by at most 
8 parts in 302. Thus Kosambi’s formula receives a simple theoretical 
justification as an approximation likely to apply in segments median 
to long chromosome arms. 

Using the methods indicated by previous work it has been shown (7) 
that there are simple expressions giving all recombination frequencies 
for 4 loci (3 segments) which also are likely to apply when all 4 loci 
are reasonably close to metrical mid arm. I gave the new relations 


J13 ~ Js eee, W 12 = WW, 
and J123 = R (www) +$21:%3(1 —e-™), 
ane =}e*sin 2x, w=ytiz, i= V—1, 
and R (wy aWs) = V1) 2)3 1% 2%3 221% 3 Ik 2 


I showed further that these relations approximate to 
I —2yy (1 —29)(1 — 2s} 
1+2)e I+4)iJ3 

and J123 = I1)23 JaJis t+ IsJ12 


These equations may be regarded as natural generalisations of 
Kosambi’s formule. While receiving their justification for short 





J13 =19s| foe 


R2 
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segments only, like Kosambi’s original formule they may often be 
found to apply in good approximation for larger segments. It will 
be seen that their form is consistent with the known disappearance 
of interference between long segments or widely separated segments. 

It may be remarked in passing that the formule for 4 loci are 
adequate for the derivation of all formule required for 5 loci (4 
segments) in virtue of identities of the type 


J1293+4 = J123 t)1204 —2) 1234 


In the present note I extend the results to give formule in y, z 
and x adequate to deal with as many as 6 loci (5 segments), and show 
how the formule may be adapted to apply to cases in which any or 
all of the pairs of segments are non-adjacent. 


2. FORMULA 
(a) 4 loct (3 adjacent segments).— 


Jig =)1J2—-%1%2, = D23 = J2I3—% Xa) 
Dis =I1I3 2% 3e™, 
J123 = I1)2.)'3 J X2%3 221% 3 32122 +$21%3(1 —e-™), 


(b) 5 loct (4 adjacent segments) —Only one formula is essentially 
new, namely that for 934, all the others resulting either directly or 
by mere change of variables from those for 4 loci. However, it is 
instructive to write out the full set. 


J12 F=I)1) 2-2 1% 2» Ira = I1 Ia 21g * 
J23 = J2)3—< Ka» Joa = JoJa—Zarse"™ 
Jia =i J3—212%3¢-™, J3a = J3)4—Z Ka» 


J123 = Ji J2I3 JDrxs3 Dr vrsJararethirZa(1 =), 
J124 =r J 2a age 4% ge — Iga +42 1%4(1 —e~M*)e-™, 
Iisa = Ii Isa — IX 94 Dah Zge* Ig ge" +4224 (1 —e-™)e-™, 
J23a = J2)3 JJ rh 3% 2% 4 aX 0%s th20%q(1 —e-™*), 
J 1234 = D1 JD 2I3)4 Di X34 D2 38 1%4 D1 3% 2% 
Di Ja 2%3 V2 Jak 123 D2 Jah 2 F212 2% 3% 
+4 1Z0%q(t —E-*) +3922 3%4(1 —e-™) 
+4 Y921%4(t —€-) +4 942123 (1 —e-™) 
—41%4(1 —e-™) (1 —e-™), 
(c) 6 loct (5 segments)—Only two essentially new formule are 
required, namely : 
J135 = Di) 3 Ise atse * 9h rZge 4 + ™ 
52,2 9¢ 8 +4242 5(1 nA) eet) 
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J12345 = D1 D2) 3 IIs — 22122) 3 a's 
—3{Z1%5 a Ja(t —e-*) +2325 Vout —e-™) 

+2325 Yo Ja(t —e-™*) +2323 4 a(1 —e-™) 
+230 Va Yo(t —e-*) +2324 V9 J5(1 —e-™) 
+2324 )2.)5(1 —e~*) +-2925.919(1 —e~™) 
+2085 )1Ia(t —e-**) +2 525.1 Yo(t —e-™)} 

—${2 1232425 (1 —€7™**) +2292 425(1 —€-™*) 
+2 1222325 (1 —e-™**) } 

+3 (21252 +Z0%s51) (1 —e-) (1 —e-™) 

+4(2325 V3 +23%5)1) (1 —e-™*) (1 —e-™) 

+4(21%5 Ya +2425.01) (1 —e-*) (1 —e-*) 

—}z,25(1 —e-) (1 —e-**) (1 —e-*), 

(d) Extension of formule to separated segments.—It is easy to generalise 
any formula relating to contiguous segments so that it applies to 
non-adjacent ones. For example, the expression for 435, relating 
to intervals x, *3, x5, x, out of the contiguous set x4, X2, %3, X4) %5) Xg» *7 
can be derived as follows : 


If xo, x4, %, were all zero so that 1, 3, 5 and 7 were contiguous we 
would have 


21357 = DiI 3I5)7 — 218357 218537 — 218135 
—£3h5 D1) 7 238715 Sh 1 S83 
+2 1% 92589 +8127 Io(1 —e-) +3212) 9(1 —e-*) 
+4225 92(1 —e—*) +32 527.1(1 —e-™) 
—3212,(1 —e-™) (1 —e-™). 


When x9, x4, %, are not zero this formula becomes 


J1357 = I1I3IsI2 2% a Io Ire —Z 25 I Jqe ht 
— 227g Vso 8 Tt +29) — 202, yy oe 
Zhi W504 — 2520 Dye 
+232 32527674 +2) 
+4232 I5(1 —E- 82) e—4% +41 +40) 4 he 2, ya(1 —e- St) eA He te) 
+432325.9(1 —e-H)e-Mm t0) 1 dee, yA (1 —e~ 4%) g— Ata +44) 
—}e4z_(1 —e- 8) (1 —e- 4) ¢—- Ar tas +) 
The rule of derivation consists merely of inserting an attenuation 


factor corresponding to each term z, z, as written and detemined 
by the total length of “ gap intervals ” between the segments j and k. 
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3. METHOD OF PROOF OF THE FORMULA 


For (k+1) loci (& contiguous segments x,, %2, ..., *,) infinitely 
distant both from the centromere and the terminus of the chromosome 
arm we have a generating function (5) which takes the value 


P(A,, Ag. «+» Ay) = $Q exp( = D;), 
)s 


where D; = 2x; (j = 1, 2,. 


~~ 


* 


and Q = 2C€,C,...C,+258,C,...C0,L,+28,8,C,...C,L3.+... 
a +++ S82... SpLie-- 4 
with C; = cosh 2VXj, S; = sinh aVA,. 
In a typical item such as S,S,C,S,C;C,S,C,L,.,, the L symbol 
L047 is to be read as the function 
ce Vr, , VAs v™) 
Vag Va. VA VA 
Similarly Loysgg in CyC,S,8,8,S,C,S,C,CoLsaseg is 
(vs Jv... tthe ) 
VaVvi | VuVNVS 
The multiple segmental recombination fractions can be calculated 


from the 2" quantities obtained from P by substituting values +1 
or —1 for the A’s (5). For example, if we set 


Ae, = An, = eee Ag = OI 


1 rT 








and all other A’s equal to +1, the resulting value of P is equal to 


en ae I —2 In, +m+ ooo Mr 


where Vn, +n,+... +n, 18 the probability that the sum of the numbers 
of exchange points in the segments x,,, X,,, -- 5 ¥n, Will be odd. 
It is easy to verify the elementary identity 


F Se Sah, ey —22 Yn +42Inn, — o. +( —)"2"Fan,... Ny? 


where Jn,n, is the probability that segments x, and x,, will show 
simultaneous recombination, and so on. Thus the segmental recom- 
bination fractions Jan Inns +++» Jig+++~% can be calculated 
inductively. 


For example, with 3 adjacent segments we have 
1—2( 91 +93) +4)13 = 1-27.43 = Pas = P(—1, 1, —1) 
= e~'P:+P)(cos D, cos D, —e~?”* sin D, sin Dg) 
= (1 —2y3)(1 — 273) —421z9e-™. 
Thus J13 = I1I3— 2130 ™. 
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Similarly 


1 —2( yy tots) +4( 1223113) —8)123 
= 1—2),40+3 = Pigg = P(—1, —1, —1) 
= e~(PitD:+D)(cos D, cos D, cos D,—sin D, sin D, cos D, 
—e~*”: sin D, cos D, sin D, —cos D,, sin D, sin Ds) 


which ultimately gives the value of 7,93. 


4. FORMULA IN TERMS OF RECOMBINATION FRACTIONS ONLY 
The substitutions 
e~ “+ x= (1 —2y,) /(1+29,) 
ZX, = (1 —2y,)(1 —29,) /(1+4.9,Is) 


give formule generating those of Kosambi (1) and Owen (7). With 
(k+1) linked loci they could be used for fitting the k parameters 
D1» 99 + + +5 Jn to Observed recombination data by maximum likelihood 
scoring. 


5. THE AUXILIARY QUANTITY z 
Z = te-™ sin 2x 
is the imaginary part of w where 
1—ow = e~™1+i), 


so that w = y+2z satisfies Trow’s formula for adjacent segments at 
mid-arm. 
1.€. Wi+9 = W1+W,—2WyWs 
Thus the effect of interference is formally removed for adjacent 
segments by the use of w, which may be called the complex recombination 
fraction or “‘ complexit”, its components y and z being named respectively 
the “ recombit” and the “ cocombit”’. 
For short segments we have the expansions 

y = $(1 —e-*cos ax) = x—$x84... 
as in Kosambi’s formula. 

z= fe-™ sin ax = x—2x2+$84... 

x= y+$8+... 
Therefore z= y—2y?+$ P+... 

= y(t —2y) /(1 —$.9*) +0( 94). 

The second order approximation to the cocombit is therefore _y(1 —2y), 


but its use is rather less effective than that of the approximation for 
Z,Z,- It may be noted, however, that a good second order approxima- 
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tion is provided by (1 —1°56y). The following table gives percentage 
values of x (cM), of» and of $ tanh 2x and also of z and of _y(1 —1-56y). 




















M 4 tanh 2x 5) z (1-156 9) 

o o 0 oO o 

5 4984 4984 4517 | 4°597 
10 9°86 9°879 8-133 8-357 
15 14°566 14°614 10°946 11°282 
20 18-998 19'I30 13°052 13°421 
25 23°106 23°386 14°540 14°854 
30 26°852 27°352 15°494 15681 
35 30°218 31010 15995 16-009 
40 33°202 34°348 16-116 15°943 
45 35°815 37°304 15'924 15°585 
50 38-080 40°062 15°478 15'025 
55 40°025 42°451 14°833 14°339 
60 41°683 44°543 14°036 13°59! 
65 43°086 49°355 13°130 12°834 
70 44°268 47°904 12°150 12105 
75 45°257 49°211 11*129 11°432 
80 46-084 50°029 10°09! 10°983 
85 46°770 51175 9:058 10°320 
go 47°340 51°894 8-049 | 9°883 
95 47°812 52462 7077 | 9°527 
100 48-201 | 52°904 | 6°153 9242 











6. LOC] PROXIMAL TO CENTROMERE ON A LONG ARM 


It is easy to derive results applying to segments on a long arm 
which are at finite distances from the centromere. We replace the 
map distance x, wherever it appears in the existing formule, by the 
interference metric ¢ and treat the centromere as having the properties 
of a localised exchange point (4, 5). 

Consider 3 adjacent segments 0, 1, 2 which have the metrical 
specification (0, ty), (to, t1), (t1, t2). We have 


Joa = Iona —Zobae-*™ 
where 
n = 4(1—e-* cos 2D), { = $e?” sin 2D 

and D is now metrical length, i.e. D, = (t, —ty). 

The probability of recombination in 2 conditional on recombina- 
tion in 0 is 

Jo2l Jo = N2—€- 4” £2(Zo/ Io). 

The limit of this probability as t)—>o is clearly the probability of 
recombination in 2 conditional on there being an exchange point 


at ¢= 0, and therefore gives the recombination over a segment 
(t;, t2) at metrical distance ¢, from the centromere as 


Jo = M—e fy 
as given by Owen (4). 
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Similarly for simultaneous recombination in separated intervals 
we have 


Jaa = MaMa Laka 0 Cag + Gaye Het ere}, 
Also 


Jre3 = NNs —(N1b0b3t+0b1bs+3bib0) +h0i03(1 —e- 4) 
—e~*{nyobst+nanssi tninase —Sibels 
—4304(1 —e~**) —nola(1 —e- 4) 
—Hoana(t —e-") —30,(1 —e-*) (1 —e-)}. 


7. SUMMARY 


Generalisations of Kosambi’s formule are obtained serving to 
give the frequencies of multiple recombination in as many as 5 adjacent 
or separated segments. The formule are expressed in terms of 
recombination fractions for single segments and auxiliary quantities z. 
It is shown how to obtain complete generalisations of Kosambi’s 
formule, expressed purely in terms of the recombination fractions. 


Ill, ELEMENTARY FORMULA. FOR MULTIPLE RECOMBINATION 
OF CLOSELY LINKED LOCI 


1. INTRODUCTION 


It has been seen that for contiguous segments 1 and 2 Kosambi (1) 
gives the useful addition formula 


Disa = (Mit) 2)/(1 +41) 2) 


connecting the recombination fractions y,, 7. and _y,,. in the segments 
1 and 2 and the joint segment 1+2. This has been found to be sub- 
stantially correct for long segments in the mouse (10), in rice (11), 
and in the sweet pea (12) and to hold for loci situated medianly in 
the arms of chromosomes VI and III in Drosophila melanogaster (2). 
Theoretical reasons have been advanced why this should be so (5, 7). 

Since the frequency of simultaneous recombination is given always 
by 

Jia = ¥(I1tI2—Ii+2) 
it follows from Kosambi’s formula that the frequency of doubles is 
given by 
J12 = 22 Ir tI2) (1 +4 I) 


It was shown above that for non-contiguous segments an appropriate 
generalisation of this result is 


I—2g (1—2)3)(1 2s) | 
I+2 2 I+4)1J3 





Pat eee IIa} 1— 
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For short segments when ),, 9, 3 are small, if we retain only terms 
of the third order we get the formule 


Jia = 22a i te), Dis = 2I1)s(J1+2I2+3s)- 


Similar formule could be got for the frequencies of triple, quadruple 
crossing over, etc. using the previous results. However, Professor 
Sir Ronald Fisher suggested to me that they could be derived rather 
simply by a completely elementary method, which may be illustrated 
by the case of 4 linked loci (3 contiguous segments). 

We note that if segment 3 is so short that y, may be replaced by 
the differential dy then the probability of simultaneous recombination 
in segment 1 and the segment dy is 


Dis = 2I1(Jit2I2)h . . 2 Gy 


Equation (1) also implies that the probability of recombination in 
1 conditional on recombination in dy is 


I13/4y = 293(91 +22) : : - (2) 


Results (1) and (2) enable the frequencies of simultaneous recombina- 
tion in any number of short contiguous segments to be calculated 
inductively provided that a further simple assumption is made. This 
assumption is in accordance with recombination theory as previously 
developed. We suppose that a cross-over point once established 
effectively divides the chromosome into non-interfering regions. 

Let now f, g, 7 be the frequencies of recombination in contiguous 
segments AB, BC, CD supposed short but not infinitesimal. Consider 
an elementary segment dy located at a point P interior to BC. Then 
to the present approximation the recombination over BP is » and 
over PC is (¢—y). By (1) the probability of simultaneous recombina- 
tion in AB and dy is 

2p(p+2y)dy. 


But by (2) the probability of recombination in CD conditional on 
recombination in dy is 


ar(r+2q—9). 
Thus the probability of recombination simultaneously in AB, dy and 
CD is 
4pr(p+2y) (r-+2q —2)dy. 
It is now easy to see that, on account of the assumed smallness 


of p, g and 7, the probability of simultaneous recombination ),, in 
AB, BC and CD is merely 


q 
[‘ser(p-+a9)(r-+20—ay)h 
which is 


4pr{(p+9)(q-+r) —4q*}. 
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The results of this method are assembled for 4, 5 and 6 loci in 
section 3. It will be seen that a parameter K has been introduced. 
This will be discussed below. 


2. THE KOSAMBI COEFFICIENT 
The Kosambi coefficient K for three linked loci is defined as 


K Be ge ee, ie ; 
4A Vi D2) 1+2 Wie 1+2 

The value of K is in general different from unity. K is equal to 
unity, if and only if, Kosambi’s formula holds exactly. The value of 
K as estimated from 3 point data may be used to test agreement with 
Kosambi’s formula. When K is significantly less than unity then 
the genetical interference which is operating must be concluded to 
be more intense than that predicted by Kosambi’s formula. The 
opposite conclusion is indicated by an observed K significantly 
exceeding unity. Reasons have been given (6) for supposing that in 
regions near the centromere K is large, and thereafter diminishes 
steadily becoming less than unity in regions near the terminus of a 
chromosome arm. At some intermediate region (actually near 
metrical mid-arm) K is approximately unity and Kosambi’s formula 
will apply there. These expectations are confirmed by the data on 
Chromosomes II and III of Drosophila melanogaster (2). Thus an 
observed value of K is some indication of the location of the triplet 
of loci on the chromosome arm. 

Inclusion of the parameter K modifies the formule for short 
segments in a simple way and confers generality upon them. 

If K is defined as above then we have identically for the frequency 
of doubles 


Dra = 2K aise = 2K a(I1 2) /(1+4Ky V2), 
which for short segments is effectively 


Ire = 2Kya(9ib2) = 2K (921 2") 
Now suppose that the segment y, is increased to be y,+dy by the 
addition of an elementary segment dy. Then the resulting increase 
in )j9 is given by 





oK 
D2 = 2K (Ii +292) +291 Ia( Jat I2) 5 dy 
2 


The second term is of the third degree in ),, y, and consequently 
may be neglected in the present approximation. Thus 
dy. = 2Ky,(91 +22) 


and the probability of simultaneous recombination in AB and dy 
becomes 


2K p(p+2y)dy, 


K appearing merely as a multiplier. 
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To the present order of accuracy therefore, K is to be regarded 
as a constant over the total region occupied by the segments and the 
formule are merely multiplied by powers of K. 


3. FORMULA 


(a) Notation.—p, q, 7, s, t, denote contiguous segments taken in 
order. It is convenient to define the following symbols. 


(p, 9, 7) = (p+9)(9+r) —4¢? 
(P 9%, 8) = (b+9)(9-+1)(r+5) —39°(r+5) —$7°(p+9) 
(P97, 5, t) = (p,q 7)(1, 5, t) —$r?(p+q)(r+#) 
(b) 4 loci (3 segments).— 
Jie = 2K pq(p+9), Jes = 2Kqr(q+r) 
Jis = 2Kpr(p+2q+7) 
Dias = 4K*pqr(p, 9, 7) 
(c) 5 loct (4 segments).— 


Jig = 2K pq(p+9), Jis = 2Kpr(p+2q+r) 
Jia = 2K ps(p+2g+r+s),  Je3 = 2Kgr(q+r), 
Jog = 2Kqs(q+5), Isa = 2Krs(r+5), 
Joss = 4K qrs(q, 1, 5), Jios = 4K*pqr(p, 9,7); 


Jive = 4K%pgs{(p+9) (g+2r+5) —3¢"}, 
Disa = 4K%prs{(r+5) (p +29 +1) —477}, 
Jri2osa = BK%pqrs(p, 9, 7, 5). 

It will be seen that with the exception of the last three, these 
expressions result mutatis mutandi from the appropriate 3-point formule. 
For instance ,, is derived from 2Kpr(p-+-2q-+r) by the substitutions 
s for r, g+-r for q (q+r being to the present degree of approximation 
the recombination in the separating segment (2+ 3)). The formule 
for 7394 and _),34 contain nothing essentially new being calculated from 
the relations 

J124 = J123+4 D123, D134 = J1+234 Vasa 
which hold in the present approximation. Their form shows them 
however to be derivable by a simple rule. We may write 
Dies = 4K%pqs(p, 9, 5)" 
where (p, 4, s)! results from (p, g, s) by inserting in the item (qg+5) 
the term 27 corresponding to the separating segment r. 


(d) 6 loct (5 segments)—Formule for two segments such as 


Jigs = 2K pt(p+2g+e2r+e2s+t) 


follow trivially from the previous results. 
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The formule for three or four segments follow either directly or 
by application of the rule for separating segments, e.g. 


Y1as = 4K %pqt(p, q, t)* = 4K*pgt{(p+q) (q+2r+25-+t) —497}, 
Ji > 4K*prt(p, ’; t) 1 4K*prt{(p+2q+r) (r+25+1) —jr*}, 
i235 = BK%pqrt(p, g, 1, t)? 

= 8K%pqrt{(p+9) (q+1) (r+25-+t) —39°(r+25-+t) —}7(p+9)}, 
J12%5 = 8K *pqst( p, q, §; t)? 

= 8K %pqst{(p+9) (q+2r+s) (s-+t) —}q?(s-+t) —}s?(p+q)}. 


The formula for 5 segments is essentially new and can be found only 
by integration. 


Jiesas = 16K *pgrst(p, 4, 7, 5, t) 


4. VERIFICATION 
The expression 
D123 = I1D2I'3 JF 2%3 2k 1%3 I 3% 1% a +42123(1 —e-™), 
given previously (loc. cit.) reduces to the form 4pqr(p, g, r) on putting 
z= y—2y?+$9%, 
e~4* — 1 —4y+8y? — 1678+ ... = (1 —29)(1 +27) 
= exp{— 4(9+39%)}+0()4), 


and retaining only terms of the fifth degree. This serves as a mutual 
check of both formule. 


A similar comparison may be made with the generalised Kosambi 
formula given by the writer (loc cit.) namely 
Ji2es = P23 —-Wis TV 12 
with )49,.%13 given by Kosambi’s formula and 
I—2qg (1—2p)(1 i) 
he = Pit ——2 — ——- 
713 p 1+29 1+4pr 
With small segments we have correct to the 5th degree 
Jias = 4pqr.2q(p+q+r). 


Consequently the generalised Kosambi formula implies rather a 
higher frequency of triples (about 50 per cent. more). This is because 
it is essentially a quadratic approximation and not a cubic one. 


5. SUMMARY 


An elementary method is developed by which the frequencies of 
simultaneous recombination of several closely linked loci can be 
specified, and formule are given for as many as 6 such loci. 
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1. INTRODUCTION 


IN studies on genetic linkage, the primary object is to form estimates 
of the frequencies of the modes of gamete formation from data on the 
gametic output of individuals of known genotype. With disomic 
inheritance one parameter, the recombination fraction, enters into 
the law of inheritance for pairs of linked loci and in the absence of 
differential viability, this can be estimated from the gametic output 
of either of the two types of double heterozygote. With polysomic 
inheritance, a multiplicity of segregating genotypes and of modes of 
gamete formation makes the analysis more complex. 


2. TETRASOMIC INHERITANCE AT TWO LINKED LOCI 


With two linked loci and tetrasomic inheritance, there are eleven 
modes of gamete formation. An example of this genetic situation 
has been discovered in Lythrum salicaria (Fisher, 1949). Fisher (1947) 
has considered this case in some detail. He has shown how, in general, 
the heterogenic genotypes can be classified into isomorphic sets such 
that all members of a set can be generated from any constituent 
member by systematic gene substitutions. It follows that genotypes 
which belong to the same set yield populations of gametes with 
equivalent frequencies. In accordance with a generalised form of 
Mendel’s law, there is determined for any genotype, a gametic series 
for each mode of gamete formation. The gametic output is then 
expressible in the form of a matrix with the rows corresponding with 
the different modes of gamete formation and the columns with the 
different gametic genotypes. In most cases, the rank of the gametic 
matrix is less than the number of modes of gamete formation so that 
the gametic output from a particular isomorphic set of genotypes is 
deficient in some of the information sought. A detailed analysis is 
required in order to find out what information becomes available on 
considering genotypes from different isomorphic sets. The results 
of the analysis for two linked tetrasomic loci are summarised in fig. 1, 
which is reproduced from a wall diagram in the Department of 
Genetics, Cambridge, by the kind permission of Sir Ronald Fisher. 

With only two alleles at each of the two loci, the frequencies of 
all eleven modes of gamete formation cannot be estimated. However, 
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six can be estimated as well as four independent linear combinations 
involving the other five. 


Fisher’s graphical analysis of the 24 gametic matrices in tetrasomic 
inheritance with loci 
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Fic. 1.—Numbers of isomorphic sets having gametic matrices of different rank, with the 
particular components of deficiency in which they differ. 


3. HEXASOMIC INHERITANCE AT TWO LINKED LOCI 


It is of some interest to consider the problem presented by linkage 
between two loci showing hexasomic inheritance. Although hexasomic 
inheritance has not yet been demonstrated in any organism, it is not 
unexpected in autohexapioids. 

It might be noted here that polyploids are conveniently classified 
both according to the origin of their germ plasm and according 
to their breeding behaviour. If polyploidy has originated from a 
cross involving only one species, the new individual is an auto-polyploid 
whereas if more than one species was involved in the cross, it is an 
allo-polyploid. Polyploids which on breeding always behave as 
diploids are properly called amphi-diploids. These are to be 
distinguished from polyploids which exhibit polysomic inheritance 
and in which, at meiosis, a gene has more than one possible partner. 
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Genetical behaviour might not be of the same nature throughout 
the whole germ plasm and it is conceivable that inheritance might 
be polysomic in one set of chromosomes and disomic in another. 

With two linked loci and hexasomic inheritance, there are forty 
modes of gamete formation. Using Fisher’s enumeration (1947, 
P- 72) these are shown in table 1 where, for each mode of formation, 








TABLE 1 
Mode of Typical Mode of Typical 
formation gamete formation gamete 

I b,/agbq/a,b3 21 1b /agbq/agb, 
2 4b, /ayb;/agbq 22 @1b/agb9/a3bq 
3 yb /agbg/agb, 23 1b /agbz/agb4 
4 1b /agb,/agbq 24 yb, /a3b9/aqbq 
2 @b,/ab,/a,b, 2 @b,/abg/agbq 
1b /ab5/a4b, a b,/a,b,/a4b, 

7 @b,/a,b3/a9b, 27 yb,/a,b9/a,bq 
8 1b /agb3/a3bq 28 a;b,/a9b,/a9b5 
9 1b, /agb3/asb, 29 ,b,/a,b,/agbs 
10 b,/agb3/agb; 30 @b;/a9b;/a3b, 
I @b,/a;b9/agby 3! @ybg/a,3/agb5 
12 bq/ayb4/agbg 32 @gb,/a5b,/a5b4 
13 ybg/aybo/asbq 33 4199/0 103/0304 
14 @1b3/agb,/agb, 34 gb,/43b;/agb, 
15 @be/ayby/agb4 35 @yb/a,b3/a4b, 
16 b9/a3b,/a5bq 36 @gb,/a,b,/a,b, 
17 b9/agb,/a9b, 3 b9/a,b3/a9b, 
18 @yb4/agb3/a4b5 3 @gb,/a3b,/a,b, 
19 be/agb5/ayb4 39 @yb_/a,b2/agbs 
20 4ybg/ab5/a3b, 40 yb, /a_b;/a3b, 




















there is given the genotype of a typical gamete that would arise from 
the parental genotype a,b, /a,b./a3b3/a4b,/a5b;/agb,. 
There are sixty hexasomic genotypes with two alleles at each of 
TABLE 2 


Numbers of isomorphic sets of digenic hexasomic genotypes classified according 
to the partitions at two loci 








Partitions (51) (42) (3?) 
(51) 2 2 : 
(42) 2 3 7 
(3?) I 2 2 




















two linked loci and these fall into seventeen isomorphic sets of which 
seven are self-conjugate, the remainder comprising five conjugate 
pairs. These sets can be classified as shown in table 2 according to 
the partitions at the two loci. 

The number of gametic genotypes to which any parental genotype 
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can give rise is determined by the partitions at both loci. These 
numbers are shown in table 3. 











TABLE 3 
Numbers of gametic genotypes to which any parental genotype can lead 
Partitions (51) (42) (3°) 
(51) 13 16 16 
(42) 16 20 20 
(3?) 16 20 20 











4. THE GAMETIC MATRICES 


Since parental genotypes belonging to the same isomorphic set 
have equivalent gametic matrices and as each of the gametic matrices 
for two conjugate isomorphic sets can readily be developed from the 
other, it is sufficient if we construct gametic matrices for one parental 
genotype from each of the seven self-conjugate sets and from one of 
each pair of conjugate sets (see tables 5-16). 

As the rank of any one of these gametic matrices cannot exceed 
two less than the number of gametic genotypes which enter into it, 
no parental genotype used alone can lead to estimates of more than 
eighteen independent linear combinations of the forty frequencies. 
The actual values of the ranks of the twelve gametic matrices are 
shown, with typical parental genotypes, in table 4. 


TABLE 4 
Ranks of gametic matrices 








(51) (42) (3*) 
(51) | AB/ablabjabjablab 11 | Ab/Ab/aB/ablablab: 13 | AB/Ab/Abjabjablab: 11 
Ab/aB/ab/ab/ab/ab 11 | AB/Ab/abjab/ab/ab: 13 
(42) me AB/AB/ab/abjablab: 16 | AB/AB/Ab/ab/ab/ab: 14 


AB|Ab/aB/ab/ab/ab: 16 | AB/Ab/Ab/aB/abjab: 6 
Ab/Ab/aB/aB]ab/ab: 16 


(3°) AK eS AB|AB/AB|ab|ab/ab : 10 
AB/AB/Ab/aB]ab/ab : 10 




















Beneath each gametic matrix there is indicated a series of gametes 
whose classification would lead to the maximum number of 
independent linear combinations available for estimation from any 
parental genotype which belongs to the given isomorphic set. The 
number of gametes in any such series is, of course, equal to the rank 
of the gametic matrix. To take an example, the self-conjugate set of 
four genotypes which includes bisimplex coupling (AB)(ab),;, can 
yield estimates of eleven independent linear combinations and at 
the bottom of table 5 there is indicated a series of eleven columns. 
In all, 147 columns are marked out in this way in the tables 5-16. 
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TABLE 5 


Parental genotype AB/ab/ab/ab/ab/ab, representing a set of 4. 
Row divisor: 6 


Composition of gamete 


























AB aB ab Ab | AB AB ab ab | AB AB ab ab | AB AB AbaB | AB ab AB AB 
AB aB ab Ab | AB ABab ab | Ab aB aB Ab | AB ab Ab aB | Ab aBaB Ab 
AB aB ab Ab | Ab aB aB Ab | Ab aB aB Ab | ab ab aB Ab | aB Abab ab 
I 3 : , ate - 
2 a's ; ; y 4 ; 
3 2 ae ‘ oy 'g ; 
4 4 ; ; ‘ +? Erk 
5 3 ’4 ‘ ae viua ; 
6 I 22 : I oe 4 
7 3 : 4 I rox 
8 3 ; ; I - 2 
9 2 a I ud 
10 3 ia ‘ a 
11 4 os ae 
12 : 3 3 r 
13 2 > xo , 
14 2 y 4 3 2 
15 2 nH ry ‘ 
16 3 op a I 
£17 4 , ; 8 
§ 18 I 22 I 
2 19 2 r 3 2 
3 20 3 : 3 
S ot: 3 I ; I , I 
© 22 3 3 : I a re 
kt 23 | 2 e of I F , 
2 24 2 oa 28 ae? I : . 
25 2 21 at , , I 
26 2 I 2 : a BM 
a7 3 I ee a ‘ I 
28 3 ee I a 
29} 3 ke , I I 
30 3 ~ 4 ; ae 
31\ I e I ‘ 
32J I 9 ee ‘ 
33\ 2 e. ‘i I 
34 . ‘oe  % I 
35 2 21  % I 
36 2 I , Ae 
37 3 I : He. I 
38 3 I ; . a I 
39 gy I e I I 
40 3 I Ref I 
* * * * * * * * *£ *& & 


$2 





270 J. H. BENNETT 


TABLE 6 


Parental genotype Ab/aB/ab/ab/ab/ab, representing a set of 4. 
Row divisor 30 


Composition of gamete 











AB aB ab Ab | AB ABab ab | AB ABab ab | AB AB AbaB | AB ab AB AB 
AB aB ab Ab | AB AB ab ab | Ab aB aB Ab | ABab AbaB | Ab aBaB Ab 
AB aB ab Ab | Ab aB aB Ab | Ab aB aB Ab | ab ab aB Ab |aB Abab ab 
I 6 9 9 : E 6 : 
2 12 4 4 4 4 2.4 ‘ ; 
3 | « 6&6 ° ie the 6 ; 
4 12 ee 4 4 ie ty ; 
5 13 + 8s 4 4 i 2-4 p ; 
6 <i — 9 F , os ~« & ; 
+ ‘ 13 4 4 g 3 uma : ; 
8 : ia : e 4 ; 
9 7 . Bes a 2 5 ‘ 
10 13 oe 4 4 : en Ree 
II . 32 a ee 4 4 ; Seca ae 
12 . 9 ~ © 'D pa 3 . was 
13 - 14 3 3 g 9 Frey aye ‘ 
14 cee ‘ 8 ‘_* - 8 Sa ‘ 
15 « 34 3 3 3 3 mie , ee 2 
16 . 29 4 4 3 3 ee ae as CE 
£17 12 on te 4 4 a , 
3 18 8 ; a oe ; ¢ 5 
é 19 7 - 8 8 3 4 
6 - 9 9 3 3 
r) et} 9 II 3 . 3 I . 
@ 22 9 3 11 S ie I 8 GaAs 
3 23\ 10 10 3 . @ I a = oo 
2 24J 10 3 10 2 I 2 25 
25\ 10 10 2 + ce Care 
26 10 2 10 ae I I o his 
27 13 — 3 4 ; : EE Wee 
28 13 . 2 ae 4 3 : ‘ : 2 wae 
29} 9 II 3 9 oe oe a ere 
go 9 © @ oS 2 ae | PR ay 
31} II ee . 8 t 4 @ tg 
32 II . 235 Qa > 3 2 Gos 
33 10 .%0 8 we i ee rs) 39 
34 10 + Sap a % een: 8 Ge 
35\ 10 - 40 2 3 eS e 4) «@ 
36) 10 - 210 § - Re Sat o 2s 
37 9 ~ Fy 3 Stas FGI ve 
38 9 3 11 ek Be R Gls 
39 13 3 4 4 3 ae : 7 
40 13 4 3 3 4 ». eae 












































—_ 
™ x) 
os 
a I3-3 a+ a a e+ 6 6 s ~ ~~ COM OO OE CO Om Oe me ts % 
393 awk a i ae a ~ ~ mm COON OE mt Om Oe Om mt tt * 
eau | te see 8 OOM MO + HOD HIND Moms MEM OH OH CH OH CY OF ene * 
* w2 gh Se tae am OO ~~ “ ~~ —~ ae * 
" > - 
5 3 933 eS Ue a Se Gal ap ee ae ee ~ oO erm 6s: ae oe WS . 6 we 4 ee om ae COCO 
s s S39 CS 8 Oe Oe ee. eS Oo Quam ss 66 Vee eb es wo Oe 
Ps, = J33 Dts SOON Om ~ OOoa ams eT a a ee ee eee 
ica 8 I33 Oh Om a. ee a ee ~~ a a 2 oo ee ee ee ee ee 
o S = 
Zz = v 
— = S 2S 
go. # Sux hs ee) _— OO: “on ao - an = + = COO me Lal co 
~ 
g «~ < .. & 339 ia ese 8 ” oF 0 ee a, ee) Om + O89 Com mt root * 
i? fs <§ 2° 338 ae ee Mees ee 6s we we ee ese Sime ea cine os Py 
nN aus a) SS ee ee Mie we Ce es on a ars ce ee a <a oe on * 
r sa. f(64 
eS S28 8 
so} a> § Sez a + one ao ane a CO HO HOD OM HAH OM 
z, a S) 339 a + oe oe ano an HHO Oe HAM mM * 
ws 3393 ~ eS ee A ee ee ee ee _ - * 
| & WN ~ ° Ore 0h) 6 Oe 8 ee M6 e 6S SM 6.6 8. « 4-6 = _ * 
= 
sd S 
oe ee Se ee ee a ae ee oe ee ee oe ee Oe ee ee ee a a cr ee er 
Z 3 I33 sa ~ e 
= 5 ce CHO WO OD + OOD ONO em mt OO OD ON OND ee OO OND . 
_ 
Ss 
nw 999 Sys ea ete S28 ed OY ea Lee a ee BO Ce ae oe oe ee ae e 
QQe ree Ree a ee oe a eee Ree ES ee a eee Er ee re a * 
DQ PPL ODO ODO em mm, cam 
= A OOPINS BO DOH A OMPHINO HO HDS = A OOMINO HO DOH A OMI HO DOS 
one ee ee. ee 
UOTFVULIOJ JO Ipoyy 
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Composition of gamete 


TABLE 8 


Row divisor : go 








Parental genotype Ab/Ab/aB/aB/ab/ab, representing a set of 4. 


























AB AB ab ab | AB AB ab ab | AB AB AbaB | AB ab AB AB 
AB AB ab ab | Ab aB aB Ab | ABab AbaB aB aB Ab 
Ab aB aB Ab | Ab aB aB Ab | ab ab aB Ab Ab ab ab 
a. aT cia 9 9 a ie - 9 9 ae ak 
6 6 6 - 12 12 - 12 12 - 12 12 «tah aie 
> ee o a a8 Piet a, 297 7 30 4 4 
6 6 6 . 12 12 - 12 12 at hime mes 
gi2 3 oe 15 15 699 tm 
. os - a5 88 2 2 6 5 5 26 6 6 
312 3 15 15 6 6 . 2 ge. oe yet oe 
ay are 9 9 9 9 -m 9 8 . 24 6 6 
Q . A eS 5 5 . BEA * Ye 
3 . ese 15 15 ster was ; 9 9 
6 - 12 12 12 12 Coa: ena < t Sa08e 
‘ 9 9 ee 312 3 3 24 6 6 
8 8 8 8 > 3 ae << 1 tae ie 
; o te ER OR 5 2% § § - 6 8.8 
I 228 8 8 3 oy ae eres 7 2 oe 
15 15 6 6 geek ue GS. as 
g . 12 12 12 12 Bee ~ ee ee ts 
3 OF 2 2 2 13 3 3 319 8 8 
Z oe 5 5 36 2.9 117 9 9 
& 9 9 9 9 SF bee Saag 18 Ong 
S 3 . at 6 . 12 6 ‘ 9 3 om 4:6 
3 3 6 a1 g 612 ' 3 9 3 15 C25 
= . 22 8 eu Oe 2 6 2 % 10:3 
ae 8 22 «£9 2 2 6 Gm 8:.. 
. a2 2 ji Gos . bese 216 8.9 
: 4 2 22 . & Be 4°83 9 210 7 5 
3 12 15 6 - 3 615 Nia oe 9 9 ©. 
3 12 6 15 3 .15 6 : : oo) 208 
. 3 ar 6 S . 39 6 6 ; 3939 
$3 6 a1 . § Ge oS 399 3 
; sms S.. £3 ¢ § 5:2 410 410 
3 21 . Se o's 2.5 41010 4 
22 8 ina 9D es 4 5 7 2:20 
8 22 - ae oe ey. & 5 710 2 
22 2 4 38 2 6 2 : 1 oe 
2 22 ~ ae tee 2 6 : 8 14 
; . 21 6 Ss . ta 6 3 9 6 3 
3 6 21 - 3 612 3 9 615 3 
3 6 15 q..%1% S 9 9 oe 
3 15 6 - 3 615 9 9 
* * * * * * 
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TABLE 9 
Parental genotype AB/Ab/aB/ab/ab/ab, representing a set of 4. 
Row divisor: 120 
Composition of gamete 
AB aB ab Ab | AB AB ab ab | AB ABab ab | AB AB AbaB | ABab AB AB | 
AB aB ab Ab | AB AB ab ab | Ab aB aB Ab | ABab AbaB | Ab aBaB Ab | 
AB aB ab Ab | Ab aB aB Ab | Ab aB aB Ab | ab ab aB Ab | aB Abab ab 
I é. ie “a > «188 lem a ae 6 18 18 18 
2 a a 4 41212 4 41212 1212 4 4 AS Rie, 2 
3 ; ; : «hee ; 4 4 220 2 2 4 24 12 12 
4 a 4 41212 4 41212 oo 4 41212 
2 I 321 3 $ 315 15 1 3 Oo @ t-te 2 2 © to es 
a ‘ 2 : 6 6 218 4 4 2 26 10 10 
7 > 32 3 s 2-9 8 S$ $ 35 35 918 7 7 wee tee %. 
8 2 is 1 * 310 a8 ; 1 oe 212 6 6 4 24 12 12 
9 «8 - +1414. 44 31755 $2799 
10 I gat 3 3 315 15 ra 60 ; fae tea 
II es 4 41212 4 12 12 red ~ - Wee 2 
12 ee oo - « 2639 . 7a: ae - 12 6 6 + 24 12 12 
13 >. “h1% 2 214 14 2 21414 8 8 8 8 Pungiew te 
14 Se * +34 = 2 2m +10 44 6 34 10 10 
15 « 416% 2 214 14 2 21414 arene 8 8 8 8 
16 I 32l 3 '$2 46% 3 31515 i, ata 913 7 7 
g 1 .. «ae os 4 412 12 4 41212 4 41212 a ete: 
31 - i «em y & 6 1139 5 5 331 9 9 
B 19 6 14 14 4 4 212 6 6 432 8 8 
20 6 18 18 $999 s97 9% 
& 
Ee) 2. 2 .30 6 Gos Ge @ 3931 5 9 715 
& 22 om. - 2 6 30 « & 6.6 a 9.2 9 5 915 7 
8 23 . 10 8 2 .28 4 4 «460 2 210 6 2 210 6 14 
= 24 210 . - 2 428 ~ °“@ Ge 210 2 6 21014 6 
25 - 10 2 2 .28 8 4.16% re 513 315 
26 210. - 2 8 8 - 4 410 $9 315 51315 3 
27 I 3.2t 3 3. 935 $ 315 9 Or we gen Oe Ae a Be 
28 I 32I 3 3 115 9 I 3 915 5 ie ae > 7199 
29 ‘So 2 . 30 6 G-. @ 6  & C4 412 412 
go «a . - 2 6 90 - 666 46 4 6 41212 4 
31 » Gea < $e oe 6 .12 8 26 6 2 412 412 
32 @6. s «(ee - 6 8 12 262 6 41212 4 
33 ‘30 2 ¢.@Bu 4 «to & a oe ae 315 513 
34 9:36 « - 2 428 - 4 810 3 S.8-9 31513 5 
35} . 10 2 2 .28 8 4.16% 2 610 2 614 210 
36 210 . - 2 8 8 . 4 o%6 2 6 210 61410 2 
37 7! oer 2 .go0 6 6 . 6 'G $$.909 935-5 9 
38 a! ee - 2 6 30 ~ & G6 1339 715 9 5 
39\;| ' 321 3 g 135 9 r-$ 9% 427 9% + (Ss eee 
40J; I 32t 3 + 3 935 $ F353 9 7 733 9 
* * * *€ * * * * *£ * * * * * * ¢ 
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| AB aB ab Ab | AB AB ab ab | Ab aB aB Ab | AB ab Ab 
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Mode of formation 
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TABLE 10 


Parental genotype AB/AB/AB/ab/ab/ab, representing a set of 2. 
Row divisor: 20 


Composition of gamete 
AB aB ab Ab | AB AB ab ab | AB AB ab ab | AB AB Ab aB | AB ab AB AB 
aB \A 


b aBaB Ab 
AB aB ab Ab | Ab aB aB Ab | Ab aB aB Ab | ab ab aB Ab aB Abab ab 





I I ‘ 9 9 I le ot 
4 4 | ‘ 6 6 a xg ay 
‘ ‘ Te oy he 4 4 gig 
4 4 } . Py is 6 « é é 
I I | s £39 a4 3 3 t-3 2 
, ; ‘ . 3 32 es 3 2244 
I I . ss 3 toe eae hee: 
‘ . . it $s. 6 6 : 
i: e.4°% . i oe sa 4% 
I I 3 3 3 ‘ i or ‘ 
‘ : 5 
+ tne I ‘ 3 3 «6S 
2 2 a as S$ ssn 2222 a 
ns 2:8 E ea 24% 22 : 
2 2 4 33 : * £2 7s 2222 
bE <2 5 3 3 3 . gs ‘ 
+ ‘ ‘ 6 6 ee ate 
‘ ; E} 2° ae ee as sa a's 
, ' 2 ee 22 4422 
ff. 2 ‘ ; a a 
I I $ - 8 i $$ 
I I - § 3 ee ee. 3S 
. . eo. #5. e236 x > in ‘ 4 4 
: ° si. ee ; 1-8 eee eG ‘ 4 4 
. : 2 x ee a S25. % a eR 4 
‘ ‘ : 2e ss : 4 2 > eae 
I I ce 3 ee ‘ ee as Lt 
I I 3 Mee ; 2 3 > ies + 2. ae 
I I 3 — ms oe 333 3 
I I _— 3 iss 2 Mis $8334 
cee + ; rea fe5e a se 26 
bs a 3 : 23% 2 S94 
: 2 2 2-3 92 2 2 2° @ te 
; : Sse. ; 2-34 22 > 3 3 
: : ot o 2 >°.8 _ 44 4 
: , ' £2 +. 2 2 , a 44 
I I S 2 2% re ; 3 3 3 3 
I I { 3 se . . 3 3 a 
I I o > ya8 . ae 3 3 3 e 4 
I 1 - ¥ - gB 3 oS 3 3 
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TABLE 11 
Parental genotype AB/AB/Ab/aB/ab/ab, representing a set of 2. 
Row divisor : 180 
Composition of gamete 
AB aB ab Ab | AB AB ab ab | AB ABab ab | AB AB AbaB | ABab AB AB 
AB aB ab Ab | AB AB ab ab | Ab aB aB Ab | ABab AbaB | Ab aBaB Ab 
AB aB ab Ab | Ab aB aB Ab | Ab aB aB Ab | ab ab aB Ab | aB Abab ab 
I - ete 999 9 i ae ie |) Saree 18 18 36 36 
aim .t%. Wm2Iz2@w2 wWI2I212 224 6 6 oo ee 
3 sane a 4 4242 444 4 1919 4 4 22 22 28 28 
At 38) 2 MH « I2 12 12 12 12 12 12 12 ea es 6 6 24 24 
5 9 6 g 6 | 15 15 15 15 6 6 6 6) at ar 1212 4. Oe 
6 ee a 4 4 4 4444 44111 18 18 30 30 
i 9 6 g 6 6 6 6 15151515 QI QI 12 12 & Ge des 
Sas ae es 9999 she ta cxeet Od IZ 121212 24 24 24 24 
9 ore aie 446464 44 4 4 116161010 25 25 22 22 
10 9 6 g 6 | 15 15 15 15 6 6 6 6 i Ae oe QI 2I 12 12 
1a oe I2I2I212 12121212 « - » & “See Ore 
12 Se ama og ra aes 9 9 9 I2IZIZIQ 24 24 24 24 
13 4:10 410 | If Sf 3P 8 Il Il Ir I! 16 16 16 16 a ah 
14 Be UE 4 4.4 4 444 4 9 91212 32 32 20 20 
I 4:10 410 | TER E802 ER EE IP 81 ae Sh ae 16 16 16 16 
I 9 6 g 6 6 6 6 6 15 15 15 15 a ee QI 21 12 12 
oe i7 |) '2 «32 . 12121212 12 12 12 12 6 6 24 24 PR EF 
~ 18 Paes 444 4 444 4 9 91515 23 23 26 26 
S 19 I I 444% 444 4 +j.%&!mt31414 go go 18 18 
E 20 9999 an & 9 91818 27 27 18 18 
a 
% ar} 3 9 - |at gat g 12 $12 9 1212 § § 121221 2 
aan @ -. -@ « 321 3 2! $12 $912 1212 39 3 1212 21 at 
3 23} 2222/16 116 1 1 416 4 1010 7 7 + 1010 22 22 
2 24 2222 116 1 16 416 416 1010 7 7 10 10 22 22 
25 962 21136 416 4 6 2196 2 nid 5 5 181815 15 
26 222 2 416 4 16 116 116 wd 5 5 1818 15 15 
27 9 6 g 6 615 615 15 615 6 a ae 12 12 QI 21 
28 9 6 9 6/15 615 6 615 615 », Ee og 12 12 QI QI 
29 §$-+«+3-«- |] 32 BB Ww gas sf 999 9 56% 15:55 
30 ree ee 321 321 3 12 12 999 9 1515 15 15 
31} © @. g& |e git St aot 8 8 8 8 16 16 16 16 
32 a Gereae g 12 312 321 3 21 8 8 8 8 16 16 16 16 
33\ 22 2°) 06 £363 16 216 & 5 511 Ir 15 15 18 18 
34 2222 116 1 16 416 4 16 5 5 1111 15 15 18 18 
35\ $222 /16 416% 16 416 7 71010 22 22 10 10 
36 222 2 416 4 16 316 1 36 7 71010 22 22 10 10 
37 G6. 3 . )S Oat. 9 te 918 2 3 31212 Qt at i212 
38 Pe ae 321 3 21 312 3 3 31212 QI aI 12 12 
39 9 6 9 6/15 615 6 615 615 12 12 21 aI + « % 
40 9 6 g 6 615 615 15 615 6 1212 21 21 
* * * * * * * * * * 
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TABLE 12 
Parental genotype Ab/Ab/aB/ab/ab/ab, representing a set of 4. An interchange of loci 
accompanied by an interchange of paired conjugate modes generates a conjugate set of 4. 
Row divisor : 60 


Composition of gamete 











| AB aB ab Ab | AB AB ab ab | AB AB ab ab | AB AB Ab aB | AB ab AB AB 
| AB aB ab Ab | 4B AB ab ab | Ab aB aB Ab | ABab AbaB | Ab aBaB aB 
| AB aB ab Ab | Ab aB aB Ab | Ab aB aB Ab | ab ab aB Ab | aB Abab ab 
1 ae ‘ g 18 or ae 9 » « 
2 12 4 : 6 12 s « €©m 4 4 ete, 
3 “aaa ‘ 7 18 . ¥ a ae - 6. ‘@ 
4 12 4 : 6 12 6 12 eis Aa . ts 
5 15 3 . 3 9 I 6 13 3 3 4 co . 
6 me i+ 3a 5 Y 2 Bs - 14 4 
7 Bs |! 6 13 ee  aaet “i . 
8 eG 2 ae . * ie fs + 2 
9 5: 7 18 7 - 5 2 - 13 3 
10 15 3 3 9 I 6 13 >t + #4 
II 12 4 6 12 6 12 . & ‘ 4 4 
12 o 4 3 18 . - 6 3 s'@ . 6 
13 18 2 I 3 10 I 3 10 so 6% a4 @ 
14 Pie ws , 7 18 ~ 9 oe oO . 
15 (0° 2 I 3 10 I 3 10 oe *$ 3's 
16 ~ Bs I 6 13 3 9 a ee 3 3 4 
£17 12 4 6 12 6 12 4 4 ee 
3 18 e ea  & - 97 2 II 5 
S 19 — oe 7 . @ 10 4 
= 20 4 - 9 9 a 9 3 
= 
oe 
6 = 9 3 co eke ® I . 2% oo 
© 22 Os “ue 3 10 i 6 8. 
3 23 12 2 j 12 10 2 : = a /, oe 6 
A 24 SB 4 ‘ 3 24 p 1 8 ew ® e686 @ 0: 
25} 12 2 ; 12 6 2 10 | es ae 
26 8 , I 24 ; 3 8 es ee Y 2 e3 
27\ 15 3 I 6 9 3 13 + we o 4 3% 
28S 15 3 oY 2s I 6 9 a oe 4 i.) Ge 
29 9 3 -15 9 I 2: es - 4: 4 
30 6 - § % 3 10 a aoe | ee a 
31\ 15 1 > ae a 3 ~ = @.F, Ss . 6 
32! 10 —s 1 6 a 2442 
33 12 2 - 12 10 2 6 ae eS % : eae 
34 8 : 2 . 1 8 if. @ ok 2:8 
35 | 12 2 < 32 2 . 10 hee ay & i i > 56 
36 8 . me 3 8 tore % 4 6 2 
37 9 3 - 15 9 I - 9 etme ae! x" 3: 7 
38 6 3 24 3 10 , 2 3 6 1 
39 15 3 - 313 I 69 4383. , 
40, 15 3 I 6 9g S45 8 2s 2 














ae 
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TABLE 13 


Parental genotype AB/Ab/ab/ab/ab/ab, representing a set of 4. An interchange of loci 


accompanied by an interchange of paired conjugate modes generates a conjugate set of 4. 
Row divisor : 30 


Composition of gamete 











AB aB ab Ab | AB AB ab ab | AB ABab ab | AB AB AbaB | ABab AB AB 
AB aB ab Ab | AB AB ab ab | Ab aB aB Ab | ABab AbaB | Ab aBaB Ab 
AB aB ab Ab | Ab aB aB Ab | Ab aB aB Ab | ab ab aB Ab |aB Abab ab 
| 
I 6 ia . g ‘ ~ * oe. , « 6 
2 12 I <a I ~ _ 4 4 ‘ qh 
3 4 F 29 ‘ ae , ie o Re 
4 12 I ae I a es a’ 
5 gq 1 I 3 7 ‘ oe asi ge, ta oe 
6 o . E 4 9 ae a Ca 4 2 
7 ® 1 ‘ - I 3 7 ie a oe aes ; 
8 OP : _ 2 : - - 3 2 « G 4 
9 4 ; 29 , a a a om 3 
10 9 I I 3 7 ; - § oe a 4 1 I 
II 12. I ae I in ee a 4 4 : 
12 a : 6 9 : - 6 . 4 - 6 . 
13 6 2 ‘ 3 6 3 6 e's 2 4 ee - 
14 4 . ‘ 2 9 - @ Pe st . @ 2 
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TABLE 14 





Parental genotype AB/Ab/Ab/ab/ab/ab, representing a set of 4. An interchange of loci 
accompanied by an interchange of paired conjugate modes generates a conjugate set of 4. 


Row divisor : 60 


Composition of gamete 














AB aB ab Ab | AB AB ab ab | AB AB abab | AB AB Ab aB | ABab AB AB 
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Composition of gamete 


Parental genotype AB/Ab/Ab/aB/ab/ab, representing a set of 2. 





An interchange of loci 
accompanied by an interchange of paired conjugate modes generates a conjugate set of 2. 
Row divisor: 180 





© ON OUP OH 


Mode of formation 
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When we take account of the gametic matrices for the five conjugate 
sets of genotypes (which can be developed from the matrices of 
tables 12-16 by an interchange of loci and of pairs of conjugate modes) 
there are 204 such columns. The rank of this 40204 matrix is 
equal to the number of independent linear combinations of the forty 
frequencies which become available for estimation on classifying the 
gametic output of hexasomic organisms with two alleles at each of 
two linked loci. 

Among the 147 columns marked out in the twelve gametic matrices, 
there are only twenty-four that are independent. Moreover, the 
gametic matrices for the five conjugate sets of genotypes contribute 


TABLE 17 


Linear combinations capable of estimation 


Column reference numbers 
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no further independent combination. It happens that estimates of 
these twenty-four linear combinations can be formed from the gametic 
output of genotypes belonging only to self-conjugate isomorphic sets. 
The following example illustrates how this can be done. The self- 
conjugate set of genotypes which includes bisimplex coupling (AB) (ab) Bs 
is competent to estimate the eleven linear combinations shown in 
the first eleven columns of table 17. The self-conjugate set of four 
genotypes which includes bisimplex repulsion (Ad) (aB)(ab),, also leads 
to estimates of eleven independent linear combinations, but of these, 
nine are held in common with those of the former set. The two 
independent combinations are shown in columns 12 and 13 of table 17. 
The self-conjugate set of four genotypes which includes biduplex 
coupling (AB),(ab),, contributes the seven independent combinations 
shown in columns 14-20, whilst the self-conjugate set of four genotypes 
which includes biduplex repulsion (Ab) ,(aB),(ab)., adds three further 
columns (21-23). The self-conjugate set of four genotypes which 
includes biduplex sub-coupling (AB)(Ad)(aB)(ab)5,* yields no linear 
combination that is independent of the first twenty-three columns 
of table 17. The remaining linear combination, shown in column 24, 
is found in the gametic matrix for the self-conjugate pair of genotypes 
one of which is bitriplex coupling (AB),(ab),;. The progression in 
rank of the combined gametic matrices is shown in table 18. 











TABLE 18 
Parental genot Rank of gametic Rank of combined 
8 init matrix gametic matrix 

(AB)(ab),_ . ; ; : Il ion 
(Ab) (aB) (ab) 4 . : : II 13 

)2(ab) 4 . ° . . 16 20 
(Ab) .(aB) (ab). , : 16 23 | 
(AB),(ab),_ - : M : 10 24 | 














The 40X16 deficiency matrix is shown in table 1g. It is seen 
that the frequencies of only three modes of gamete formation, namely 
12, 31 and 32 can be estimated singly. This deficiency matrix merits 
further consideration as it must contain the essentials of a system of 
classification of the modes of gamete formation. 


5. SUMMARY 


Gametic matrices are given for twelve isomorphic sets of hexasomic 
genotypes with two alleles at each of two linked loci. It is found that 
from a classification at two linked loci of the gametic output from 
digenic hexasomic organisms, we can estimate only twenty-four 


* There are four neutral genotypes. These are duplo-triplex and comprise a conjugate 
pair of isomorphic sets each of two genotypes. A typical genotype is (AB) (Abd) (4B) (ab) »- 
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independent linear combinations of the frequencies of the forty modes 
of gamete formation. The deficiency of information is therefore much 
greater than with two tetrasomic loci where, from doubly digenic 
genotypes, we can estimate ten independent linear combinations of 
the frequencies of the eleven modes of gamete formation. 


TABLE 19 


Deficiency matrix 





16 —I ; : : : ._-!I ‘ : ; P ; ; ‘ .--! 
—I ._-!I ._—!I 
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Mode of gamete formation 
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1, INTRODUCTION 


THE gametic output from polysomic organisms of various genotypes 
can be expressed in terms of a single parameter, namely, the frequency 
of double reduction. This is the frequency with which a gamete 
receives two alleles derived from the same chromosome at meiosis. 

The term “ double reduction ”’ was first introduced by Darlington 
(1929) and defined as a case “‘ where both chromatids of one chromo- 
some pass into the same nucleus at the second telophase”. This 
offered an explanation for the appearance of one recessive in a 
progeny of 137 from a reciprocal cross between a triplex and nulliplex 
Dahlia obtained by Lawrence (1929). Similar results had been 
obtained by Blakeslee, Belling and Farnham (1923), who got 2 per 
cent. recessives from a triplex Datura when backcrossed to the 
recessive, and 1 per cent. when crossed to the simplex genotype. 
These early workers, however, expected to find an abnormal number 
of chromosomes in the recessives. Blakeslee et al. examined their 
recessives but found the numbers to be normal. 

Haldane (1930) pointed out that though the passing of two identical 
chromatids to the same gamete in a diploid organism implied an 
abnormal chromosome complement in the resulting zygote, this was 
not the case in polyploids. He calculated the ratios to be expected if 
the eight chromatids of a tetraploid were distributed at random into 
the gametes. In the mating between a triplex and a nulliplex, 
one recessive in 28 could be expected on this hypothesis. Darlington 
(1931), however, showed that cytologically there were a number of 
restrictions to Haldane’s theory, and concluded that intermediate 
conditions between random segregation of chromatids and diploid 
segregation were to be expected, and had in fact been observed by 
Lawrence. 

In four crosses involving one triplex with four different nulliplex, 
Lawrence (1935) observed double reduction at the Y locus in tetraploid 
Dahlia. He obtained 7 recessives out of 385 plants and explained 
this as “‘ probably due to chromatid segregation”. At none of the 
three other loci studied did double reduction appear. It should be 
noted, however, that it was never certain that a triplex was involved 
in the matings. 

An important contribution was made by Mather (1936) when he 
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pointed out that “the amount of double reduction is a function of 
the cross-over distance between the locus and the spindle fibre,” and, 
therefore, will vary from locus to locus. This has since been demon- 
strated in Lythrum salicaria by Fisher (1949), who gave preliminary 
estimates of the frequency of double reduction at the purple/pink 
locus as 6-5 and at the short/non-short locus as 2:5. These estimates 
were obtained incidentally in an attempt to ascertain the frequencies 
of the eleven modes of gamete formation now known to be characteristic 
of linkage in tetrasomic organisms. They are based on the proportion 
of duplex plants obtained from simplex parents backcrossed to the 
recessive long-styled pink. This method involves making a second 
backcross to determine the genotype—whether duplex or simplex— 
of the dominant offspring. 

An estimate based on the proportion of recessives obtained on 
backcrossing the triplex genotype is more direct, requiring no second 
backcross and can, in consequence, be based on large numbers with 
greater facility. ‘This latter method has been used to estimate the 
frequencies of double reduction both in pollen and ovules at the 
mid locus in the work described here. It should be noted that for 
an estimate at the short locus this method cannot at present be used, 
since a triplex short could only be obtained by illegitimate fertilisation. 
All attempts to obtain such a plant have so far failed. 


2. DOUBLE REDUCTION IN OVULE FORMATION 


When Fisher and Mather (1943) established the existence of 
polysomic inheritance in L. salicaria, they did so by backcrossing to 
the long-styled form, 16 Shorts and 12 Mids obtained from a cross 
between a Short and a Mid, both of which had one and the same 
dominant mid allele. The resulting progenies showed that one 
Short was nulliplex for mid, 11 were simplex and 4 duplex ; 8 of the 
Mids were simplex, 3 duplex and 1 plant, M4o-2F, gave no Longs 
out of 41 non-shorts. This plant was tested the next year with a 
progeny of 1034 (Fisher, 1944). It gave 22 Longs out of a total of 
1075 in two years, and must, therefore, be a triplex and have arisen 
by double reduction. This first evaluation of the frequency gave an 
estimate of 8-19 +1-73 per cent. 

In an attempt to obtain a Short triplex for mid, two of the Shorts 
that were found to be duplex for mid were crosszd to a simplex Mid 
(Ss; Mm, = s, Mm), and twenty-five Shorts obtained from this 
cross were tested in 1945. The twenty-five progenies are given in 
table 1, from which it will be seen that five plants triplex for mid were 
identified. Fifteen of these plants were again backcrossed, this time 
to pale pinks to provide triplex pales and data on the linkage of the 
Short and purple loci. The test progenies grown in 1946 are given 
in table 2, in which the Long, Mid and Short segregations alone are 
given. There is evidence that some stray pollen was present from the 
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2 Mids that appeared in the progeny of the nulliplex plant. The data 
from the four triplex plants will nevertheless be included in our estimate 
of double reduction, for, as will be shown later, a small percentage of 
foreign pollen does not alter the estimate appreciably. The fifth 
triplex plant was crossed to the triplex Mid, M4o-2F, in an attempt 
to obtain both Shorts and Mids quadruplex for mid. In 1947 eleven 


























TABLE 1 
Test progenies from 25 Shorts (1945) 
| | 

Genotype | Plant L M | Ss | Total 

‘a hee Ge | 
Nulliplex . $44°2G15 105 89 | 194 
Simplex ‘ $44°2F1 30 63 105 198 
$44°2F2 39 62 ao 198 
$44°2F3 47 | 47 | 107 201 
S44°2F4 38 60 101 199 
S44°2F5 - I Ss | 9 | 202 
$44°2F6 39 43 114 | 196 
$44°2G1 41 57 102 200 
$44°2G5 55 | 62 84 | 201 
S44°2G8 30 51 123 204 
S44°2G16 57 | 51 88 196 
$44:2G18 36 42 84 162 
| Total . ‘ 464 | 593 | Ir00 | 2157 

| ay ae: 
Duplex 3 S44°2F7 12 94 95 201 
$44°2G3 14 go g6 200 
$44°2G4 8 | gI 102 201 
S$44°2G6 13 85 105 203 
S44°2G7 16 68 | 114 198 
S$44°2G9 21 71 109 | 201 
S44°2G10 18 85 | 98 | 201 
$44°2G14 10 83 | 96 189 
Total . 112 66, —| 815 | 1504 

= et) See a | 

Triplex : S44°2G2 2 | 95 105 202 
$44°2Gr1 I | 87 III 199 
$44°2G12 4 | 89 106 199 
S$44°2G13 a 93 | 100 | 197 
$44°:2G17 I | 88 109 198 
Total . ; 12 | 452 531 995 














Shorts and twelve Mids obtained from this cross were backcrossed to 
Longs and progenies grown in 1948. It will be seen from the square- 
root chart in fig. 1 that five Shorts and six Mids triplex for mid 
were found. 

The possibility that the eleven plants classified as triplex were not 
all that should have been so classified, i.e. that one or more had been 
classified as quadruplex, may be examined as follows. The number 
of Long progeny observed from nearly equal progenies with the same 
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frequency of double reduction, should be distributed nearly in a Poisson 
Series. If m is the expectation, then the frequency of cases in which 
none occur will be e-”, and the frequency of cases with at least one 






































TABLE 2 
Test progenies from 15 Shorts (1946) 

Genotype Plant | L | M S Total 
| re | | 

Nulliplex . S$44-2Gi5, | 48 2 63 | 113 
Simplex ° S$44°2F3 42 44 101 187 
S$44°2G1 32 | 41 93 166 

$44°2G5 35 | 51 95 181 

$44°2G8 29 55 g2 176 
$44°2G16 34 40 108 182 | 
$44:2G18 42 44 101 187 | 
Total . ; 214 275 590 1079 | 
Seen MCILEED ane (Cuesialen 2 hese 82 Mar 
Duplex : $44°2G7 io 68 95 72) (| 

$44°2G9 13 | 71 108 192 
$44°2G10 9 72 108 189 
S$44:2G14 6 | 60 | 69 1395 | 
Total . ‘ 37 | 271 380 688 

Triplex ; S$44°2G11 3 47 46 96 
$44°2G12 2 64 110 176 

$44°2G13 3 | 80 98 181 

S44°2G17 ~ | 79 70 149 

Total . ‘ 8 | 270 324 502 














will be 1—e-”. Among these latter the average number of Long 
plants observed must be m/(1—e-”). Actually the average number 
observed among the 11 progenies classified as from triplex seed parents 








TABLE 3 
Double reduction in ovule formation 

l l l 

Year grown | Progenies | L M Total 
a. eee S| a See eee 5 Sees ees 

1942 and 1943 .| Mg4o-2F | 22 1053 1075 
1945. : of 5 Shorts | 12 452 464 
1946. ‘ .| 4Shorts | 8 270 278 
satel | f 5 Shorts 16 418 434 
=? : : 6 Mids | 15 458 | 473 
1946-1951 ‘ ‘ S$45°2 5 267 | 272 
| Total . ; 78 2918 | 2996 

| | 











is 31/11 or 2-81. This corresponds with an estimated m of 2-61114 
for which e~” is 007345. The number of triplex progenies correspond- 
ing with 11 observed is, therefore, 11/0-92655, or 11°872. If, therefore, 
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we add 7/8 of the average number of plants in a progeny, ¢.g. about 
70, to the number of plants not Long observed in the progenies, 
i.e. if we take an estimate based on 31 out of 977 instead of 31 out 
of 907, we shall have made a reasonable allowance for the possibility 
that a triplex progeny has been excluded. The correction arrived 
at is so trifling, that more elaborate calculations are not necessary. 
Another triplex, S45-2, was identified from a repetition in 1944 
of the cross, Ss; M,m, = s, Mmsg, made the previous year. Progenies 
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from this plant have been grown repeatedly since then. The results 
are shown in the fifth line of table 3, in which all the data bearing 
on ovule formation are summarised. The percentage of Longs in 
the progenies—2-6035 per cent. over all—represents a quarter of a, 
the frequency of double reduction, in ovule formation, which is, 
therefore, estimated at 10-414 per cent. If we include the allowance 
made on the 1948 progenies, the percentage of Longs is 2-544, giving 
an estimate for a = 10-176 percent. Both the crude and the corrected 
data are homogeneous. In the first case yx? = 4-393 and in the 
second x? = 3-304 for 4 degrees of freedom. 
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3. SHORTS AND MIDS QUADRUPLEX FOR MID 


Amongst the plants tested in 1948, four Shorts and five Mids 
gave no Longs in their progenies (fig. 1). Four of these plants, two 
Shorts and two Mids, were tested again in 1949 with larger progenies. 
The combined data from both years are given in table 4. 











TABLE 4 
Test progenies from 2 Shorts and 2 Mids quadruplex for Mid 
Plant mee ° M | S Total 
- — —| H 
$4716. 7a os 552 575 1127 
S847:23 . : ae | 57! 582 1153 
M478 . : I 546 cas 547 
M47"12 . ; a 507 2 569 











The presence of two Shorts in the progeny of the Mid, M47:12, 
must be due to stray pollen at the crossing plot or to contamination 
in the handling of the seeds, which are very small. The Long observed 
in the progeny of the other Mid must also be a contaminant, since, 
if this Mid were a triplex 15-65 Longs would be expected from 547 
plants. It was, therefore, concluded that all four plants are quadruplex 
for mid. 


4. DOUBLE REDUCTION IN POLLEN FORMATION 


To obtain an estimate of a in pollen formation it was necessary 
to use the triplex Mid as pollen parent and the recessive Long as seed 
parent. As demonstrated by Fisher and Mather (1943), the advantages 
of the technique of open pollination, which we have used for all 
tests in this investigation, far outweigh the disadvantage of the risk 
of contamination by foreign pollen, and they found this to be less 
than 1 per cent. in their progenies. It should be pointed out, however, 
that the risk of error in the estimate of double reduction due to the 
same percentage of foreign pollen is about fifty times greater when the 
seed parent is a Long than when it is a Mid or Short. This becomes 
apparent when it is realised that all the ovules from a Long will, if 
fertilised by pollen carrying neither mid nor short, increase the 
proportion of long-styled individuals in the progeny. Since calculation 
shows that 42 per cent. of such pollen exists in a theoretical population 
(Fisher, 1941), even 5 per cent. stray pollen at the crossing plot would 
almost double the percentage of Longs observed in the progeny. On 
the other hand, when the triplex Mid is the seed parent, only those 
ovules not carrying mid, that is about 2 per cent., can if fertilised 
by pollen carrying mid, result in a deficiency of Longs in the progeny. 
If 5 per cent. foreign pollen were present, the deficiency would be less 
than 0-05 per cent., a difference that could make no serious difference 
to the estimate of a. 

















DOUBLE REDUCTION IN LYTHRUM SALICARIA 2g 


With this very real source of error in mind, great care was taken 
in 1948 to find a plot far from any known Lythrum populations 
for the pollination of a long-styled plant by the triplex Mid, M4o-2F. 
To make plenty of the right pollen available, the Mid was divided 
and planted in a ring round the Long. Seed taken from the Long 
was grown first in 1949 and again in 1951. The data are given in 
table 5. 




















TABLE 5 
Double reduction in pollen formation 
Year L SS | Total 
on ae! Ei # 
19949 + 10 374 | 384 
195! ‘ ; 21 943 eee 964 
Total . J 31 1317 ‘saa | 1348 














The percentage of Longs—2-300 per cent.—is slightly less than 
that observed in ovule formation. This, and the fact that no Shorts 
were observed in either year, is convincing evidence that foreign 
pollen was absent at the crossing plot. The estimate of a is 9-199 per 
cent. The data for both years are homogeneous—y*,,. = 02225. 

The data from all years, 1943-1951, from both ovules and pollen 
combined are homogeneous, x? = 5-110 for 6 degrees of freedom, 
indicating that seasonal variations do not affect the frequency of 
double reduction. Taking the two estimates together we get 
2°509+°237 per cent. Longs and a = 10-:037+-949 per cent., and 
x? = 0-351 for 1 degree of freedom, on the crude data. If we include 
a twelfth progeny from the 1948 families, the percentage of Longs 
is 2°469-+-+234 per cent., giving a = 9-878+-936 per cent., and 
x? = 0-232 for 1 degree of freedom. Our estimate is only slightly 
higher than that forecast by Fisher and Mather in 1943. 

I am much indebted to Professor Sir Ronald Fisher, not only for his constant 


guidance and help in every aspect of the work, but for supplying the material 
and suggesting the investigation. 


5. SUMMARY 


About the time when it was realised that the gametic output from 
polysomic organisms of various genotypes could be expressed in terms 
of a single parameter, viz. the frequency of double reduction, a, an 
attempt was first made to estimate this frequency at the mid locus in 
Lythrum salicaria. When Fisher and Mather (1943) established poly- 
somic inheritance in L. salicaria, they obtained a triplex mid that could 
only have arisen by double reduction. This plant was backcrossed 
to the recessive long-styled form and in two seasons gave 22 Longs 
out of 1075. The proportion of Longs represents a quarter of a in 
ovule formation, which was then estimated at 9-18 per cent. Similar 
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data has been obtained in subsequent years. In 1948 triplex and 
quadruplex plants were distinguished only by their performance in 
producing or not producing Longs. Allowance has been made for 
the possibility that one or more triplex had been classified as a 
quadruplex by adding 7/8 of the average number of plants in a 
progeny, 70, to the Mids observed. The corrected total of 78 out 
of 3066 plants were obtained in five years. The proportion of Longs 
is 2°544 per cent. giving a = 10'176 per cent. 

To obtain an estimate of a in pollen formation the reciprocal cross 
was made in 1948. In such a cross the risk of error from contamination 
by foreign pollen is about 50 times greater than when the triplex 
plant is the seed parent. Since the open pollination technique was 
used throughout, special care was taken in selecting a site for the 
pollination plot for this cross. In two years 31 Longs were obtained 
out of 1348 plants ; no Shorts were observed, thus providing evidence 
that foreign pollen was absent. The percentage of Longs is 2-300 per 
cent., giving an estimate of a = 9-199 per cent. in pollen formation, 
very close to the value for ovules. Combining the data from ovules 
and pollen the proportion of Longs is 2-469+-°234 per cent., and 
a = 9°878+-936 per cent. The data obtained from 1942 to 195! 
are homogeneous, x24) = 5*110, indicating that a is not affected by 
seasonal differences. 

Of the Shorts and Mids likely to be quadruplex for mid, two of 
each type have been tested with adequate numbers and their genotype 
verified. All four were morphologically normal showing the com- 
pleteness of dominance at the mid locus. 
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QUANTITATIVE INHERITANCE. Papers read at a colloquium held at the Institute of 
Animal Genetics, Edinburgh University, under the auspices of the Agricultural Research 
Council, 4th to 6th April 1950. Her Majesty’s Stationery Office, London, 1952. 
Price 20s. net. 


It was enterprising of the Agricultural Research Council to organise 
and to publish the contributions to this symposium on Quantitative Inheritance. 
Nearly all the contributors stress the importance of quantitative characters 
in evolution and in animal and plant improvement, and this represents a 
real advance from the attitude of the first generation of geneticists, about 
1910, though it is no novelty for the present generation. 

Like most such occasions, this discussion suffered to some extent from 
the reiteration by the older hands of opinions of which they have ceased 
to be critical, and by the naive misapprehensions of the newcomers, 
Nevertheless, some material of interest has been gathered. 

Attention may be drawn to the short paper by L. L. Cavalli, “ An 
analysis of linkage in quantitative inheritance” using the methods and 
material of K. Mather. Instructive also is Mather’s gentle but effective 
reproof of the extravagances of a purely negative contribution by B. Woolf. 

The paper, however, from which most can be learnt, and which is the 
chief justification for the publication, is that of Mather and Vines on “‘ The 
inheritance of height and flowering time in a cross of Nicotiana rustica.” 
This sets out fully and usefully the analysis up to the fourth generation, 
with a valuable discussion of those factors, such as scaling, which help to 
make the analysis more accurate. R. A. F. 


HEREDITY IN UTERINE CANCER. By Douglas P. Murphy. Harvard University Press 
(London : Geoffrey Cumberlege). 1952. Pp. xi+128. 16s. 


Although geneticists on the whole incline to the view that a predisposition 
to cancer is inherited, the opinions held by many authorities who have 
studied the subject are somewhat diverse. With the possible exception 
of the work of Brobeck, the existing literature seems to lack any really 
satisfactorily comprehensive data on the familial occurrence of cancer, 
involving in particular an adequate series of controls. The present book 
describes the methods and results of an enquiry whose purpose was to 
make good this deficiency. It is concerned with the occurrence of cancer 
among the female relatives of 201 women suffering from cancer of the 
uterine cervix and of 215 controls. In the first section of the book the 
literature is reviewed and discussed, and emphasis is laid on the need for a 
new survey. The second section deals with the materials and methods 
employed, and is a model of clarity and precision. Full details are given 
of the purpose of the investigation, the use made of a pilot study, the selection 
of both the cancer and control probands, the methods of obtaining and 
handling the clinical and social information required, personnel employed, 
and so on. The results achieved are elaborated in the third section and 
discussed in the fourth. Certain details of the administrative background 
and the questionnaires used are presented in appendices. 

The general conclusions of the survey are that hereditary factors do 
influence the frequency of occurrence of uterine cancer, though the effect 
is small. Thus the ratio of cases of uterine cancer to all cases of cancer 
was 26 per cent. in the female relatives of cancer probands, compared 
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with about 15 per cent. in the relatives of the control probands. Again, 
3°2 per cent. of the mothers and aunts of cancer probands were recorded 
as having cancer of the uterus, compared with 1-4 per cent. of the mothers 
and aunts of the controls. The total prevalence of cancer at some site 
did not differ between the cancer and control families. Finally, excluding 
the proband-mother relationships, it is shown that g per cent. of the mothers 
of non-cancerous daughters had cancer in some site, compared with a 
figure of 17 per cent. for mothers of daughters with cancer irrespective of 
site ; whereas only g per cent. of the mothers of the cervical cancer probands 
had cancer in some site. This seems to imply that genetic factors other 
than those concerned with uterine cancer affect the occurrence of cancer 
elsewhere. 

The chief criticism which must be made of the treatment is that most 
of the tables presented are simple one-way analyses of material relating 
either to the cancer family groups or to the control family groups. No 
attempt has been made to examine two- or three-way tables. This is 
important as there are, for example, marked differences in the economic 
and educational levels of the cancer probands compared with the controls 
(and perhaps their families as well). Thus some of the apparent differences 
in the prevalence of cancer referred to above might be ascribable to purely 
environmental factors. However, in spite of this qualification, it is true 
to say that Dr Murphy’s book does make an extremely useful contribution 
to the study of the etiology of cancer, and is well worth the serious attention 
of clinicians and geneticists alike. Norman T. J. BarLey. 


STATISTICAL METHOD IN BIOLOGICAL ASSAY. By D. J. Finney. Griffen. 1952. 
xix+661 pp. 68s. 

Frequently it happens, whether from necessity or convenience, that we 
wish to measure the potency of some substance or stimulus, which may be 
chemical, physical, psychological, etc., by means of the effects produced 
in living subjects. This is the field of biological assay, which, at least so 
far as quantitative methods are concerned, is extremely susceptible of 
considerable statistical development. The great advances of recent years 
in both the scope of bioassay itself and in the statistical design and analysis 
of the experimental procedures have been evolved mainly in relation to 
pharmacology. Thus one of the most convenient accounts of the statistical 
methods available was, until the publication of the present book, the third 
chapter (occupying more than a third of the whole volume) by Dr Finney 
in the second edition of Biological Standardization by Burn, Finney and 
Goodwin. However, there are immediate applications to other subjects 
like psychology or physiology, where, using similar experimental methods, 
interest is centred in the living subjects, rather than in the stimuli them- 
selves. Again, methods of analysis may be used elsewhere because of 
similarities in the underlying mathematical model. There is for example 
the application of the modified form of probit analysis, used when there 
is an independent source of natural mortality, to the problem of estimating 
the average age of onset of a chronic disease from its prevalence in different 
age groups. 

In the present book Dr Finney has, in his own words, “‘ attempted to 
provide a comprehensive account of designs and statistical analyses for 
biological assays, both as a textbook for the student of statistics and as 
a work of reference for the practitioner of bioassay.” This aim has been 
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largely fulfilled. Assuming a certain acquaintance with basic statistical 
methods, a wide variety of problems peculiar to assay work is introduced, 
and the appropriate analytical procedures are explained, illustrated and 
discussed. Thus there are chapters dealing with parallel line assays ; 
slope ratio assays; quantal responses ; symmetrical designs; efficiency, 
reliability and sensitivity ; incomplete block designs ; cross-over designs ; 
use of concomitant information; combination of estimates, etc., etc. 
There is also a fairly extensive bibliography and seventeen appendix tables 
giving the distributions of ¢, x, Fisher-Behrens ratio and variance ratio, 
as well as all the usual tables of transformations, working values, ranges, 
and weighting coefficients, for probits, logits, angles and loglogs. A great 
many general statistical techniques are introduced into the text and full 
explanations are given of their application and interpretation. As a result, 
the book is a more or less self-contained handbook on the subject, valuable 
as a work of reference both to the practitioner of bioassay and to the student 
of the primarily statistical aspects. One useful feature is the inclusion of 
four suggested alternative selections of sections from the book designed 
to meet the needs of the professional statistician desiring a short course 
in biological assay, the non-mathematical user of assay, the reader requiring 
a general survey of the function of statistics in the subject, and the reader 
who is mainly interested in work with quantal responses. 

The general excellence of Statistical Method in Biological Assay justifies, 
in the opinion of the reviewer, the somewhat high price. Because of its 
comprehensive treatment and lucid exposition, Dr Finney’s book deserves 
to be widely read and to become a standard work of reference in the field 
of bioassay. Norman T. J. Batey. 
EVOLUTION IN THE GENUS DROSOPHILA. By J. T. Patterson and W. S. Stone. 

York : The Macmillan Company. 1952. Pp. 1-599. $8.50. 

The evidence of evolutionary processes is at once more comprehensive 
and more significant in Drosophila than in any other group of animals or 
plants. It is derived from comparative anatomy, geographical distribu- 
tion and ecology, combined with breeding and chromosome studies of 
many hundreds of species and of races in all parts of the world. This book 
is an account of the evidence accumulated, chiefly since 1910, to which 
the authors and their collaborators have so largely contributed. 

Each of the different aspects of Drosophila evolution is considered in 
some detail and particularly later in relation to the bases of genetic isolation. 
In general the authors keep their eyes fixed on one problem at a time 
and they do not allow themselves to be distracted, by what seem irrelevant 
clues, from their appointed programme. This is a successful policy up to 
a point. But, beyond that point, the desire to solve problems which were 
envisaged at the beginning suppresses the tendency for new problems to 
push their way into the discussion. 

Moreover, we find that, in keeping with this attitude, the text has been 
put together in detail as though from an index of notes or reprints classified 
chiefly under the names of authors. Nearly half the paragraphs begin 
with these names (sometimes three or four of them): the connexion 
between such paragraphs will not be clear to the reader unacquainted 
with the hundreds of interesting personalities mentioned. On a larger 
scale the order of the book is confused in another way. For example there 
are 56 pages of sporadic references to salivary gland chromosomes before 
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the authors take the plunge and tell us how the meaning of these chromo- 
somes was discovered. Even so they do not tell us the part played in this 
by Heitz and Bauer who arrive on the scene a year late. 

Such discontinuities of style and arrangement are not unrelated to the 
theoretical character of the book. ‘‘ Natural Selection” appears (p. 264) 
as a sub-heading between, and on an equal footing with, “‘ Gene 
Homologies ’’ and ‘‘ Genic Balance”. ‘“‘ Genetic system” is equated with 
*“ genome” (p. 112) : whether the same equation holds later (on p. 537) 
is doubtful. But perhaps most serious of all is the authors’ inevitable 
failure to take advantage of the specific opportunities of discovering 
fundamental principles from these very experiments with Drosophila. 

Of the 121 diagrams which illustrate types of chromosome complement 
in the male, about 30 per cent. reveal no difference between X and Y 
chromosomes. In some instances as an afterthought the authors call 
attention to this apparent identity, e.g. ‘“‘ The largest V is the X chromosome ; 
and the Y is also a V, which is indistinguishable from the X ” (p. 147). 
Now evolution in the sex chromosomes is clearly rapid in Drosophila: four 
species are represented as having no Y at all in the male; others have 
two X’s. Certain questions therefore arise. Are the two apparently 
identical sex chromosomes a newly arisen system with a single gene 
difference? Or are there large but equal-sized differential segments in 
X and Y? Or has the female become the heterozygous sex in these races 
or species ? 

To all these questions the salivary gland chromosomes and the distribu- 
tion of heterochromatin might provide answers. But these matters are 
not brought into relation. Indeed, in the enormous index, sex chromosomes, 
or chromosomes : sex, or crossing-over, or linkage, or differential segments, 
do not appear. 

This book, it is to be feared, will not be attractive to the general student : 
at least not so attractive as the importance of its theme suggests or deserves. 
Nevertheless it should fulfil a special purpose. For one who knows what 
he wants to find out from Drosophila, one who is already familiar with its 
genetical history and uses, and with the problems of evolutionary genetics 
in general, for such a one the references, tables, figures, diagrams and 
maps will suggest any number of new enquiries. Cc. D. D. 


GENETICS FOR SCHOOLS. By Prof. K. Mather, F.R.S. Modern Science Memoirs, No. 31. 
London: Published for the Science Masters’ Association by John Murray. Pp. 36. 
1s. 6d. 


A valuable first year’s course in all aspects of genetics: plant and 
animal, cytological and physiological, practical and theoretical. 


THE BIOLOGY OF MENTAL HEALTH AND DISEASE. With 108 Contributors. Foreword 
by Stanley Cobb, M.D. London: Cassell. Pp. i-xvii+654. 75s. 


This book contains 38 articles, two of which, by F. J. Kallmann and 
R. R. Gates, and a discussion of mongolism, are of genetic interest. 
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